Formulario para Máquinas Axiales-TIM 42-Curso 2020

1. Ecuación del teorema de Euler:

𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 − ⃗⃗⃗⃗
𝑢1 . ⃗⃗⃗⃗
𝑣1
𝐻𝑡 =
𝑔
⃗⃗⃗⃗1 . ⃗⃗⃗⃗
a. Variante para el caso de ingreso de la velocidad sin prerrotación (𝑢 𝑣1 = 0) :
𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 𝑢2 . 𝑣2 . 𝑐𝑜𝑠(𝛼2 )
𝐻𝑡 = =
𝑔 𝑔

Notación: 𝑣2 . 𝑐𝑜𝑠(𝛼2 ) = 𝑣2𝑢 = 𝑣2𝜃

𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 𝑢2 . 𝑣2𝑢
𝐻𝑡 = =
𝑔 𝑔

Por otra parte, aún en el caso que no se cumpla ⃗⃗⃗⃗

𝑢1 . ⃗⃗⃗⃗
𝑣1 = 0, sigue valiendo:

𝑢2 = ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑢1 = 𝑢
⃗
Entonces:

𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 − ⃗⃗⃗⃗
𝑢1 . ⃗⃗⃗⃗
𝑣1 𝑢 ⃗⃗⃗⃗2 − ⃗⃗⃗⃗⃗⃗
⃗ . (𝑣 𝑣2 )
𝐻𝑡 = =
𝑔 𝑔

Además:

⃗⃗⃗𝑣𝑖 = ⃗⃗⃗⃗⃗
𝑣𝑖𝑟 + ⃗⃗⃗⃗⃗
𝑣𝑖𝑢

Siendo ⃗⃗⃗⃗
𝑣𝑖𝑗 la componente de la velocidad 𝑖-esima (𝑖 = 1,2 )
en la dirección 𝑗 (𝑗 = 𝑟, 𝜃, 𝑧 ).

Por las características de los rotores axiales la velocidad radial a la salida y a la entrada
son iguales, por lo tanto:

𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 − ⃗⃗⃗⃗
𝑢1 . ⃗⃗⃗⃗
𝑣1
𝐻𝑡 =
𝑔
𝑅⃗ = 𝐿 ⃗ +𝐷 ⃗
Operando:
⃗⃗⃗⃗
𝑅𝑢 = ⃗⃗⃗⃗
𝐿𝑢 + ⃗⃗⃗⃗
𝐷𝑢 = 𝐿. 𝑠𝑒𝑛(𝛽∞ ) + 𝐷. 𝑐𝑜𝑠(𝛽∞ )
1
𝐿 = . 𝐶𝐿 𝜌∞ 𝑣∞ 2 𝑙
2
 1
𝐷= . 𝐶 𝜌 𝑣 2𝑙
2 𝐷 ∞ ∞
1 1
⃗⃗⃗⃗
𝑅𝑢 = ⃗⃗⃗⃗
𝐿𝑢 + ⃗⃗⃗⃗
𝐷𝑢 = . 𝐶𝐿 𝜌∞ 𝑣∞ 2 𝑙. 𝑠𝑒𝑛(𝛽∞ ) + . 𝐶𝐷 𝜌∞ 𝑣∞ 2 𝑙. 𝑐𝑜𝑠(𝛽∞ )
2 2
1
⃗⃗⃗⃗
𝑅𝑢 = . 𝜌∞ 𝑣∞ 2 𝑙. (𝐶𝐿 𝑠𝑒𝑛(𝛽∞ ) + 𝐶𝐷 𝑐𝑜𝑠(𝛽∞ ))
2
1 𝐶𝐷 𝑐𝑜𝑠(𝛽∞ )
⃗⃗⃗⃗
𝑅𝑢 = . 𝜌∞ 𝑣∞ 2 𝑙. 𝐶𝐿 𝑠𝑒𝑛(𝛽∞ ) (1 + )
2 𝐶𝐿 𝑠𝑒𝑛(𝛽∞ )

Se tiene por definición:

𝐶𝐷
≝ tan 𝜆
𝐶𝐿

mientras que por trigonometría sabemos que:

𝑐𝑜𝑠(𝛽∞ )
= tan 𝛽∞
𝑠𝑒𝑛(𝛽∞ )
Queda entonces:

1 tan 𝜆
⃗⃗⃗⃗
𝑅𝑢 = . 𝜌∞ 𝑣∞ 2 𝑙. 𝐶𝐿 𝑠𝑒𝑛(𝛽∞ ) (1 + )
2 tan 𝛽∞

Por otra parte:

𝑃𝑡 = 𝛾𝑄𝑡 𝐻𝑡 = 𝜛𝑀𝑡

Pero para un diferencial de pala (𝑏 = Δ𝑟), donde estamos calculando la carga 𝐻𝑡 (𝑟)

La potencia disipada por unidad de ancho de pala (b o r), análoga al caso de fricción
entre sólidos se puede plantear como:

𝑃𝑑 𝐹𝐷
= .𝑣
𝑏 𝑏 ∞

Por otra parte, la potencia hidráulica disipada, por unidad de ancho de pala, es de la
forma usual (P=QH) que, para este caso particular, resulta:

𝑃𝑑
= 𝛾Δ𝑄𝑡 𝜉𝐻
𝑏
La potencia de la fuerza de sustentación, al ser perpendicular a la velocidad, es nula
⃗⃗⃗𝐿 . ⃗⃗⃗⃗⃗
(𝐹 𝑣∞ = 0).

Se destaca que se trabaja con una sección de la pala de perfil constante, de ancho r por
lo que se toma allí una parte del caudal Δ𝑄𝑡 .
Se puede aproximar:

Δ𝑄𝑡 = 𝑣𝑎 𝑡

Entonces:

𝑃𝑑
= 𝛾𝑣𝑎 𝑡𝜉𝐻
𝑏

Igualando:

𝑃𝑑 𝐹𝐷
= 𝛾𝑣𝑎 𝑡𝜉𝐻 = .𝑣
𝑏 𝑏 ∞

Notaremos las fuerzas por unidad de ancho de pala:

𝐹𝐷
𝐷≝
𝑏
Y

𝐹𝐿
𝐿≝
𝑏
Queda así:

𝑃𝑑 𝐹𝐷
= 𝛾𝑣𝑎 𝑡𝜉𝐻 = . 𝑣 = 𝐷. 𝑣∞
𝑏 𝑏 ∞

Entonces, despejando de la última ecuación obtenida:

𝛾𝑣𝑎 𝑡𝜉𝐻
𝐷=
𝑣∞
Por lo tanto se ve la relación entre el arrastre y las pérdidas.

Usando la definición de los coeficientes de fuerzas, en particular la de arrastre o “drag” se

puede plantear, para la grilla:

𝛾𝑣𝑎 𝑡𝜉𝐻 1
𝐷= = . 𝐾. 𝐶𝐷 𝜌∞ 𝑣∞ 2 𝑙
𝑣∞ 2
Para trabajar con el rendimiento hidráulico recordamos su definición:

𝐻 𝐻𝑡 − 𝜉𝐻 𝜉𝐻
𝜂𝐻 = = =1−
𝐻𝑡 𝐻𝑡 𝐻𝑡
 Usando las relaciones anteriores se puede obtener una expresión ( teniendo en cuenta
que 𝑣𝑎 = 𝑣∞ . 𝑠𝑒𝑛(𝛽∞ ) )

1 𝐾𝐶𝐷 𝜌∞ 𝑣∞ 3 𝑙 1 𝐾. 𝐶𝐷 𝑣∞ 2 𝑙
𝜉𝐻 = . = .
2𝑡 𝛾𝑣𝑎 2𝑔𝑡 𝑠𝑒𝑛(𝛽∞ )

Por otra parte, recordando la ecuación hallada para la carga:

𝑙 tan 𝜆
𝐻𝑡 = 𝑢. 𝑣∞ 𝐾. 𝐶𝐿 (1 + )
2𝑔𝑡 tan 𝛽∞

Sustituyendo en la expresión del rendimiento hidráulico:

𝑙 𝐾. 𝐶𝐷 𝑣∞ 2
𝐻 𝐻𝑡 − 𝜉𝐻 𝜉𝐻 2𝑔𝑡 . 𝑠𝑒𝑛(𝛽∞ )
𝜂𝐻 = = =1− =1−
𝐻𝑡 𝐻𝑡 𝐻𝑡 𝑙 tan 𝜆
2𝑔𝑡 𝑢. 𝑣∞ . 𝐾. 𝐶𝐿 (1 + tan 𝛽∞ )

Se tiene así:

𝐶𝐷 𝑣∞ 2
𝑠𝑒𝑛(𝛽∞ ) 𝐶𝐷 𝑣∞ 2
𝜂𝐻 = 1 − =1−
tan 𝜆 tan 𝜆
𝑢. 𝑣∞ . 𝐶𝐿 (1 + ) 𝑢. 𝑣∞ . 𝑠𝑒𝑛(𝛽∞ ). 𝐶𝐿 (1 + )
tan 𝛽∞ tan 𝛽∞

𝐶𝐷
Recordando que = tan 𝜆 y 𝑣𝑎 = 𝑣∞ . 𝑠𝑒𝑛(𝛽∞ ) se puede simplificar a las siguientes
𝐶𝐿
(indistintas) expresiones:

𝑣∞ tan 𝜆 𝑣∞ 2 tan 𝜆
𝜂𝐻 = 1 − =1−
𝑢. 𝑠𝑒𝑛(𝛽∞ ) (1 + tan 𝜆 ) 𝑢. 𝑣𝑎 (1 + tan 𝜆 )
tan 𝛽∞ tan 𝛽∞
𝑣∞ 2 tan 𝜆 𝑣𝑎 tan 𝜆
𝜂𝐻 = 1 − =1−
𝑢. 𝑣𝑎 (1 + tan 𝜆 ) 𝑢. [𝑠𝑒𝑛(𝛽∞ )]2 (1 + tan 𝜆 )
tan 𝛽∞ tan 𝛽∞
Procedimientos de cálculo para problemas de máquinas axiales: “diseño”

Se trata de determinar, en la sección significativa, la pareja (𝛽𝑃 , 𝑙𝑃 ) dados Q y H (o bien H

calculado a partir de Q, conocida la instalación, o sea dado Q, H=H(Q)).

Entonces:

𝐻 = 𝐻(𝑄) ≈ 𝐾. 𝑄 2 , por ejemplo

𝜂𝑉 ∶ 𝑑𝑎𝑡𝑜

𝜂𝐻 ∶ 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑟 𝑜 𝑝𝑢𝑒𝑑𝑒 𝑑𝑎𝑟𝑠𝑒 𝑐𝑜𝑚𝑜 𝑑𝑎𝑡𝑜 (𝑒𝑥𝑐𝑒𝑝𝑐𝑖ó𝑛)

Se suele utilizar el diagrama de Cordier para definir o verificar el diámetro a seleccionar

de una máquina que se desea diseñar.

∆= ∆(Ω)

Otros datos usuales:

Z,
perfil,
% de D en sección de diseño: 𝐷 = 𝐷𝑖 +0,6(𝐷𝑒 − 𝐷𝑖 )

Hipótesis comunes:
 trabajar a tan 𝜆𝑀𝐼𝑁 (optimiza el rendimiento): importante, establece un ángulo de
ataque, 𝛼 , fijo.
 𝐷𝑑 ≈ 𝐷𝑒
 Para iterar: suponer para arranque 𝑙𝑃 ≈ 0,1 m (no va a ser muy distinto de ese
orden)

Ejemplo de secuencia de iteración:

𝑒𝑀𝐴𝑋 𝑒𝑀𝐴𝑋 (%).𝑙𝑃 /100

Sup.: 𝑙𝑃 ⇒ 𝐶 = 1 − = 1−
𝑡 𝑡
Y también suponemos un 𝛽𝑃

Por el teorema de Euler:

𝑢2 . ⃗⃗⃗⃗
⃗⃗⃗⃗ 𝑣2 − ⃗⃗⃗⃗
𝑢1 . ⃗⃗⃗⃗
𝑣1 𝑢 ⃗⃗⃗⃗2 − ⃗⃗⃗⃗⃗⃗
⃗ . (𝑣 𝑣2 )
𝐻𝑡 = =
𝑔 𝑔

Usaremos la versión simplificada:

𝑢. ∆𝑣𝑢 𝑢. 𝑣2𝑢 𝑢. 𝑣2𝜃
𝐻𝑡 = ≈ ≝
𝑔 𝑔 𝑔
 Pero también podemos usar:

𝐻 𝐾. 𝑄 2
𝐻𝑡 = ≈
𝜂𝐻 𝜂𝐻

O bien:

𝐻 ≈ ∆𝑧 + 𝐾. 𝑄 2 = ∆𝑧 + ∆𝐻 = ∆𝑧 + 𝜉𝐻
Entonces:
𝐻𝑡 𝑔
∆𝑣𝑢 =
𝑢
Por lo visto en teórico:
Δvu

Δvu/2
v2 v∞
vR2 v1
vR1

u=u1 = u2 v∞, ABS

Además usaremos 𝑣𝑎 = 𝑣𝑎,𝐶

Velocidad axial corregida:

4𝑄
𝑣𝑎,𝐶 =
𝜋(𝐷𝑒 − 𝐷𝑖 2 ). 𝐶. 𝜂𝐻
2

Velocidad en ∞ (observar el triángulo dibujado más arriba y el teorema de Pitágoras):

∆𝑣𝑢 2
𝑣∞ = √𝑣𝑎,𝐶 2 + (𝑢 − )
2

De aquí podemos calcular:

𝑣𝑎
𝛽∞ ya que 𝑣∞ . sen 𝛽∞ = 𝑣𝑎, por lo tanto sen 𝛽∞ = 𝑣𝑎/𝑣∞ y 𝛽∞ = 𝑎𝑟𝑐𝑠𝑒𝑛𝑜(𝑣 )
∞

Impuesto 𝛼 , hallamos 𝛽𝑃 = 𝛼 + 𝛽∞ y también con ese valor de 𝛼 :

𝐶𝐷 = 𝐶𝐷 (𝛼) y
𝐶𝐿 = 𝐶𝐿 (𝛼)
O bien 𝐶𝐿 = 𝐶𝐿 (𝛼) y 𝐶𝐷 = 𝐶𝐿 . tan𝜆
En general aceptaremos, cuando estemos diseñando, utilizar: tan𝜆 = tan 𝜆𝑀𝐼𝑁
Entonces podemos calcular ya:
𝑙 tan 𝜆
𝐻𝑡 = 𝑢. 𝑣∞ 𝐾. 𝐶𝐿 (1 + )
2𝑔𝑡 tan 𝛽∞
Y
𝑣∞ tan 𝜆
𝜂𝐻 = 1 −
𝑢. 𝑠𝑒𝑛(𝛽∞ ) (1 + tan 𝜆 )
tan 𝛽∞

Métodos de “cierre de iteración”

Se puede cerrar de muchas formas siempre que todas las ecuaciones se interrelaciones y sean
compatibles. Una de las variables o varias van a converger en la iteración a un valor final que
aceptaremos como la solución de este sistema de ecuaciones no lineales.
Para el práctico aceptaremos un criterio de convergencia del 2%, por ejemplo:
𝑙∗ − 𝑙
< 0,02 (2%)
𝑙
Ejemplo de método de “cierre” de iteración:
Una vez hallado
𝑙 tan 𝜆
𝐻𝑡 = 𝑢. 𝑣∞ 𝐾. 𝐶𝐿 (1 + )
2𝑔𝑡 tan 𝛽∞
Podemos plantearlo como
tan 𝜆
𝑙 𝑢. 𝑣∞ 𝐾. 𝐶𝐿 (1 + tan 𝛽∞ )
𝐻𝑡 =
𝑡 2𝑔
Se puede reescribir como:

2𝑔𝐻𝑡 𝑡 𝑡
𝐾=( ) ≝𝑚
tan 𝜆
𝑢. 𝑣∞ . 𝐶𝐿 (1 + ) 𝑙 𝑙
tan 𝛽∞
Donde, evidentemente:

2𝑔𝐻𝑡
𝑚=( )
tan 𝜆
𝑢. 𝑣∞ . 𝐶𝐿 (1 + )
tan 𝛽∞
 Este planteo permite dibujar en el ábaco de Weinig: dos puntos y al unirlos una recta que al cortar con
la curva para el 𝛽𝑃 que corresponda nos permite obtener un nuevo valor del espaciamiento entre palas
𝑡 ∗
específico (llamémosle (𝑙 ) ).

2m

m m

𝑡 ∗
( )
𝑙

Finalmente (previo a volver a iterar) para este caso particular:

𝑡 ∗ 𝑡
( ) = ∗
𝑙 𝑙
debido a que 𝑡 está fijado por el número de palas y la sección de diseño por:
𝑡 = 𝜋𝐷/𝑍
Entonces:
𝑡
𝑙∗ =
𝑡 ∗
( )
𝑙
“Método del artillero” para mejorar convergencia de iteraciones (no deberíamos tener más de
3 vueltas completas):
-Empezamos a iterar con un valor supuesto: x => x*
-Volvemos a iterar con el valor con el que salimos e la iteración x= x* => (x*)* = x*
-Si no convergió la iteración planteamos un esquema con el del gráfico siguiente: si está bien
armado en el 99% de los casos converge a lo sumo (puede que antes) en este tercer paso.
Esquema gráfico para el “Método del artillero”

y=x

x*

x*

x x*

x (solución)
2. Procedimientos de cálculo para problemas de máquinas axiales: “operación”

Sería el camino opuesto al de diseño. En este caso se trata de determinar, en la sección

significativa, Q y H(r) dada ahora la pareja (𝛽𝑃 , 𝑙𝑃 ).
Sigue siendo válido obtener H calculado a partir de Q, conocida la instalación, o sea dado
Q, H=H(Q)).

Entonces, ahora el caudal es desconocido, pero podemos suponer uno y empezar a iterar
hasta que converja:

Supuesto 𝑄, 𝐻 = 𝐻(𝑄) ≈ 𝐾. 𝑄 2 , por ejemplo

𝜂𝑉 ∶ 𝑑𝑎𝑡𝑜 (siempre)

𝜂𝐻 ∶ 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑟 𝑜 𝑝𝑢𝑒𝑑𝑒 𝑑𝑎𝑟𝑠𝑒 𝑐𝑜𝑚𝑜 𝑑𝑎𝑡𝑜 (𝑒𝑥𝑐𝑒𝑝𝑐𝑖ó𝑛)

Podemos suponer también un valor de 𝜂𝐻 para calcular:

𝐻
𝐻𝑡 =
𝜂𝐻
De allí podríamos obtener:

𝐻𝑡 𝑔
∆𝑣𝑢 =
𝑢

Otros datos usuales:

Z,
perfil,
% de D en sección de diseño: 𝐷 = 𝐷𝑖 +0,6(𝐷𝑒 − 𝐷𝑖 )

Hipótesis comunes:
 𝑄 tal que velocidad en ducto sea por ejemplo 10 m/s o 𝑄 ≈ 1 𝑚3 /𝑠 si estamos
cerca de ese valor para acelerar las cuentas.
 Trabajar cerca de tan 𝜆𝑀𝐼𝑁 , pero teniendo en cuenta que eso NO tiene por que
cumplirse al cerrar la iteración. Con ello también se pueden acortar cuentas ya que
establece un ángulo de ataque, 𝛼 , fijo al menos para la primera vuelta de cálculos.

Ejemplo de secuencia de iteración:

Sup.: 𝑄 ⇒ 𝐻 ⇒ 𝐻𝑡 ⇒ ∆𝑣𝑢 ⇒ 𝑣𝑎,𝐶 ⇒ 𝑣∞ (etc)

Ahora conocemos 𝛽𝑃
Igual que antes la velocidad axial corregida es calculable por:
4𝑄
𝑣𝑎,𝐶 =
𝜋(𝐷𝑒 − 𝐷𝑖 2 ). 𝐶. 𝜂𝐻
2

𝐶 ahora es fijo porque el largo de pala y sección significativa son datos.

Velocidad en ∞:

∆𝑣𝑢 2
𝑣∞ = √𝑣𝑎,𝐶 2 + (𝑢 − )
2

De aquí podemos calcular:

𝑣𝑎
𝛽∞ ya que 𝑣∞ . sen 𝛽∞ = 𝑣𝑎, por lo tanto sen 𝛽∞ = 𝑣𝑎/𝑣∞ y 𝛽∞ = 𝑎𝑟𝑐𝑠𝑒𝑛𝑜(𝑣 )
∞

Se puede calcular 𝛼 , prefijado 𝛽𝑃 , como: 𝛼 = 𝛽𝑃 − 𝛽∞ y tambien con ese valor de 𝛼 :

𝐶𝐷 = 𝐶𝐷 (𝛼) y 𝐶𝐿 = 𝐶𝐿 (𝛼)
O bien 𝐶𝐿 = 𝐶𝐿 (𝛼) y 𝐶𝐷 = 𝐶𝐿 . tan𝜆
Entonces podemos calcular ya:
𝑙 tan 𝜆
(𝐻𝑡 )∗ = 𝑢. 𝑣∞ 𝐾. 𝐶𝐿 (1 + )
2𝑔𝑡 tan 𝛽∞
Y
𝑣∞ tan 𝜆
(𝜂𝐻 )∗ = 1 −
𝑢. 𝑠𝑒𝑛(𝛽∞ ) (1 + tan 𝜆 )
tan 𝛽∞
𝐻
Por lo tanto, se puede hallar H* (a partir de 𝐻𝑡 = 𝜂 ) y entonces obtener un nuevo Q* para volver a iterar
𝐻
ya que si fuera:

𝐻 ≈ 𝐾. 𝑄 2
Entonces se puede despejar directamente:

𝑄 ∗ = √𝐻 ∗ /𝐾

Y se puede volver a iterar entrando con ese valor y así hasta que converja.
Métodos de “cierre de iteración”
Al igual que antes, pero por ejemplo:
𝑄∗ − 𝑄
< 0,02 (2%)
𝑄