I.

RESUMEN:
En el presente proyecto, se realizó un estudio de una turbina hidráulica Michell Banki, que
según la teoría de turbomáquinas hidráulicas es una turbina de acción de flujo transversal
y admisión parcial. Al pasar dos veces se obtiene una eficiencia elevada para flujos
variables, además de limpiar el rotor de residuos. La máquina es de velocidad baja, ideal
para bajas alturas pero elevados caudales, gracias a su simplicidad constructiva suelen ser
máquinas de coste reducido. Para realizar el análisis de la turbina se la modeló usando el
programa Solidworks con sus dimensiones reales para posteriormente mediante el
complemento Flow Simulation comenzar el estudio fluidodinámico. Se determinó las fuerzas
que intervienen en la turbina, además del torque y las velocidades que se producen en la
salida.

La modelación computacional del flujo en el interior de turbinas Michell-Banki permite

analizar y evaluar la incidencia del flujo en álabes para diferentes situaciones. A través de
estos resultados se pueden plantear modificaciones en el diseño analizando la energía
transferida por parte del fluido a la turbomáquina y su eficiencia hidráulica. Los resultados
muestran que se puede obtener un incremento en la eficiencia, modificando el tamaño de
los álabes para obtener la mayor eficiencia de la máquina.

II. OBJETIVOS
 Construir y ensamblar en SOLIDWORKS una turbinada MICHELL BANKI.

 Simular la turbina MICHELL BANKI, usando el programa flow simulation, para

poder conseguir valores de diferentes variables físicas que su conocimiento es de
importancia en la turbina.

III. IMPORTANCIA Y/O JUSTIFICACIONES:

La turbina de flujo transversal o turbinas Michell Banki es una maquina utilizada
principalmente para pequeños aprovechamientos hidroeléctricos. No obstante esto no
impide que la turbina se utilice en grandes instalaciones. Aunque la turbina de flujo
transversal se conoce como una turbina de pequeña escala, existen actualmente máquinas
de este tipo de hasta 6 MW.
Las principales características de esta máquina son las siguientes:
 La velocidad de giro puede ser seleccionada en un amplio rango.
 El diámetro de la turbina no depende del caudal.
 Se alcanza un aceptable nivel de rendimiento con pequeñas turbinas.
 Se puede regular el caudal y la potencia por medio de un alabe ajustable.
El desarrollo de este trabajo se justifica por adquirir conocimiento sobre la Turbina Michell
Banki, realizar su simulación y ver cómo es que se comporta un flujo de agua, analizando
velocidades, fuerzas, y caudal.
IV. FUNDAMENTO TEÓRICO:

Turbina Hidráulica

Las turbinas hidráulicas son la evolución natural de la sencilla rueda hidráulica de Poncelet;
sin embargo se diferencian en su construcción y en las mayores velocidades obtenidas en
las turbinas de forma que su eficiencia es mucho más elevada que el de la rueda hidráulica
de Poncelet. Así como una bomba absorbe energía mecánica y transfiere dicha energía al
fluido: una turbina absorbe energía del fluido y transfiere dicha energía en algún movimiento
mecánico.
La función de una turbina y de toda máquina hidráulica es efectuar un cambio de energía
entre un sistema mecánico y un sistema fluido.

Según la variación de la presión estática a través del rodete:

Turbinas de impulsión o de acción.

Actualmente se le llama Turbina de impulsión a aquella en la que la energía de presión o
potencial del agua se convierte en energía cinética antes de que esta agua incida sobre
una limitada porción periférica de un elemento rotativo, sin que haya un cambio posterior
de presión
Entre las Turbinas de Impulsión tenemos:

 Turbina Michell Banki (Cross - flow)

 Turbina Pelton de uno o más inyectores
 Turbina Turgo

Turbinas de reacción.
La conversión inicial presión - velocidad se realiza sólo parcialmente, de forma que el agua
entra en el elemento rotativo por toda la periferia y el flujo pasa por todos los espacios libres.
A diferencia de las bombas, la mayoría de las turbinas han de funcionar a potencia distinta
de la normal durante considerables periodos de tiempo, haciéndose frente a las
variaciones de carga, mediante la regulación de la cantidad de agua, pero manteniendo
constante la velocidad de sincronismo.
Dentro de las Turbinas de Reacción tenemos:
  Turbinas Francis en sus variantes: lenta, normal y rápida.
 Turbina Deriaz.
 Turbina Kaplan y de hélice.
 Turbinas axiales, en sus variantes: tubular, bulbo y de generador
periférico.

Campos de acción en términos de caudal y altura de turbinas hidráulicas

La turbina de flujo transversal o Michell-Banki

Es una máquina utilizada principalmente para pequeños aprovechamientos hidroeléctricos.

Sus ventajas principales están en su sencillo diseño y su fácil construcción lo que la hace
atractiva en el balance económico de un aprovechamiento a pequeña escala. No obstante
esto no impide que la turbina se utilice en grandes instalaciones. Aunque la turbina de flujo
transversal se conoce como una máquina de pequeña escala, existen actualmente
máquinas de este tipo de hasta 6 MW.

Las principales características de esta máquina son las siguientes:

 La velocidad de giro puede ser seleccionada en un amplio rango.

 El diámetro de la turbina no depende necesariamente del caudal.
 Se alcanza un aceptable nivel de rendimiento con pequeñas turbinas.
 Se puede regular el caudal y la potencia por medio de un álabe ajustable.

Principio de funcionamiento
La turbina consta de dos elementos principales: un inyector y un rotor. El agua es restituida
mediante una descarga a presión atmosférica. El rotor está compuesto por dos discos
paralelos a los cuales van unidos los álabes curvados en forma de sector circular.

El inyector posee una sección transversal rectangular que va unida a la tubería por una
transición rectangular - circular. Este inyector es el que dirige el agua hacia el rotor a través
de una sección que toma una determinada cantidad de álabes del mismo, y que guía el
agua para que entre al rotor con un ángulo determinado obteniendo el mayor
aprovechamiento de la energía. La energía del agua es transferida al rotor en dos etapas,
lo que también da a esta máquina el nombre de turbina de doble efecto, y de las cuales la
primera etapa entrega un promedio del 70% de la energía total transferida al rotor y la
segunda alrededor del 30% restante. (Los ensayos realizados por distintos investigadores
sitúan el rendimiento hidráulico de esta máquina entre un 65-70%, otros autores mencionan
un 61% aclarando que la segunda etapa entrega un 17%, y en general muchos autores
indican un 70% hasta un 84%.) Una característica atractiva de esta máquina es la forma
aplanada de su curva de rendimiento. Esto se logra con un diseño de la turbina con
admisión parcial. Por ejemplo: si se divide el rotor en 3 partes iguales y la admisión del agua
se puede realizar por 1/3, 2/3 o la totalidad del rodete.
 Turbina De Flujo Transversal De Admisión Parcial

Como se mencionó anteriormente, la turbina de flujo transversal es especialmente

apropiada para ríos con pequeños caudales. Estos generalmente llevan durante varios
meses muy poca agua, por lo que en su diseño debe considerarse para el mínimo caudal
que será el parcial y para épocas de abundancia de agua, se considerará el caudal total
que será utilizado para usos productivos. De la curva característica de rendimiento de cada
turbina depende si durante este tiempo se sigue produciendo energía eléctrica.

Descripción de los componentes

A continuación de describen los componentes de la turbina Michell-Banki como se

puede observar en la siguiente figura:

1. Carcasa

La carcasa de la turbinas Michell-Banki es totalmente construida en acero

extremadamente robusto y resistente a los impactos, para evitar posibles
deformaciones en todo el conjunto. Para obtener una estructura rígida las chapas
son ajustadas y abulonadas antes de la soldadura.
2. Inyector

El inyector posee una sección transversal rectangular que va unida a la tubería por
una transición rectangular - circular. El inyector es el que dirige el agua hacia el rotor
a través de una sección que toma a una determinada cantidad de álabes del mismo,
y que guía el agua para que entre al rotor con un ángulo determinado obteniendo el
mayor aprovechamiento de la energía.
Esta conducción deberá poseer una buena aceleración y una distribución de
velocidades uniforme en la sección de salida, así como un bajo nivel de pérdidas de
carga, para lograr la mayor transformación posible de energía potencial en energía
cinética.
El inyector puede tener distintas geometrías diferenciadas fundamentalmente por el
ángulo de admisión y el órgano de regulación que posea, si es que éste existe.
Rotor
Se trata de la pieza más significativa de esta turbina. Está equipado con una serie
de álabes de acero lijado y perfectamente pulido. Los álabes están soldados a unos
discos laterales mediante un proceso que precisa un especial cuidado. Existen
estudios realizados que presentan un número óptimo de álabes que están en
función del diámetro del rotor, éstos varían ente 22-28(Instituto Nacional de Energía,
1986). El perfil curvado de estos álabes, hace que se produzca una fuerza axial
limitada, de forma que se puede prescindir de un eje axial reforzado, que podría
acarrear una complicada fijación.
En rotores anchos los álabes son reforzados con varios discos. Antes del montaje
final se equilibran los rodetes cuidadosamente y se realiza un control
defectoscópico.

3. Eje

Las turbinas Michell-Banki están equipadas con un eje horizontal pasante, anclado
bilateralmente a dos cojinetes de empuje. Es decir que el flujo pasa a través de la
parte interna del rotor. Esta construcción, por más que se tomen cuidados, no
elimina la interferencia que causa el eje al flujo del agua, lo que disminuye el
rendimiento.
Los cojinetes absorben todos los esfuerzos axiales derivados de la turbina, para
posteriormente transmitirlos a la cimentación. Al mismo tiempo el rodete está
centrado en el eje, y por tanto, entre ambos cojinetes.
Otro tipo de eje es el que se llama, eje bridado, que exige un cuidado mayor durante
la fabricación, tanto en el posicionamiento y como en la fijación de las bridas a los
discos laterales del rotor. Estas deben estar alineadas para el montaje de los
rodamientos, con el objeto de evitar posibles vibraciones indeseables para la
estabilidad del conjunto. En este tipo de eje es importante hacer una comprobación
de la rigidez de la pieza, ya que ésta se ve perjudicada por la falta del eje pasante
(Gonçalves Chávez, y otros, 2007).
Campo de aplicación

La turbina Michell-Banki permite aprovechar económicamente saltos de baja

potencia siendo por esto recomendable para su uso en pequeñas centrales.

Fundamentalmente su aplicación se destina a la producción de energía eléctrica en

pequeña escala, o en otros casos, su eje se acopla por correa a otros dispositivos
mecánicos, y la energía mecánica obtenida se utiliza directamente en trabajos de
taller.

Análisis del diagrama de velocidades de la turbina Michell Banki.

Los perfiles de los alabes del rodete de la turbina Michell Banki, se determinan en base a
los diagramas de velocidades en cada punto del rodete. Para determinar estos diagramas,
es necesario definir la velocidad de salida del agua del inyector, la que se determinan en
base a la ecuación de Bernoulli aplicada entre la superficie del reservorio, donde la
velocidad del agua es aproximadamente cero y la salida del inyector.
P0 C02 P1 C12
+ + Z0 = + + Z1 + ∆Ht + ∆H1
γ 2g γ 2g
En donde:
 Co y C1: representa la velocidad de una partícula de agua en la superficie del
reservorio y a la salida del inyector, respectivamente.
 Po y P1: representa las presiones en la superficie del reservorio y en la salida
del inyector, respectivamente. En este caso ambas presiones se pueden
 considerar iguales a la atmosférica.
 Zo y Z1: representa los niveles topográficos, en la superficie del reservorio y la
posición del inyector, respectivamente, y su diferencia es igual al salto bruto.
 ɣ y g: representan el peso específico del agua y la aceleración de la gravedad,
respectivamente.
 ΔHt: es la pérdida de presión por efecto de la fricción del agua con las paredes
de la tubería de presión.
 ΔH1: es la pérdida de presión por efecto de la fricción del agua con las paredes
del inyector.
Con todas estas consideraciones se determina que la velocidad del agua a la salida del
inyector es:

∆H1
C1 = √1 − . √2gH
H

Dónde:
H: es el salto efectivo o neto, obtenido de la diferencia entre el salto bruto y las pérdidas de
presión en la tubería.

 Se define Kc, como coeficiente de velocidad del inyector. Este coeficiente de velocidad
Kc afecta a la velocidad absoluta de entrada y tiene en cuenta las pérdidas que se generan
en el escurrimiento dentro del inyector. De acuerdo a esto se puede decir que el coeficiente
Kc afectará de manera directa al rendimiento hidráulico de la turbina.
Cuando Kc se aleja de la unidad decreciendo su valor (lo que equivale a un inyector
ineficiente) se hace necesario un incremento en el ángulo de admisión de la turbina.
El coeficiente Kc es un valor determinado generalmente en forma experimental. Al igual que
en la turbina Pelton donde Kc toma valores entre 0,97 y 0,98, en una turbina de doble efecto
toma valores menores a la unidad.
Los distintos investigadores han obtenido valores de kc que van desde 0,95 a 0,98.

Por lo tanto el Kc, coeficiente de velocidad del inyector queda expresado por:

∆H1
K c = √1 −
H

Con lo cual la velocidad de salida del agua del inyector queda expresada por:
C1 = K c . √2gH
La velocidad del agua a la salida del inyector es igual a la velocidad de ingreso del agua al
rodete. Este chorro de agua a su vez se orienta hacia el rodete con un ángulo promedio
denominado α2, el cual posee valores prácticos de que se encuentran alrededor de los 16
grados.

Fig. 3. Triangulo de velocidades del rodete de una Turbina Michell Banki

De acuerdo a la Ecu. 𝐂𝟏 = 𝐊 𝐜 . √𝟐𝐠𝐇 , procederemos analizar el triángulo de velocidades

que se muestran en la Fig. 4 y se deducen las siguientes fórmulas.

Fig. 4. Triángulo de velocidades unificados de 1° y 2 ° etapa.

𝐖𝟏𝟐 = 𝐂𝟏𝟐 + 𝐔𝟏𝟐 − 𝟐𝐂𝟏 . 𝐔𝟏 . 𝐜𝐨𝐬 𝛂𝟏

 𝑚
W = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 ( )
𝑠
𝑚
C = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 ( )
𝑠
𝑚
u = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑡𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙 𝑑𝑒 𝑙𝑎 𝑡𝑢𝑟𝑏𝑖𝑛𝑎 ( )
𝑠
α = 𝑎𝑛𝑔𝑢𝑙𝑜 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 𝑦 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑡𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙 (𝑔𝑟𝑎𝑑𝑜𝑠°)

Por lo tanto la máxima utilización de la energía:

Cul 2 Cul
W12 = C12 + ( ) − 2C1 . ( ) . cos α1
2 2
𝑚
Cu = 𝑐𝑜𝑚𝑝𝑜𝑛𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 𝑒𝑛 𝑙𝑎 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛 𝑡𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙 ( )
𝑠

Y observando el triángulo de velocidades tenemos

C1 . cos α1 2 Cul
W12 = C12 + ( ) − 2C1 . ( ) . cos α1
2 2

Por lo tanto obtenemos:

3 cos 2 α1
W12 = C12 (1 − ) … (a)
4
Además se cumple:
Cm1 = C1 . sin α1 … (b)
Cm1 = W1 . sin β1 … (c)
Cm = 𝑐𝑜𝑚𝑝𝑜𝑛𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 𝑒𝑛 𝑙𝑎 𝑑𝑖𝑟𝑒𝑐𝑐𝑖ó𝑛 𝑚𝑒𝑟𝑖𝑑𝑖𝑎𝑛𝑎.

Relacionando las ecuaciones (a) y (b) con (c) se obtiene:

sin α1
β1 = arcsin
√1 − 3 . cos2 α1
4

β = á𝑛𝑔𝑢𝑙𝑜 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑎 𝑦 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑡𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙 (𝑔𝑟𝑎𝑑𝑜𝑠 °).

α = á𝑛𝑔𝑢𝑙𝑜 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 𝑦 𝑙𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑡𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙 (𝑔𝑟𝑎𝑑𝑜𝑠 °).

El ángulo α1 varía generalmente entre los 14º y los 17º. Muchos autores utilizan este ángulo
con valores que no sobrepasan los 16º. De acuerdo a esto el ángulo β1 variará de acuerdo
a los valores mostrados en la Tabla 2.
El ángulo influye β1 en la construcción del rotor lo que lleva a tratar de utilizar valores que
faciliten dicha construcción. De acuerdo a esto será conveniente adoptar un ángulo α1 de
16,102º con lo se obtendrá de la expresión (15) un ángulo β1 de 30º.

Tabla 2

Como se mencionó anteriormente la velocidad en el inyector está dada por la ecuación:

𝐂𝟏 = 𝐊 𝐜 . √𝟐𝐠𝐇.

Reduciendo esta expresión se obtiene:

C1 = 4.429 K c . √H

Considerando la hipótesis de impulso y de acuerdo al triángulo de velocidades a la entrada

se tiene:
Cul C1 4.429 K c . √H. cos α1
U1 = = . cos α1 =
2 2 2

Luego:
U1 = 2.214. K c . √H. cos α1

Y si α1 =16°, la velocidad tangencial será:

𝑈1 = 2.127. 𝐾𝑐 . √𝐻
Por lo tanto la velocidad relativa a la entrada del rotor será:

W1 = 2.458. K c . √H
Después de realizar el análisis de los triángulos de velocidades el cual es un punto
importante para describir el perfil de los alabes del rotor, pasaremos a mencionar las
formulas y consideraciones que se tendrá para el diseño de la turbina Michell Banki.

Diseño del Rodete.

A. Rendimiento y Potencia de la Turbina.

1. Rendimiento interno teórico: Asumiendo ángulos de ataque del inyector (1) y su
coeficiente de pérdida (Kc).
0.95 ≤ 𝐾𝑐 ≤ 0.98
14° ≤ 𝛼1 ≤ 18°
𝛈𝐢𝐭 = 𝐊 𝐜 𝟐 ( 𝟏 − 𝐓𝐚𝐧𝟐 𝛂𝟏 )

2. Rendimiento interno: Considerando el coeficiente de fuga de caudal: 𝑎 = 0.95

𝛈𝐢 = 𝐚 ∗ 𝛈𝐢𝐭
3. Rendimiento de la turbina: Asumiendo que el rendimiento mecánico 𝜂𝑚 es entre (0.85
– 0.92) según la calidad constructiva:
𝛈𝐭 = 𝜼𝒊 ∗ 𝛈𝐦

4. Potencia nominal al eje: Se considera caudal (Q) y altura (H) de diseño.

𝐏 = 𝐠. 𝐐. 𝐇. 𝛈𝐭

B. Determinación de la Rotación.

1. Criterio de la Velocidad de Rotación.

Matemáticamente la velocidad de rotación óptima de una Turbina Michell Banki esta

expresada por la siguiente fórmula:

√H
N = 40 (rpm)
D
Dónde:
H: se refiere a la altura neta.
 D: Se refiere al diámetro del rodete.

2. Cálculo de la Velocidad especifica del caudal.

3 N Q1/2
nq =
H 3/4

N: rpm deseada de la turbina:

Donde  q debe salir dentro del siguiente rango: 40 ≤ 𝑛𝑞 ≤ 180

3. Cálculo de ángulos de ataque entre el alabe y el rotor.

Este ángulo se entre la velocidad tangencial y la relativa a la entrada del rotor.

β1 = 2 ∝1

C. Cálculo de velocidades a la entrada del rodete, diámetro del rodete

1. Velocidad absoluta a la entrada del rotor.

𝑚
𝐶1 = 4.429 𝐾𝑐 √𝐻 ( )
𝑠

2. Velocidad tangencial a la entrada del rotor.

𝑈1 = 2.127 𝐾𝑐 √𝐻

3. Diámetro externo del rotor.

𝑈1
𝐷 = 19.08 (𝑚)
𝑁
4. Diámetro interno del rotor.
𝐷𝑖 = 0.58 𝐷 (𝑚)

DISEÑO DE LOS ÁLABES

Números de álabes del rotor.

La selección del número de álabes se realizará en base al diámetro, espesor del material y
las condiciones de funcionamiento de la turbina, es decir, altura y caudal.
Se deberá tener en cuenta que un reducido número de álabes provocará pulsaciones en la
generación de la potencia, y un número elevado producirá una aceleración de la vena fluida
con el consiguiente aumento de las pérdidas y el efecto de reja.
 𝐷
𝑍 = 0.10 𝜋
𝑒𝑝

Dónde: ep = Equivale al espesor del material a utilizar.

En la tabla 3 se transcriben los resultados de diferentes investigaciones reflejadas en la

literatura. De esta tabla se desprende que el número óptimo de álabes está entre 24 y 30.

Tabla 3: Parámetros característicos en turbinas ensayadas en diferentes investigaciones.

Paso externo de los álabes.

D
tt = π (m)
Z

Radio medio de los álabes.

Tomando X4 como r ó 1/3 D

(𝐷2 + 𝐷𝑖2 − 2𝐷𝑖 ∗ √𝐷2 − 4 𝑋42 )

𝑟= (𝑚)
8𝑋4
DETERMINACIÓN DE LAS DIMENSIONES DEL INYECTOR Y ALABE DIRECTRIZ.

Para definir la geometría del inyector es necesario considerar en el diseño una buena
conducción y aceleración del flujo de agua, así como una adecuada orientación y regulación
de éste flujo hacia los álabes del rodete.
Para las dimensiones del inyector y del alabe directriz se tendrá en cuenta la tabla que a
continuación se adjunta.
COTA DIAMETRO DEL RODETE EN mm
(mm) 200 300 400 500
a 261 348 435 522
b 195 260 325 390
c 31 41 52 62
d 102 136 170 204
e 85 113 142 170
f 55 73 92 110
RI 168 224 280 336
RII 151 201 252 302
RIII 28 37 47 56
R1 133 177 138 161
R2 75 100 96 111
R3 110 110 138 160
R4 70 80 83 126
 R5 35 50 52 66
R6 45 77 52 29
R7 100 90 101 126
x 40 47 48 63
y 120 136 147 211
L1 98 111 122 179
L2 60 70 90 136
L3 24 27 25 64

El chorro entra al rotor con un ángulo α2 que es constante en toda la admisión y tangente a
la periferia de la rueda. El flujo que abandona las paredes sólidas del inyector es definido
como un chorro libre. La velocidad a la salida del inyector tiene un valor un poco más
pequeño que el valor de diseño lo que provoca un incremento en el arco de entrada.
Como ya se mencionó, la diversidad de diseño en la geometría del inyector hace que se
adopten distintos ángulos de admisión. A través de las diversas investigaciones que se han
realizado sobre esta máquina los ángulos de admisión del inyector van desde los 30º hasta
los 120º.
Gran parte de la bibliografía existente parece coincidir en que el ángulo de admisión θ ad
óptimo para este tipo de turbina es de alrededor de los 90º.
Además se deberá tener en cuenta el "efecto de reja" que causa el espesor de los álabes
en la entrada.

Arco de admisión:
πD
Lad = . Zad + Zad . e
Z
Luego el ángulo de admisión θad será:
 360°
θad = . Lad
πD
Por último nos queda hallar la función que representará la curva envolvente del inyector
(cara superior). El modelo matemático de la entrada y salida del caudal en el inyector puede
definirse como un flujo potencial.
En general para cualquier ángulo entre 0º y θad, la curva envolvente del inyector está dada
por:
θ
1−
θad
[ ]Q
B∗C
τθ = R ∗ e
En el que:
𝐷𝑒 . √𝐻
𝐶 = 2.3 𝜂ℎ
𝐾𝑐

Luego la altura del inyector en cada punto de la envolvente será:

h t = τθ − R
Dónde:
De
R=
2
En el cual:
τθ = Radio de la curva envolvente del inyector (m)
𝜃 = Ángulo de la envolvente en un punto cualquiera
𝜃𝑎𝑑 = Ángulo de admisión.
Q = Caudal de diseño (m3/s)
B = Ancho del rotor, (m)
𝜂𝑡 = Eficiencia hidráulica.
De = Diámetro del rodete o rotor (m)
H = Altura neta, (m)
K c = Coeficiente de velocidad del inyector.
e = espesor de los álabes

El área de la admisión entonces como se sabe será, m2:

Q
𝐴𝑎𝑑 =
C2
En donde:
C2 = Velocidad del agua a la salida del inyector, (m/s)
Ancho del inyector

𝑄
𝐵𝑖 = 0.96 (𝑚)
(√𝐻 ∗ 𝐷)
Ancho del rotor.

𝐵 = (1.2 − 1.5) 𝐵𝑖 (𝑚)

Diámetro primitivo.

Dp = 2 √r 2 + (0.5 Di )2 (𝑚)

INSTALACIONES COMPLEMENTARIAS AL DISEÑO (EJE)

Diseño del eje del rotor:

 Distancia del filete interno a la línea del centro del rotor.

−1
D
L = 0.5Di [√1 + sen2 β1 ( )] (m)
Di

 Diámetro mínimo del eje del rotor.

dmin = 0.5 L (m)

V. MODELO COMPUTACIONAL DE LA TURBINA
MICHELL-BANKI
Existen dos métodos fundamentales para diseñar y analizar sistemas de ingeniería que son:
la experimentación con la construcción de modelos y el cálculo, que implica la resolución
de ecuaciones. Ambos métodos se complementan entre sí, obteniendo propiedades
globales (potencia, presión, fuerza de arrastre, etc.) experimentalmente y mediante el
cálculo, teniendo detalles acerca del campo del flujo como perfiles de presión, velocidad,
líneas de corrientes y otros. Las cantidades globales obtenidas mediante el cálculo son
posteriormente verificadas con los datos experimentales.

La Dinámica Computacional de Fluidos, (CFD por sus siglas en inglés) es el cálculo de

ecuaciones diferenciales del flujo de fluidos mediante el uso de computadoras que a la vez
permiten hacer la representación del flujo y predecir su comportamiento. Usando CFD es
posible construir un modelo computacional que represente el desarrollo del flujo en un
sistema o en un equipo que se quiera estudiar. Después se especifican las condiciones
físicas y químicas del fluido al prototipo virtual y el software entregará la predicción de la
dinámica del fluido. La herramienta CFD ofrece la capacidad de simular flujos de gases,
líquidos, transferencia de masa y calor, cuerpos en movimiento, física multifases,
reacciones químicas, interacción fluido-estructura y acústica a través de la modelación en
el computador.

Procedimiento para el análisis

1º. Se define un dominio computacional y se genera una malla o red de nodos. Éstas definen
las celdas, para el caso de un domino bidimensional son áreas y para un dominio
tridimensional son volúmenes.

2º. Se especifican Condiciones de frontera en cada lado del dominio (flujos 2D) o en cada
cara del domino (flujos 3D).

3º. Se especifica el tipo de fluidos junto con las propiedades del mismo.

4º. Se especifican los valores para una primera aproximación de las variables de entrada.

5º. Establecer un criterio de convergencia de la solución, especificando un número máximo

de iteraciones y un valor mínimo para el residuo que se genera en la solución de las
ecuaciones discretizadas.

6º. Cuando converge la solución se expresan las variables de campo mediante gráficas,
archivos Excel, que mejoran la visualización de resultados.
VI. SIMULACIÓN
a. PARA EL ALABE DIRECTRIZ UBICADO A 0°

Presión

Temperatura de fluido
Velocidad en x

Velocidad en y
 b. PARA EL ALABE DIRECTRIZ UBICADO A 15°

Presión

Temperatura de fluido
Velocidad en x

Velocidad en y
 c. PARA EL ALABE DIRECTRIZ CERRADO

Presión

Temperatura de fluido
Velocidad en x

Velocidad en y
 VII. RESULTADOS
a) PARA EL ALABE DIRECTRIZ UBICADO A 0°

MICHELLE BAN.SLDPRT.SLDASM [ANALISIS DE TURBINA M.BANKI [0]]

Goal Name Unit Value Averaged Value Minimum Value Maximum Value Delta Criteria
GG Av Temperature (Fluid) 1 [K] 293.216809 293.2168007 293.2167886 293.2168098 2.1171E-05 0.002588134
GG Mass Flow Rate 1 [kg/s] -0.000833472 -0.000159962 -0.006849341 0.00795454 0.002051966 5.619611701
GG Av Velocity 1 [m/s] 2.982730081 2.980364408 2.976412236 2.982924075 0.006428779 0.413372468
GG Av Velocity (X) 1 [m/s] 0.848972306 0.848941452 0.847107029 0.850241922 0.000530668 0.122608402
GG Av Velocity (Y) 1 [m/s] -1.831540455 -1.829668194 -1.833717145 -1.825754575 0.004284877 0.253021782
GG Av Velocity (Z) 1 [m/s] -0.000824482 -0.000167661 -0.001516374 0.002139977 0.000272251 0.000318956
GG Force 1 [N] 10024.53419 10023.39884 10016.16723 10031.47623 5.456673634 258.3593234
GG Force (X) 1 [N] 9830.620192 9831.611848 9827.260938 9838.966932 5.108255337 324.1684822
GG Force (Y) 1 [N] 1962.189934 1951.365837 1933.038459 1967.398002 15.28155891 250.7770771
GG Force (Z) 1 [N] -1.724710333 -0.298215947 -3.570615499 2.47369766 0.409218665 1.560905132
GG Torque (X) 1 [N*m] -589.4961406 -586.1234864 -591.5626536 -580.724276 4.875404093 74.36773774
GG Torque (Y) 1 [N*m] 2949.186628 2949.328613 2947.791363 2952.381923 1.648572287 97.09764009
GG Torque (Z) 1 [N*m] 790.1034493 789.6889734 787.4732356 792.2394864 1.024732688 75.98857233
SG Normal Force 1 [N] 7601.650011 7610.28042 7600.146935 7622.455858 22.30892333 421.2082635
TRABAJO [J] 7.88134E+12 7.89926E+12 7.87823E+12 7.92454E+12 46318235797 7.19784E+11

Iterations [ ]: 73
Analysis interval: 21
 b) PARA EL ALABE DIRECTRIZ UBICADO A 15°

MICHELLE BAN.SLDPRT.SLDASM [ANALISIS DE TURBINA M.BANKI(1) [15]]

Goal Name Unit Value Averaged Value Minimum Value Maximum Value Delta Criteria
GG Av Temperature (Fluid) 1 [K] 293.2176807 293.2171951 293.2150606 293.2176807 0.00262012 0.002689071
GG Mass Flow Rate 1 [kg/s] 0.007747301 0.0091271 0.003229382 0.016238869 0.003970195 3.862472954
GG Av Velocity 1 [m/s] 2.199477755 2.201635945 2.198560154 2.206742633 0.008182478 0.304638113
GG Av Velocity (X) 1 [m/s] 0.545203223 0.544804475 0.543487919 0.546427039 0.00293912 0.084054134
GG Av Velocity (Y) 1 [m/s] -1.267098286 -1.268605208 -1.271774843 -1.266402031 0.005372812 0.175223685
GG Av Velocity (Z) 1 [m/s] -0.000666746 -0.000133966 -0.000917782 0.000350608 7.11529E-05 0.000121252
GG Force 1 [N] 10107.57272 10106.15797 10103.31019 10111.38397 8.073781324 222.9118719
GG Force (X) 1 [N] 10055.11229 10052.91096 10050.46448 10057.92906 7.464579974 292.8338944
GG Force (Y) 1 [N] 1028.466367 1036.041546 1027.664963 1044.835485 17.07346874 306.687151
GG Force (Z) 1 [N] -0.015580658 -0.053035739 -0.899544514 0.784322895 0.613322197 1.15990987
GG Torque (X) 1 [N*m] -308.957535 -311.4239043 -313.9513579 -308.6842204 5.074590723 91.34289922
GG Torque (Y) 1 [N*m] 3016.782516 3016.058608 3015.180073 3017.447139 2.267065843 87.82718135
GG Torque (Z) 1 [N*m] 621.2072146 620.6431901 619.224256 622.1017572 2.877501163 84.48138479
SG Normal Force 1 [N] 4296.015446 4294.703768 4288.429757 4300.466544 12.03678648 268.8151682
TRABAJO [J] 2.51719E+12 2.51566E+12 2.50831E+12 2.52241E+12 14100431370 2.83131E+11

Iterations [ ]: 54
Analysis interval: 21
 c) PARA EL ALABE DIRECTRIZ CERRADO

MICHELLE BAN.SLDPRT.SLDASM [ANALISIS DE TURBINA M.BANKI(2) [CERRADO]]

Goal Name Unit Value Averaged Value Minimum Value Maximum Value Delta Criteria
GG Av Temperature (Fluid) 1 [K] 293.2188864 293.2188801 293.2188564 293.2188868 3.04122E-05 0.00198691
GG Mass Flow Rate 1 [kg/s] 0.000877869 0.000708353 -7.5089E-05 0.001527475 0.000975657 0.967511351
GG Av Velocity 1 [m/s] 0.673152066 0.673969485 0.67256698 0.675312783 0.002029168 0.086831982
GG Av Velocity (X) 1 [m/s] 0.133862658 0.133582036 0.133319438 0.133862658 0.000374084 0.018761584
GG Av Velocity (Y) 1 [m/s] -0.332738463 -0.332524455 -0.332944292 -0.332162258 0.000280246 0.040109525
GG Av Velocity (Z) 1 [m/s] -0.000365485 -0.000198112 -0.000447768 3.66963E-05 4.5169E-05 4.583E-05
GG Force 1 [N] 9951.59707 9951.218176 9950.605444 9952.016087 0.321787118 97.2885241
GG Force (X) 1 [N] 9951.385142 9951.00617 9950.397021 9951.804771 0.334340602 120.2994545
GG Force (Y) 1 [N] 64.94546422 64.95452668 63.97747661 66.03580977 1.892426455 243.0640735
GG Force (Z) 1 [N] -0.284801902 -0.076378075 -0.339486509 0.152840383 0.048218641 0.65643204
GG Torque (X) 1 [N*m] -19.6601218 -19.58569509 -19.95683293 -19.19491013 0.584112221 72.70455298
GG Torque (Y) 1 [N*m] 2985.512129 2985.33417 2985.099036 2985.628281 0.10322985 36.0472856
GG Torque (Z) 1 [N*m] 433.5320295 433.4725444 433.1923821 433.7187804 0.472304412 69.27743205
SG Normal Force 1 [N] 526.2462113 524.6535491 523.3083552 526.2462113 2.937856129 323.3914488
TRABAJO [J] 37771344780 37543165226 37350792463 37771344780 420552316.5 1.23914E+11

Iterations [ ]: 88
Analysis interval: 21
VIII. CONCLUSIONES

 Debido al aprendizaje del uso del programa flow simulation, podemos

comprobar ,al realizar esta simulación y al saber la parte teórica, que
estamos capacitados para analizar de forma real cualquier máquina fluido
dinámica y poder diseñarlas.

 El análisis de la turbina que se ha hecho es de forma bastante general,

pero con la programación y la máquina ya hechas se puede calcular
valores más precisos y puntuales para algún tipo de aplicación que se
desee.

IX. BIBLIOGRAFÍA
 file:///C:/Users/USER/Downloads/Practica%20Supervisada%20%20Andres
%20Alejandro%20Leal.pdf
 http://www.biblioteca.usac.edu.gt/tesis/08/08_0550_M.pdf
 http://tesis.ipn.mx/dspace/bitstream/123456789/98/1/TESIS%20APLICACIoN%20.
pdf
 http://www.uned.es/ribim/volumenes/Vol17N1Abril2013/V17N1A06%20Marchegiani
.pdf
 Marchegiani, A. R., Nigro, N. M., & Stortti, M. A. [2002]. Modelación numerica del
flujo en el inyector de una turbina hidraulica de flujo Transversal
X. ANEXOS

Componentes en solidworks
Exterior o carcasa

Rodete
Álabe

Ensamble
 Ubicación del álabe
 a 0°

 a 15°
 Cerrado