Compresores axiales.pdf
- 1. 32 UNIDAD 4 Compresores axiales. Teoría alar. Efectos viscosos. Pérdidas. Performance 1. Introducción Los compresores axiales tienen ciertas ventajas y desventajas con respecto a los compresores centrífugos. Entre las ventajas se pueden citar menor área frontal (importante para usos aeronáuticos) y mayores relaciones de compresión y eficiencias, aunque estas últimas ventajas no son tan grandes si se consideran compresores centrífugos de varias etapas y de diseño moderno. Las principales desventajas del compresor axial son su costo y su relativamente menor robustez, dada la fragilidad de los álabes (comparando con el rotor centrífugo de una sola pieza). El compresor axial consiste en un rotor de forma cilíndrica que gira dentro de una carcasa o estator. El fluído de trabajo circula por el espacio anular entre el rotor y el estator, pasando por hileras de álabes fijos y móviles (Figura 4.1) Figura 4.1: Compresor axial (Kováts) El rotor está generalmente compuesto de discos en cuyas periferias se montan los álabes móviles (Figura 4.2): Figura 4.2: Rotor de compresor axial típico (P&WA) Los álabes, tanto fijos como móviles, de los compresores axiales, son en su mayoría del tipo de reacción. Por lo tanto, se estudian y diseñan en base a la teoría de perfiles alares de la aerodinámica, por lo que se presenta una breve reseña de la misma. 2. Teoría alar 2.1 General El ala de un avión, y los álabes de las turbomáquinas axiales de reacción, presentan una sección característica denominada sección o perfil alar. Existen un gran número de formas para estas secciones, para aplicaciones específicas, cuyas características geométricas y fluidomecánicas se encuentran en varios textos, manuales e informes de ensayo.
- 2. 33 2.2 Perfiles alares La Figura 4.3 indica los parámetros que definen a un perfil alar: Figura 4.3: Perfil alar, nomenclatura (Gannio) El perfil presenta al flujo un borde de ataque redondeado y uno de salida agudo. Las caras convexa y cóncava definen el espesor del perfil y su línea media. La distancia entre las intersecciones de la línea media con el perfil definen la cuerda t. El largo del álabe (envergadura del ala) es b. 2.3 Teoría del ala Cuando este perfil enfrenta una corriente de velocidad C formando un ángulo se genera sobre el perfil una fuerza que puede descomponerse en dos componentes, una perpendicular y otra paralela a la dirección del flujo libre. Estas fuerzas se denominan sustentación y resistencia al avance. La aparición de la sustentación puede explicarse en fluídos ideales por medios puramente matemáticos mediante la teoría de la circulación, mientras que la resistencia al avance es debida a la viscosidad del fluído y no aparece si se consideran fluídos ideales. Existe un ángulo de ataque −α1 tal que la sustentación desaparece (ángulo de sustentación nula) pero la resistencia al avance en fluídos reales nunca se anula totalmente. 2.3.1 Circulación En el estudio de flujos potenciales se trata el caso de un cilindro rotante inmerso en un flujo potencial y se encuentra que se genera una fuerza perpendicular a la dirección del flujo (efecto Magnus). El valor de esta fuerza L está dado por: Γ = c L ρ , donde Γ es la Circulación, dada por la integral de superficie de la velocidad tangencial del flujo sobre la periferia del cilindro. Este caso del cilindro con circulación puede transformarse mediante una transformación conforme en otra forma geométrica, muy similar a un perfil alar que forma un ángulo α∗ con el flujo. Este es el llamado teorema de Kutta-Joukovsky, y en él se demuestra que la circulación en el plano transformado se puede calcular como: * senα π t C = Γ , por unidad de longitud de envergadura. 2.3.2 Sustentación y resistencia al avance La teoría de la circulación nos permite obtener una primera aproximación al valor de la sustentación. Si la presión dinámica del flujo es 2 2 1 C ρ y la superficie sobre la que actúa L es b.t, definimos el coeficiente de sustentación CL como:
- 3. 34 bt C L CL 2 2 1 ρ = , y, con los resultados de la teoría de circulación, * sen 2 α π = L C Dado que los ángulos de ataque son generalmente pequeños podemos poner * 2πα = L C Como se ha dicho, la resistencia al avance, que se indica con D, es debida a la viscosidad del fluído y no hay una fórmula teórica para evaluarla. No obstante, se define el coeficiente de resistencia al avance bt C D CD 2 2 1 ρ = 2.3.3 Coeficientes de sustentación y resistencia. Diagrama polar De lo visto, la gráfica del coeficiente de sustentación en función del ángulo de ataque es una línea recta con pendiente 2π. Esto se verifica experimentalmente pero sólo hasta ciertos valores máximos de ángulo de ataque (Figura 4.4): Figura 4.4: Coeficiente de sustentación experimental (Abbott) Se nota que la sustentación se anula para un ángulo negativo (-α1). El coeficiente de resistencia también puede representarse en función del ángulo de ataque, aunque es más usual representarlo en función del coeficiente de sustentación, gráfica que se conoce como la Polar del perfil (Figura 4.5):
- 4. 35 Figura 4.5: Polar del perfil (Abbott) Se nota que, como se ha dicho, la resistencia al avance no se anula para ningún ángulo de ataque. 2.4 Efectos viscosos En lo que respecta a la sustentación, la viscosidad causa, en primera instancia, un apartamiento de la teoría tal que la pendiente de la curva de C L es menor que 2π. Luego, se produce el efecto de entrada en pérdida (stall) que causa la terminación de la gráfica con el brusco descenso de la sustentación. La resistencia al avance aumenta aproximadamente como el cuadrado de la sustentación hasta que se produce la entrada en pérdida, cuando aumenta bruscamente. 2.4.1 Distribuciones de presiones Tanto la sustentación como la resistencia al avance son fuerzas resultantes de la presión que ejerce el paso del fluído sobre el perfil alar. Debido a la forma del perfil y al ángulo de ataque la presión en la cara cóncava, inferior, es más alta que la del flujo libre, mientras que en la cara convexa, superior, es menor que la del flujo libre. La Figura 4.6 ilustra la composición de la fuerzas de presión: Figura 4.6: Fuerzas de presión sobre el perfil alar (Kermode)
- 5. 36 La Figura 4.7 ilustra la distribución de presiones sobre los perfiles alares de una rueda fija (estator) que actúa como difusor, aumentando la presión estática y manteniendo constante la presión de estagnación: Figura 4.7: Distribución de presiones en una rueda fija (Balje) 2.4.2 Efecto de ángulo de ataque Tanto la distribución como la magnitud de las presiones dependen del ángulo de ataque. La Figura 4.8 ilustra los cambios típicos de la distribución de presiones: Figura 4.8: Distribución de presiones y ángulo de ataque (Kermode)
- 6. 37 2.4.3 Pérdida de sustentación Para ángulos de ataque positivos el gradiente de presión sobre la cara superior es negativo, es decir, la presión aumenta (es menos negativa) al pasar del borde de ataque al de fuga. El flujo sobre la cara superior enfrenta un incremento de presión y se decelera, particularmente cerca de la superficie, lo que provoca el engrosamiento de la capa límite y su eventual separación de la superficie. Esto es lo que se denomina entrada en pérdida, situación en la que en la parte superior del perfil se forman movimientos turbulentos y el flujo principal se separa del perfil. La depresión en la cara superior del perfil desaparece debido a que el flujo ya no sigue la forma aerodinámica del perfil, y la sustentación decae (Figura 4.9): Figura 4.9: Perfil alar en pérdida (Kermode) 3. Compresores multietapa 3.1 Etapas, diagramas de velocidades Los compresores axiales comprenden un número de etapas, cada una de las cuales consiste en una rueda móvil y una fija. En la primera se le entrega al fluído cantidad de movimiento y la segunda actúa como difusor, recobrando presión estática. Los diagramas de velocidades se construyen en la forma usual, componiendo la velocidad periférica U y la velocidad absoluta C1 para formar la velocidad relativa W1 según la relación 1 1 W U C r r r + = . El ángulo de ataque de W1 se asume habitualmente tal que coincida con la tangente a la línea media del perfil alar, como en el caso del compresor centrífugo. La Figura 4.10 ilustra los diagramas de velocidades de una etapa: Figura 4.10: Diagramas de velocidades de una etapa (Shepherd)
- 7. 38 En la figura se pueden hacer varias observaciones: • La velocidad absoluta de salida de la etapa es igual a la de entrada, por lo que la próxima etapa puede ser idéntica. Lo que cambia es la presión de estagnación en la rueda móvil. • Las fuerzas de sustentación y resistencia al avance son normal y paralela respectivamente a la velocidad relativa de entrada, que es la velocidad de flujo desde el punto de vista del perfil. Tanto L como D tienen componentes en las direcciones tangencial y axial. Las componentes tangenciales son la causa del torque que resiste el movimiento del rotor y absorbe la potencia en el eje. Las componentes axiales causan un empuje axial en el rotor (hacia delante), vale decir, el rotor se comporta en este sentido como una hélice. • La velocidad relativa de salida se asume tangencial al perfil, lo que no es exacto, como se discute en el apartado siguiente. Se debe notar que la velocidad periférica U cambia en magnitud desde la base o raíz del álabe al extremo debido al cambio de radio. Asumiendo la velocidad absoluta constante en magnitud y dirección, la velocidad relativa también cambia desde la raíz al extremo tanto en magnitud como en dirección. Esto hace necesario cambiar el ángulo físico del álabe en el borde de ataque para obtener el correcto ángulo de ataque, es decir, el álabe debe ser alabeado en su longitud. El alabeo involucra un número de cuestiones de diseño de difícil tratamiento, entre las que se mencionan: • Cómo generar el alabeo (forjado, maquinado, tallado químico, alabeo mecánico) • Efecto sobre las características físicas (tensiones residuales, corrosión, fatiga) • Efecto sobre las características mecánicas (distribución de presiones, momentos flexores y torsores, fuerzas centrífugas excéntricas, modos de vibración) • Efecto de las temperaturas (creep, cambio en las características mecánicas) 3.2 Grado de reacción Bajo ciertas condiciones se puede encontrar un grado de reacción óptimo para la etapa axial. La siguiente deducción se puede encontrar en el texto de Shepherd. Haciendo referencia a la Figura 4.10, llamamos L R y L S a la sustentación en el rotor y estator respectivamente (no confundir con LR y LA de la figura, que son las componentes tangencial y axial de la sustentación del rotor) y DR y DS a la resistencia al avance en rotor y estator (nuevamente, no confundir con DR y DA). La potencia mecánica está dada por: ) sen cos ( m R m R D L U N α α + = Las pérdidas en el rotor y estator son el producto de la componente de D en la dirección de la velocidad media relativa al álabe y ésta última: rm R R p V D N = , m S S p V D N = , donde hemos supuesto que DR es aproximadamente paralela a la velocidad relativa media V rm y que DS es aproximadamente paralela a la velocidad (absoluta, en el estator) Vm. Formamos la eficiencia de la etapa como ) sen cos ( 1 potencia pérdidas 1 m R m R m S rm R D L U V D V D α α η + + − = − = Para continuar se hacen las siguientes hipótesis: • Que la componente de DR es despreciable comparada con la de L R en el cálculo de la potencia mecánica (ver Polar del perfil, Figura 4.5). • Que las fuerzas tangenciales en el rotor y en el estator son iguales, ya que una proviene de darle rotación al fluído y la otra de eliminarla; luego: m S m R L L β α cos cos =
- 8. 39 • Que ambos álabes tienen la misma relación D/L. Con ayuda del diagrama de velocidades de la Figura 4.10 y lo anterior se llega a: + − = a m rm UV V V L D 2 2 1 η Para una etapa dada el caudal de masa, la velocidad de rotación y la geometría son dadas, por lo que D, L, U y Va son fijos, y sólo varían las velocidades medias. Además, por construcción, U rm V m V r r r + = . Si Va divide a U en dos fracciones aU y (1-a)U, podemos poner Vrm 2 + Vm 2 =2Va 2 +(aU)2 +((1-a)U)2 . El mínimo de esta expresión es a=1/2, con lo que obtenemos la condición que maximiza el rendimiento de la etapa: m rm V V = Esto es, se trata del diagrama simétrico de la Figura 4.11: Figura 4.11: Diagrama simétrico (Shepherd) En el diagrama simétrico se cumple también que 2 2 2 2 + = U V V a rm con lo que la eficiencia resulta: + − = U Va U a V L D 2 1 2 2 1 η Esta última expresión se puede maximizar con respecto a Va/U, obteniéndose 2 1 = η max a U V Con esto el rendimiento de la etapa resulta L D 2 1− = η y el grado de reacción: 2 1 2 2 2 1 2 1 2 2 2 2 2 1 = − + − − = r r r r V V V V V V ρ Resumiendo, con las hipótesis simplificadoras que se han adoptado, el mayor rendimiento corresponde a la etapa con diagrama de velocidades simétrico y grado de reacción 0.5.
- 9. 40 3.3 Desviación Así como en el caso del compresor centrífugo, el flujo que abandona el álabe móvil resulta ligeramente desviado en dirección opuesta al movimiento del rotor. Este fenómeno se conoce como desviación, y puede considerarse causado por la diferencia de presiones que hace que el flujo tienda a pasar de la cara de presión a la de succión, rodeando el borde de fuga. La magnitud del ángulo de desviación puede estimarse en base al ángulo de salida del flujo y la relación cuerda/espesor del perfil, según la gráfica de la Figura 4.12: Figura 4.12: Estimación del ángulo de desviación (Stepanoff) 3.4 Pérdidas Las pérdidas por rozamiento en la rueda fija se pueden estimar con el coeficiente ϕ utilizado para toberas. Las pérdidas en la rueda móvil se estiman en base a un coeficiente ψ que reduce el valor de la velocidad relativa. Este coeficiente de velocidad se determina en base al ángulo total de giro del vector velocidad al pasar por el perfil (Figura 4.13): Figura 4.13: Coeficiente de pérdida en la rueda móvil (Vivier)
- 10. 41 3.5 Diagrama i-s El diagrama de la compresión en una etapa del compresor axial no difiere del ya visto para el compresor radial (Figura 3.16). 3.6 Relación de compresión La relación de compresión que se obtiene en una sola etapa del compresor axial es relativamente baja; en los primeros diseños era del orden de 1.2, por lo que se necesitaba un elevado número de etapas para alcanzar relaciones altas de compresión. Los compresores modernos alcanzan relaciones de compresión por etapa más altas, del orden de 1.4 a 1.6, incluso haciendo uso de los fenómenos de flujo supersónico y aprovechando las ondas de choque para comprimir el fluído. La Figura 4.14 ilustra la relación entre las relaciones de compresión por etapa y total: Figura 4.14: Relaciones de compresión Se aprecia que las mayores relaciones de compresión por etapa permiten reducir sensiblemente el número de ruedas. 3.7 Mapa de compresor Para cada etapa del compresor axial multietapa puede generarse un mapa como un compresor individual, con su línea de bombeo, área de atoramiento y líneas de velocidad y eficiencia constante. Notar sin embargo que las condiciones de entrada de cada etapa son las de salida de la anterior, por lo que los ejes π1 y π2 y las velocidades reducidas π3 no son las mismas. Para generar el mapa del compresor multietapa se deben adicionar los mapas individuales de las etapas con la observación anterior y la condición de caudal de masa constante. Esta es una tarea sumamente dificultosa por la necesidad de iterar el diseño para lograr un adecuado acoplamiento (matching) entre las etapas, a cada velocidad del rotor. El resultado es un mapa esencialmente similar a los ya vistos para compresores centrífugos (Figura 4.15). Se deben notar sin embargo ciertas diferencias: • Como el fluído es comprimido al pasar por las etapas, su volumen se reduce y, si todos los álabes fueran del mismo largo, la velocidad del flujo en la entrada de la primera etapa sería máxima, y se presentaría el atoramiento. Para aminorar este problema los álabes son de distinto largo, aumentando la sección de entrada en las primeras etapas y reduciendo así la velocidad. Aún así, el atoramiento es más marcado en estos compresores que en los centrífugos de una etapa. El alargamiento de los álabes también magnifica los problemas causados por el alabeo.
- 11. 42 Figura 4.15: Mapa de compresor axial multietapa (Stepanoff) • La variación de largo puede hacerse reduciendo el diámetro de la carcasa o aumentando el del rotor (Figura 4.1). La elección depende de varios factores de detalle tales como flujos secundarios y la distribución de velocidades axiales. • Para cada etapa el resto del compresor se comporta como un reservorio de fluído a presión, por lo que cada etapa está sujeta al riesgo de bombeo. Dado que la compresión por etapa es baja y el acoplamiento correcto entre los diagramas de velocidad es muy difícil de asegurar, el riesgo de bombeo es alto y la línea límite de bombeo se especifica más baja que para los compresores centrífugos. • Como resultado el área de operación es más restringida, con márgenes de bombeo y atoramiento muy estrechos. Con respecto al fenómeno de bombeo, se debe notar que en los compresores axiales multietapa, aparte del caso conocido de bombeo del sistema, pueden presentarse fenómenos de bombeo internos. Esto puede deberse, por ejemplo, a pequeñas diferencias (de manufactura) de ángulo entre álabes, depósitos (polvo y aceite) o daños mecánicos (impacto de objetos succionados por el compresor), que causen que uno o varios álabes de una rueda caigan en la zona de bombeo, o bien entren en pérdida de sustentación. Estos fenómenos localizados pueden en ciertos casos propagarse o desplazarse a velocidades varias alrededor de la rueda, lo que se conoce como “rotating surge” y “rotating stall”. Esto causa altas solicitaciones mecánicas alternativas a los álabes, y es sumamente peligroso.