Facultad de Ingeniería - UNLP

ESTUDIO DE COEFICIENTES ROTODINÁMICOS EN TURBINAS TIPO FRANCIS

QUE PRESENTAN FENÓMENOS DE WHIRLING
Autores: Villar Juan I., Scarabino Ana E.

G.F.C.
Grupo de Fluidodinámica Computacional
Facultad de Ingeniería – Departamento de Aeronáutica
Universidad Nacional de La Plata
Calle 116 y 48
(1900) La Plata
Buenos Aires – Argentina
PALABRAS CLAVE
Whirling, Coeficientes Rotodinámicos, Estabilidad de Rotores, Sellos Anulares
Turbulentos, Efecto Lomakin, Sellos Laberínticos.

RESUMEN
El objetivo de este trabajo de investigación es determinar los coeficientes estáticos y
dinámicos de los laberintos del rodete través de una simulación numérica tridimensional en
CFD, basado en los Papers de referencia y en función de los datos de entrada y condiciones
de funcionamiento. Estos enfoques dan la tendencia creciente a que se comiencen a utilizar
herramientas CFD para calcular las fuerzas presentes en el sistema sujeto a diversas
condiciones de operación. Cada modelos permiten en mayor o menor medida estimar la
respuesta dinámica del rotor frente al fenómeno de whirling, debido a las limitaciones y
ventajas que este impone, estabilizantes o no, convirtiéndolos críticos para el diseño.
INDICE
Resumen e Indice ............................................................................................................................... 1
1 Introducción...................................................................................................................................... 2
2 Parte experimental ........................................................................................................................... 3
3 Resultados y Discusión ..................................................................... ¡Error! Marcador no definido.
Conclusiones ....................................................................................................................................... 6
Bibliografía .......................................................................................................................................... 6

1. Modelo estructural - TURBINA DE BELO MONTE

En las siguientes figuras se puede observar el diseño de la Turbina Francis, y en detalle el

diseño de los anillos laberintos en Corona.

Figura 1. Modelo de Turbina tipo Francis Figura 2. Turbina Detalle - Laberinto Lado Corona

Página 1 / 6 172
Segundas Jornadas de Investigación y Transferencia - 2013

1 Introducción

El objetivo de este trabajo de investigación es determinar los coeficientes estáticos y

dinámicos de los laberintos del rodete, en banda y corona, a través de una simulación
numérica tridimensional en CFD, basado en los Papers de referencia disponibles y en
función de los datos de entrada y condiciones de funcionamiento indicados en el presente
documento. El estudio de Coeficientes roto dinámicos en turbinas Tipo Francis que
presentan fenómenos de whirling, deviene del estudio de las fuerzas roto dinámicas a que
se ve sometido el eje de rotación de la turbo máquina como consecuencia del movimiento
excéntrico de dicho eje. Las fuerzas que el fluido imparte al rotor en el plano
perpendicular al eje de rotación son las llamadas fuerzas roto dinámicas o fuerzas radiales
(referidas a una terna polar local) y son debidas al desplazamiento y movimiento del eje de
rotación (consideramos dicho desplazamiento lo suficientemente pequeño como para que un
modelo de perturbación lineal sea adecuado). Típicamente el eje rotará con una
excentricidad media que produce la aparición de fuerzas fluido dinámicas que contrarrestan
el peso del rotor y otras fuerzas radiales. Es importante reconocer que las fuerzas que el
fluido impone no solo resultan en fuerzas restitutivas cuando el flujo en el anulo está
dominado por efectos viscosos, que predominan a números de Reynolds . Dónde
, y . Lamentablemente no solo esto sucede, en caso
de que el signo de las fuerzas fluido dinámicas se revierte, tiende a magnificar la
excentricidad en lo que se suele llamar efecto Bernoulli o inercial. La presión en el fluido
decrece donde el huelgo entre rotor y estator o alojamiento es más pequeña ( ),
consecuentemente habrá una fuerza neta en el rotor en la dirección del desplazamiento.

En un sello corto o rodamiento el efecto predominante de la perturbación es causar

variaciones circunferenciales en la velocidad axial del fluido, que causa el efecto Lomakin,
donde es causada por y es la velocidad media promediada sobre el
huelgo. El efecto Lomakin es causado por variaciones circunferenciales en las perdidas de
entrada del flujo, en el lado con menor huelgo las perdidas serán menores debido a que
allí es menor. Consecuentemente la presión será mayor en el lado en que el huelgo es
menor resultando en una fuerza restitutiva debido a la distribución de presiones
circunferencial. El efecto Lomakin resulta en una fuerza fluidodinámica inducida que
conduce a una rigidez positiva .Los efectos presentes de Lomakin y Bernoulli
compiten entre si para definir la rigidez K del sistema mientras la geometría cambia. .

Mientras el rotor se aparta de su posición nominal, las fuerzas reactivas generadas por el
fluido pueden ser vinculadas al desplazamiento, la velocidad y la aceleración usando un
sistema de 2 grados de libertad, mediante un modelo masa resorte amortiguador.

Donde, Matriz masa; Matriz Amortiguamiento; Matriz rigidez; e movimiento

del rotor relativo al estator; y son componentes de la fuerza reactiva aplicada sobre el
rotor. Para cada condición de funcionamiento se requiere obtener la rigidez de los sellos
para la situación cuasi estática, caracterizado por los coeficientes y , y el comportamiento

Página 2 / 6 I 173
Departamento de Aeronáutica
 Facultad de Ingeniería - UNLP

dinámico caracterizado por los coeficientes de rigidez de amortiguamiento y de

inercia En casi todas las mediciones experimentales hechas en turbo máquinas las
matrices resultan ser antisimétricas, y por tanto invariantes ante rotaciones. Conocer Las
fuerzas del fluido variantes para las distintas condiciones de operación y de movimiento del
rotor nos permite describir las matrices , y a partir de la ecuación (1).

2 Parte Experimental

El enfoque convencional de este problema cuando es utilizar la ecuación

aproximada de Reynolds de fluido dinámica para una película delgada, que para el presente
contexto resulta en un fluido contenido entre dos cilindros circulares, de modo que dicha
ecuación resulta en

Donde son coordenadas radial, circunferencial y axial respectivamente, esta

ecuación una vez resuelta tiene por solución la distribución de presiones del fluido sobre el
rotor (promediada en dirección radial sobre el huelgo, para una separación entre
rotor y estator dada por , y una velocidad U del cilindro interior .
Una excentricidad, , a una frecuencia de whirling lleva a una separación dada por
de modo que substituyéndola en la anterior queda una expresión
para de forma tal que integrando la misma sobre el rotor se tienen las fuerzas roto
dinámicas que el mismo experimenta y a partir de estas los roto dinámicos vinculados a
estas por (1).

Este modelo simple que considera perturbaciones lineales en el cual se asume para el
primer término de la ecuación de Reynolds es predominante dado que se asume que las
perturbaciones mayormente ocurren en la dirección circunferencial más que en la axial que
puede despreciarse junto con los términos de orden superior en , mientras que si se asume
que las perturbaciones mayormente ocurren en la dirección axial más que en la
circunferencial dará como resultado aquellos otro set de valores para mientras
los valores .

Se ha encontrado que ocurren violentos movimientos del eje cuando la velocidad del eje del
rotor alcanza dos veces su valor critico del eje, esta es la respuesta dinámica del sistema
cuando se vuelve positiva, Condiciones similares ocurren a elevados . Es decir que si
, entonces será rotodinámicamente desestabilizante, mientras que si se
da el caso contrario el sistema tenderá a la estabilidad. Este primer modelo lineal prueba ser
típicamente estable pero resulta estar lejos del rango de operación de la Turbina Francis,
para lo cual se considero un otro modelo.

Modelos Semi analíticos: Sellos Anulares turbulentos.

En sellos anulares el flujo es usualmente turbulento debido al alto número de Reynolds

en el cual operan. Para identificar y modelar la rotodinámica de sellos turbulentos se usaron
modelos similares a las ecuaciones de lubricación de Reynolds. Estas emplean las
velocidades y que son promediadas sobre el huelgo. Fue Childs en 1983
quien publicó una versión revisada del modelo que surgía de estas consideraciones y utilizó
una generalización geométrica del mismo modelo para examinar la roto dinámica de las
corrientes de fuga de flujos descarga-succión alrededor de los impulsores de bombas
centrífugas y resultó más conveniente incluir una forma más general de análisis.

Página 3 / 6 174
Segundas Jornadas de Investigación y Transferencia - 2013

En la siguiente figura se representa un corte lateral del rotor centrífugo, de forma tal que se
describe le movimiento del fluido en coordenadas intrínsecas según el eje meridional del
huelgo dadas por y con donde la coordenada está medida a lo largo
del meridiano del huelgo.
La variación del huelgo se describe en términos de cuyo valor no perturbado es
. Las ecuaciones que gobiernan el flujo son promediadas sobre el huelgo, esto conduce
a las ecuaciones de continuidad de la forma y a las ecuaciones de cantidad de movimiento

El enfoque usado por Hirs (1973) es empleado para determinar las tensiones de corte
turbulentas y aplicadas al estator por el fluido en direcciones s y respectivamente y
las tensiones y aplicadas al rotor.

; 

Donde el valor del número de Reynolds en dirección meridional y las constantes

, , y son elegidos para ajustar la información y los datos disponibles sobre las
tensiones de corte turbulentas. Childs elije los siguientes típicos valores para dichas
constantes .Luego, el huelgo, la presión y las
velocidades son modeladas dividiéndolas en componentes medias ( de subíndice 0) que se
corresponderían con la ausencia de whirling, mas una pequeña componente de perturbación
lineal (de subíndice 1) debida a la excentricidad ε, y rotando a una frecuencia de whirling ω:

;

;
Estas ecuaciones sustituidas en las ecuaciones de continuidad y cantidad de movimiento
resultan en un set de ecuaciones para el flujo medio para las cantidades y un
segundo set de ecuaciones para las cantidades ; aquellos términos de orden
cuadrático u orden superior en ε son despreciados. Luego es necesario resolver ambos sets
de ecuaciones numéricamente para poder evaluar las presiones, y luego las fuerzas sobre
el rotor. Para una geometría simple de un sello anular plano, se tiene
y en el cual .

Donde es el caudal volumétrico; Childs (1983) fue capaz de obtener soluciones analíticas
para ambos sets de ecuaciones tanto las de valor medio como las de la perturbación. La
evaluación de las fuerzas roto dinámicas, conduce a los siguientes coeficientes roto
dinámicos:

; ;

Efecto Lomakin Efecto Bernoulli

Página 4 / 6 I 175
Departamento de Aeronáutica
 Facultad de Ingeniería - UNLP

Donde es el coeficiente de flujo ,y es la caída de presión total a

través del sello y dónde , están dadas por

; ; ; ;

; ; ;

Y Donde es el coeficiente de perdida dado por Yamada (1962).

3. Resultados y Discusión :

Considerando nuestros datos geométricos y de operación los siguientes

9,424777961 10 62,83185307 405300 70,3901 3 0,000875 0,0035

Caudal y
Cel 1+Cel + 2λ2
 ] [Kg/s] Área trasversal del sello
0,000001003 1013,85 0,065973446 0,1 181,3614728

Así Pues Definido , luego quedan definidos y luego , estando dado

por las condiciones de operación, lo mismo que y pues.

Parámetros independientes

Parámetros dependientes

En definitiva, para nuestras condiciones de operación

0,543515778 0,004482 90,13073638 1,406222 0,004302 0,481752 10,77318 0,084966

Y los coeficientes roto dinámicos serán

0 -0,038684342 1121,151017 2242,302033 763,4369971 763,4369971

Los mismos resultan ser desestabilizantes en dirección radial a la luz del valor , lo cual
indica un predominio del efecto de Bernoulli sobre el Efecto de Lomakin. Este es el resultado
de aplicar el modelo descripto semi-analítico para las condiciones de operación dadas, que
en rigor consideran no las fuerzas rotodinámicas sino las tensiones corte sobre el rotor y
estator así como la las presiones existentes entre entrada y salida.

Modelo numérico

Siguiendo el método de cálculo de coeficientes propuesto Por Athavale &Hendricks (1996) la

perturbación en la ubicación del rotor es usada para generar fuerzas reactivas en el fluido y
subsecuentemente coeficientes rotodinámicos. El modelo considera variaciones del flujo
tanto en direcciones axiales como circunferenciales pero los valores en dirección radial son
integrados para obtener los valores medios. Luego las ecuaciones de flujo son perturbadas
para resultar en un set de ecuaciones para el flujo estático y otro para las perturbaciones en
el flujo. Con un movimiento prescripto pequeño del rotor las variables de perturbación son
calculadas y las presiones integradas sobre la superficie del rotor para generar una fuerzas

Página 5 / 6 176
Segundas Jornadas de Investigación y Transferencia - 2013

reactivas del fluido mientras que curvas apropiadas se ajustan al grafico de Fuerzas vs
Frecuencia de whirling para obtener los coef rotodinámicos.
Se empezó por analizar el diversos casos cuasi-estáticos de desplazamiento del rotor en
dirección radial, pudiendo verificar algunos resultados obtenidos experimentalmente y
mediante simulación numéricas por otros investigadores citados validando las modelización
y consideraciones realizadas.

Conclusiones
Queda por modelizar el proceso completo de whirling a diversas relaciones entre y .
El proceso de análisis todavía está lejos de concluir dado que solo considera casos de matrices
antisimétricas exclusivamente lo cual si bien constituye el caso más general no es el único
posible.
Puede verse que también se obtienen valores principales de la componente Hidrostática de la
matriz rigidez que son desestabilizantes para el sistema redundando en la posterior
interferencia mecánica entre rotor y estator para algunas configuraciones operativas.

6. Referencias

1. M.M.ATHAVALE, R.C.HENDRICKS: “A Small Perturbation CFD Method for

Calculation of Seal Rotordynamic Coefficients,1996.
2. K. Graf & T. Staubli, “Comparison of two id models for rotordynamic coefficients of
labyrinth seals”. ETH Eidgenössische Technische Hochschule Zürich, 1990.
3. R. Gordon Kirk, “Fluid-Film Bearings, Seals and Rotor Stability” Mechanical
Engineering department, Virginia Politechnic Institute and State University; 2000?
4. C. Brennen, Radial and Rotordynamic Forces, Hydrodinamics of pumps, ISBN 0-19-
856442-2 (Oxford University Press)

<< Volver al Indice Página 6 / 6 I 177

Departamento de Aeronáutica