UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

Facultad de Ciencias Físicas, E.P. Ingeniería Mecánica de Fluidos,

Cl. Germán Amézaga N° 375 – Lima 01, Telf. 619-7000 Anexo 3819, e-mail:

SIMULACIÓN FLUIDODINÁMICA EN UNA BOMBA

CENTRÌGUGA Y CÀLCULO DE LA CURVA ALTURA VERSUS
CAUDAL

AUTORES
Alfaro Vidal Lenin Jharol
Arèvalo Cotrina Aaron
Sayas Palacios Josè Marìa
Yanac Calvo Paolo

RESUMEN
Con la finalidad de estudiar las interacciones y propiedades físicas del agua dentro de una bomba centrífuga, en este
trabajo se lleva a cabo el análisis numérico del campo de flujo a través del dominio fluidizado impulsor-voluta en
una bomba centrifuga; por medio del cual se valida la confiabilidad del método de dinámica de fluidos
computacional (CFD), por la aproximación de los resultados arrojados, lo cual es de gran ayuda para analizar sus
características funcionales mediante el uso de computadoras, permitiendo afinar el diseño geométrico de la bomba
mediante la ubicación de las zonas con presencia de fenómenos como recirculación o cavitación, que pueden afectar
de forma permanente el impulsor de la bomba y otros elementos constitutivos, como procurar el aumento de su
eficiencia global implicando un gran ahorro económico ya que se evita la necesidad de construir modelos a escalas
que cumplan con las leyes de semejanzas de las bombas.

Palabras clave: Bomba centrífuga, CFD, Métodos numéricos, Ensayos, Curvas características, Presión, Flujo.

ABSTRACT

In order to study the interactions and physical properties of water within a centrifugal pump, in this work the
numerical analysis of the flow field is carried out through the fluidized impeller-volute domain in a centrifugal
pump; through which the reliability of the computational fluid dynamics (CFD) method is validated, by the
approximation of the results obtained, which is of great help to analyze its functional characteristics through the use
of computers, allowing to refine the geometric design of the pump by locating the areas with the presence of
phenomena such as recirculation or cavitation, which can permanently affect the impeller of the pump and other
constituent elements, such as ensuring the increase of its overall efficiency involving a great economic saving since it
avoids the need to build scale models that comply with bomb similarity laws.

Keywords: Centrifugal pump, CFD, Numerical methods, Tests, Characteristic curves, Pressure, Flow.
 INTRODUCCIÓN de los fenómenos de incremento de energía y de
difusión.
La mecánica de fluidos es una ciencia que estudia
el comportamiento de los gases y líquidos así como Uno de los campos más beneficiado por la mejora
las fuerzas que los provocan. Al contrario de las de los métodos de CFD es precisamente el de
demás ciencias y sus ramas, es relativamente diseño de turbo-maquinaria, seguido de cerca por
nueva. Los primeros principios básicos del la aerodinámica. Se podría decir que estos son los
movimiento de los fluidos fueron desarrollados campos que más utilizan la mecánica
durante los siglos XVI y XIX. Nombres como computacional en sus diseños y pruebas, pero no
Pascal, Bernoulli, Euler, Navier, Kelvin o son los únicos, puesto que poco a poco va ganando
Reynolds desarrollaron los principios de la ciencia, popularidad en campos como bioingeniería o
que continuaron otros como Hagen, Pouseuille o ingeniería de construcción
Darcy durante los últimos siglos.
El trabajo de estos científicos fue el comienzo de la Debido a la complejidad de las partes que
Mecánica de Fluidos computacional componen una turbo-máquina, no se puede exponer
(Computational Fluids Dynamics, CFD). Los a cometer errores en su fabricación ni en su diseño,
primeros métodos fueron desarrollados durante el con lo que la mecánica computacional es un
año 1930, cuando los ordenadores comenzaron a elemento que permite mejorar y comprobar el
mejorarse de forma que pudieran resolver simples funcionamiento correcto de una parte de la
ecuaciones lineales. A medida que iba avanzando máquina, e incluso de detectar errores de diseño de
la tecnología lo hacía también la resolución de éstas.
problemas con CFD hasta el punto que el primer
problema resuelto se llevó a cabo en 1967. Los
estudios realizados por la NASA, hicieron que
poco a poco se llegara a la tecnología y capacidad
de resolución que se tiene hoy en día.
Los principios de todos los problemas resueltos con PROBLEMÁTICA
CFD están basados fundamentalmente en las
ecuaciones de Navier y Euler, con interfaz gráfica
o sin ella, esto depende del tipo de programa que se
esté utilizando. OBJETIVOS
La turbomaquinaria que es una parte aplicativa de
la mecánica de fluidos, clasifica como
turbomaquinas hidráulicas a las maquinas que Objetivo Generales
funcionan con agua en estado líquido, dentro de
estas están las bombas y las turbinas, el presente 1. Encontrar la curva de la bomba H vs
trabajo consiste en realizar la simulación numérica Q.
de una bomba centrifuga. 2. Aprender a utilizar la herramienta
Las bombas centrıfugas, que en cualquiera de sus computacional (ANSYS).
tipos (axial, radial, de flujo mixto), son parte de
numerosas redes hidráulicas ya que aportan la Objetivos específicos
energía necesaria al fluido para su desplazamiento
y uso de acuerdo a la aplicación a la que este 1. Interpretar la curva característica de
destinada. una bomba centrifuga.
Las propiedades del fluido que se desplaza en el 2. Saber tomar puntos por medio de una
interior de una bomba centrıfuga pueden variar simulación con ANSYS.
significativamente, conduciendo a la generación de 3. comparación de los resultados
fenómenos físicos como por ejemplo la cavitación, experimentales de una bomba
recirculación, desprendimiento de la capa limite. A centrífuga con los obtenidos mediante
través de los ensayos experimentales en bancos de un análisis computacional.
pruebas de bombas centrıfugas, se puede 4.
determinar sus curvas características y parámetros
principales, además de encontrar la eficiencia de
trabajo. Así también fenómenos básicos presentes
en las bombas centrıfugas, como la visualización
 MARCO TEÓRICO limpias, es decir, geometrías reparadas y
simplificadas. Al usar este algoritmo, el
Como se sabe el uso de una bomba conlleva a la mallado se realiza a partir de las caras y
previa observación y conocimiento de sus partes y
bordes, finalizando en los volúmenes grandes.
funcionalidades.
En términos coloquiales, al usar este algoritmo
Con la bomba a trabajar es una centrifuga la cual
posee un conjunto de paletas que se encontraran se consigue que los elementos finitos se
rotando encerradas dentro de una voluta o coraza. adapten a la geometría, y no al revés.
El funcionamiento del flujo dentro de nuestra  Patch Independent: En este algoritmo, en
bomba a través del centro o ojo del rodete y nuestro contraste con el anterior descrito, primero se
fluido ganara energía a medida que los alabes o genera una malla desde el volumen central y
paletas del rodete lo transporten o aumente su hacia los bordes y caras, de tal manera que no
velocidad hacia afuera en la dirección radial. El necesariamente se respeta la geometría, en los
aumento de aceleración por parte del fluido detalles. En términos coloquiales, al usar este
producido por las paletas nos mostrara un aumento algoritmo se consigue que la geometría se
de energía de presión y cinética, gracias a la forma adapte a la forma de los elementos finitos, y no
de caracol que posee la voluta o carcaza que genera al revés.
gracias a su geometría un aumento gradual del área
que recorrerá el flujo.  Mallado hexagonal: Es un método de mallado que
La fuerza centrífuga producida depende tanto de la genera elementos hexaédricos. La malla hexaédrica
velocidad en la periferia del impulsor como de la se caracteriza por presentar un número reducido de
densidad del líquido, la energía que se aplica por
elementos, con lo cual se mejora el tiempo de
unidad de masa del líquido es independiente de la
densidad del líquido. Por tanto, en una bomba dada convergencia del modelo. Asimismo, se consigue
que funcione acierta velocidad y que maneje una reducción del error numérico con lo que se
un volumen definido de fluido, la energía que se mejora la precisión del cálculo. Estas
aplica y transfiere al fluido, es la misma para comparaciones son válidas para mallas tetraédrica
cualquier líquido sin que importe su densidad. Por y hexaédrica de primer orden (lineales). Existen
tanto, la carga o energía de la bomba se debe dos métodos de mallados hexaédricos, los cuales se
expresar en metros y es por eso por lo que se describen a continuación.
denomina "altura".
 Hexaedros dominantes: El método hexaedro
dominante se caracteriza por crear una malla
de elementos hexaédricos dominante a los
MÉTODOS Y MATERIALES elementos tetraédricos. En términos
Métodos de mallado en ANSYS coloquiales, este método de mallado crea
elementos hexaédricos en su mayoría (donde
El software ANSYS tiene diversos métodos de la geometría del modelo permite la creación de
mallado 3D, los principales se detallan a elementos hexaédricos), y crea elementos
continuación: tetraédricos (donde la geometría no permite
que se creen elementos hexaédricos).
 Mallado automático: Es un método de mallado  Mallado de barrido: El mallado de barrido
que combina el mallado tetraédrico patch genera únicamente elementos hexaédricos los
conforming con el mallado de barrido, los cuales se cuales se forman a partir de elementos
explican posteriormente. Este método identifica cuadriláteros que se forman en la superficie a
todas las partes que pueden ser malladas con el partir de la cual se barre. Posteriormente, se
método de barrido y procede a mallar de tal manera. genera una malla 3D de elementos hexaédricos
Luego, todas las partes restantes son malladas con que se barren desde la superficie inicial. Los
el método de tetraedros patch conforming. elementos de la superficie de llegada deben
tener la misma topología que los elementos de
 Mallado tetraédrico: Es un método de mallado
la superficie inicial. Cabe mencionar que
que genera elementos tetraédricos, para los cuales
circunstancialmente se pueden crear elementos
existen dos algoritmos de mallado:
prismáticos en el barrido. La siguiente figura
 Patch Conforming: Este algoritmo es muestra lo antes mencionado.
preferentemente usado para geometrías
 continuación, una breve descripción de los
principales parámetros en la configuración local:
Estos métodos combinan los diferentes tipos de
elementos finitos, de acuerdo a las necesidades del  Sizing: Es una herramienta que sirve para definir
problema. Cada método presenta ventajas y tamaños de elementos finitos para zonas
desventajas frente a los otros métodos. La elección determinadas de la malla global.
depende del usuario y del tipo de problema a
resolver. Asimismo, además de configurar el  Contact sizing: Es una herramienta que sirve para
método de mallado a utilizar y los tipos de definir tamaño de elementos para zonas de
elementos, se debe definir las configuraciones contacto.
globales y locales del mallado.
 Mapped face meshing: Es una herramienta que
Comandos globales sirve para estructurar elementos en una cara
determinada.
Las configuraciones globales de malla se utilizan
para el ajuste global en la estrategia de mallado, que  Inflation: Sirve para crear elementos por capas,
incluye todas las funciones de sizing, inflation, como se explicó anteriormente. Se configura esta
smoothing, parámetros de entrada, etc. Son muy herramienta en la configuración local para
útiles para resolver las principales características superficies específicas.
del modelo con un mínimo de insumos. A
Una vez definido todos estos parámetros, se puede
continuación, una breve descripción de los
pasar a generar la malla.
principales parámetros en la configuración global:

 Physics based settings: Se establece la física de

preferencia, es decir, si se trata de un problema
CFD o un problema de mecánica de sólidos.

 Global sizing controls: Se establecen

configuraciones globales de tamaño de elemento.

 Inflation: Inflation es una herramienta que sirve

para crear un mallado por capas cercano a una
superficie o superficies de interés. Esta herramienta
es usada principalmente para problemas de
interacción fluido-estructura donde se requiere
información en la capa límite del fluido. En este
punto, se establece la configuración global de la
herramienta en el mallado. En la siguiente figura se
aprecia con detalle lo antes mencionado.

 Patch confirming options: Se establecen

configuraciones globales sobre el algoritmo patch
conforming para mallados tetraédricos.

 Statistics: Es una herramienta que sirve para

visualizar las estadísticas de mallado obtenidas
luego de generar la malla. Posteriormente, se
explica a detalle este punto

Comandos locales

En las configuraciones locales se establece

configuraciones específicas de malla para ciertas
zonas o superficies de interés que las requiera. A
PROCEDIMIENTO DEL DESARROLLO DE  Todos los elementos finitos deben estar
LA SIMULACIÓN correctamente interconectados.
 Por ningún motivo debe existir elementos
Geometría superpuestos.
La geometría del modelo se realiza tomando en  La distribución de los elementos finitos debe
consideración el dominio de interés (parte del ser la adecuada para capturar las características
sistema físico delimitado donde el software geométricas.
resuelve las ecuaciones gobernantes). Es
importante la correcta definición del dominio de Mallado
interés, ya que únicamente sobre este dominio se El siguiente paso es generar el mallado. Con toda
realiza la simulación numérica. La definición del la configuración realizada en el paso anterior, solo
dominio depende principalmente del sistema físico queda por generar rápidamente el mallado. Es
delimitado y del objetivo del problema. El diseño importante mencionar en este punto, que se puede
CAD del modelo se realiza con ayuda de un visualizar y evaluar el tipo, la calidad y las
software CAD. ANSYS trae algunos modeladores propiedades del mallado, luego de la realización
CAD 3D tales como SpaceClaim o Design del mismo. En caso la malla no sea la adecuada, se
Modeler, los cuales son muy sofisticados. La puede refinar y dar ajustes para llegar a la calidad
elección del modelador CAD depende del dominio de malla adecuada.
del usuario. Sin embargo, para fines de la presente
tesis, se eligió usar Autodesk Inventor como Condiciones de Frontera
software CAD, ya que presenta grandes facilidades Parámetros de control de la simulación
al dibujante al modelar en 3D. No existe ninguna (condiciones de borde) Los parámetros de control
restricción al crear el modelo 3D en Autodesk son los que determinan el tipo de modelamiento,
Inventor y exportarlo a ANSYS. Al realizar el los alcances que tiene la simulación y el grado de
modelo 3D, se recomienda simplificar al máximo precisión del cálculo. En otras palabras, los
la geometría del dominio de interés. De esta parámetros de control de la simulación son las
manera, se aligera los posteriores procesos durante condiciones de frontera o borde que se aplican
la simulación como es la etapa de solución, por tanto en los límites del dominio de interés, como
ejemplo. Asimismo, una vez que el modelo 3D es dentro del mismo. Condiciones de borde más
exportado a ANSYS, es recomendable “limpiar” la utilizadas en problemas CFD
geometría del mismo, es decir, corregir errores de  Material del dominio de interés
dibujo tales como cruces de aristas, etc. Esto se  Tipo de turbulencia
realiza con el fin de reducir la complejidad de las  Condiciones de ingreso, salida, etc.
mallas y, al mismo tiempo, hacer de la geometría lo
 Velocidad inicial
más cercana a la realidad. Asimismo, la facilidad
 Temperatura inicial
en el manejo del software CAD es imprescindible
en esta etapa.  Tipo de condición de borde, etc.

Condiciones de borde más utilizadas en problemas

Malla
Luego de tener definido y creado el dominio de de mecánica de sólidos  Presión  Momentos
interés, se procede al mallado. Se define mallado  Fuerzas puntuales o distribuidas
como el proceso de dividir el dominio de interés en  Velocidad
subdominios mucho más pequeños (elementos  Contactos
finitos), donde a su vez se resolverán las ecuaciones  Apoyos
gobernantes del modelo. Sin embargo, antes de  Condiciones de temperatura
realizar el mismo, se debe definir los parámetros
del mallado, tales como tipo de mallado (geometría La correcta configuración de los parámetros de
del elemento a usar), método de mallado a utilizar, control de la simulación van de acuerdo al objetivo
etc. En general, durante el mallado se debe cumplir del problema e influyen totalmente en la
con las siguientes indicaciones: simulación. Por ello, esta etapa debe ser de sumo
 La selección de la forma del elemento cuidado. Cada parámetro elegido debe ser
(tetraédrico, hexagonal, etc.) debe ser estudiado previamente. Como se mencionó
completamente entendido por el solver del anteriormente, que la simulación converja no
ANSYS. asegura la correcta solución del problema.
 Los elementos finitos deben cubrir al dominio
de interés en su totalidad.
Para la obtención de la curva característica se
procede a simular utilizando como condición de ´ DISCUSION DE RESULTADOS
borde de entrada la presión de vapor para evitar la
cavitación, y como condición de ´ borde de salida,
el flujo másico. ´ La obtención de los puntos 1. Curva característica H vs Q:
pertenecientes a la ´ curva característica, se realiza Se presenta en la Tabla 1 todos los flujos másicos
mediante la realización de 10 simulaciones, ´ (caudales (Q)) simulados, con las presiones de
utilizando caudales desde 1m 3 /h hasta 20m 3 /h; entrada (Pe) y salida (Ps) la cual fueron elegidos al
con ayuda de los valores obtenidos del ensayo a la azar.
bomba centrífuga, los puntos simulados en este
trabajo, los cuales dieron resultados satisfactorios Después se halla muestra la curva característica,
en cuanto a la convergencia. disponibilidad en función del caudal (H vs Q), ´
obtenida mediante simulación.

Flujo màsico Presiòn de Como se puede observar en la Figura XXXX, la

Punto tendencia de la curva H vs Q, encontrada mediante
(𝒌𝒈/𝒔) Salida
1 0.1 el análisis numérico, tiene una tendencia
2 0.2 decreciente, como sugiere la teoría de bombas
3 0.25 centrífugas.
4 0.3
5 0.4 2). Visualización de los Resultados:
6 0.5
7 0.6 Ya tenemos los valores previamente escogidos
8 0.9 como valor de frontera para ingresarlos al Ansys.
9 1
10 1.33 Figura XYZ
Tabla 1: Condiciones de borde de las simulaciones
realizadas. Se observa en la Figura xyz dla presión en el tubo
de entrada es negativa porque la bomba absorbe el
El flujo másico se calcula en función del ´ caudal fluido, en este caso usamos la presión de vapor
y de la densidad, la cual es constante, por para evitar la cavitación en el rodete de la bomba.
considerar que el fluido tiene una temperatura
constante de 15◦C. Se observa que la Figura xyz la presión es mucho
Los resultados se validaran, comparando la ´ mayor, esta energía de presión se ha dado por la
tendencia de la curva característica, Altura en bomba.
función del Caudal (H vs Q) obtenida ´ mediante La energía cinética se mantiene constante porque
el análisis numérico el diámetro de entrada de la bomba y el diámetro
de salida de la bomba se mantiene constante al
igual que el caudal.

La variación de energía potencial se mantiene

invariable (es casi nulo). Sólo se ve el cambio de
Energía de presión

CONCLUSIONES Y RECOMENDACIONES

En este artículo se ha ilustrado el papel de la

Dinámica de Fluidos Computacional (CFD, por sus
siglas en inglés) en el diseño, optimización y
análisis de turbomáquinas hidráulicas.
 REFERENCIAS

Marden Marshall, Elizabeth and Bakker, André.

“Computational fluid mixing”. Fluent Inc.
Lebannon, New Hampshire, USA.
www.fluent.com

Bakker, André. “Overview balance equations

isothermal single-phase flow and turbulence
models”.
Simulación de un elemento de una turbina
hidráulica

Versteeg, H.K and Malalasekera, W. “An

introduction to Computational Fluid Dynamics:
The Finite Volume Method”. Longman Scientific
& Technical, 1995.

White, Frank M. “Mecánica de Fluidos”. McGraw-

Hill, 2000. Traducido por Rodríguez Fernández,
Manuel profesor adjunto del departamento de
mecánica de fluidos de la Universidad Politécnica
de Madrid

.Mataix, Claudio. Turbomáquinas Hidráulicas 2ª

Edición. Madrid. Publicaciones Universidad
Pontificia Comillas, 2009. (MATA09)

· Álvaro Ruiz de Galarreta López. Diseño de una

Bomba-Turbina Integral. Trabajo Fin de Grado (U.
P. Comillas), 2014. (RUIZ14)

· Cayetana Urbina Soguero. Prototipado y ensayo

de un modelo a escala de bomba-turbina integral.
Trabajo Fin de grado (U. P. Comillas), 2015.
(URBI15)

· Santiago Laín Beatove. Simulación numérica del

flujo en turbomáquinas hidráulicas. Estado del arte
y fuentes de error. REVISTA Universidad EAFIT.
Vol. 44. No. 152, 2008. (LAIN08)