Modelamiento Hidráulico Rev 1 101218

PROYECTO CONSTRUCTIVO DE LA PLANTA
DE TRATAMIENTO DE AGUAS RESIDUALES

TABOADA - PTAR TABOADA
MODELAMIENTO - HIDRAULICO
PLANTA DE TRATAMIENTO DE
AGUAS RESIDUALES “LA TABOADA”
CALLAO – LIMA-PERU
Noviembre 2010
ANEXO 2
MODELAMIENTO HIDRÁULICO
REALIZADO POR EMPRESA KSB
PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 1 de

27
CAMARA DE BOMBEO DE AGUA CRUDA

CBD 1
Índice
2.3 Base de cálculo ............................................................................................. 5
2.1 Dibujos disponibles..................................................................................... 20
2.2 Bombas............................................................................................................ 20
2.3 Base de cálculo ........................................................................................... 21
4. Resumen ............................................................................................................ 26

27
1. Introducción
La planta de aguas servidas PTAR Taboada consta de dos estaciones de bombeo las cuales
deben ser equipadas con bombas axiales tipo Amacan P. Para estudiar las condiciones de
flujo a las bombas, se llevarán a cabo cálculos computacionales de dinámica de fluidos. La
meta era realizar un análisis, de las distribuciones de velocidad a nivel de suelo y de las
condiciones de flujo a las bombas. La estación de agua bruta está equipada con 9 bombas
tipo AmacanP 1200-870.
El plano disponible de la planta consideraba una gran profundidad en la sentina. Dadas las
pequeñas velocidades que se producen con esto, este diseño tiende a la sedimentación.
Para evitar la sedimentación el nivel de instalación de las bombas será elevado en forma
importante.
El riesgo de la aparición de vórtices es eliminado mediante el uso de la cámara cubierta de

KSB.
Para poder lograr la atención del cliente, los cálculos computacionales de dinámica de
fluidos necesarios debían ser llevados a cabo in situ. En la preparación de estos trabajos se
modeló la sentina modificada. La Fig.1 muestra el modelo completo de la estación AGUA
CRUDA.
Fig.1: Modelo completo de la estación AGUA BRUTA
Para el cálculo in situ, el modelo fue modificado debido a la capacidad limitada de memoria
del equipo computacional utilizado.

27
En el cálculo con menos de 5 bombas el modelo es completo considerando las bombas en

operación.
En el caso de las 8 bombas operando, con el fin de disminuir los requerimientos del
modelo, las bombas se reemplazaron por una condición de borde.
Esta condición de borde era posterior al rectificador de flujo de la cámara de succión de las
bombas. Para este cálculo se analizó especialmente la distribución de velocidad en la
antecámara.
Con la utilización del rectificador de flujo y la cámara de succión cubierta, las bombas son
desacopladas de la antecámara. Dado esto, este procedimiento de simplificación es
admisible. Por lo cual el modelo fue descompuesto en las siguientes variantes:
bombas 5 y 7
velocidad de entrada
bombas 2 a la 8
velocidad de entrada
a) Cálculo con 2 bombas en operación (Modelo I) b) Cálculo con 8 bombas en operación

(Modelo II)
Fig. 2: Modelos estación AGUA CRUDA
Mientras que en la variante a) se evaluaron los flujos en la antecámara y en la cámara

cubierta para el nivel mínimo, en la variante b) se analizaron principalmente los flujos para
el nivel máximo.
2. Datos técnicos

27
2.1 Dibujos disponibles
En la planta de la estación de AGUA BRUTA, en azul se indica el límite de la zona de

cálculo.
El nivel de piso de la sentina en la zona de las bombas se fijo en h=-0.2m.
2.2 Bombas:
9 AmacanP A4 1200 870/19012UTG1
Cantidad: 9
Punto de Operación
Caudal / Capacidad: 2538 l/s
Altura de elevación (man.): 5 m
Potencia: 156.8 kW
Velocidad: 588 min-1
2.3 Base de cálculo

Los cálculos realizados son resumidos en el concepto CFD (Computational Fluid Dynamics).
Bajo esto se entiende la determinación de flujos, intercambio de calor y materia, reacciones
químicas y fenómenos similares mediante soluciones numéricas de ecuaciones
matemáticas. Estas se obtienen de ecuaciones de, conservación, continuidad, movimiento,
movimiento rotatorio y energía.

27
En particular son:
 Ecuación de continuidad
 
 ui   0
t xi
 Ecuación de la conservación del movimiento
 T
ui    ui u j   ij  f i
t xi x j
 Ecuación de la conservación de la energía
u

e   ui e  Tij j  hi  q
t xi xi xi
Donde  es la densidad del fluido, u i la componente de velocidad, Tij la componente de

tensión, f i la componente de fuerza por unidad de masa, e la energía interna específica,
h el vector del flujo de calor por unidad de área y q la fuente de calor.

27
Las tres ecuaciones se resumen como las ecuaciones de Navier-Stokes y que con la elección
de condiciones de borde adecuadas se pueden acoplar en un sistema no lineal y parcial de
ecuaciones diferenciales que pueden ser resueltas numéricamente. Para flujos
incompresibles, isotermales en un volumen se obtienen consecuentemente cuatro
ecuaciones con cuatro incógnitas (tres componentes de velocidad y la presión) que pueden
ser determinadas.
Para el cálculo de flujos turbulentos se utilizan los llamados modelos turbulentos. En los
cuales para el estudio se utiliza el modelo k    Modell de viscosidad. La viscosidad se
calcula de:
k²  u iu  i
t  C  , 
  x j x j
Donde C  es una constante empírica, k la energía cinética y  la tasa de disipación.
Para poder utilizar las descritas ecuaciones de Reynolds – Navier-Stokes en flujos, estos
deben estar divididos en una finita cantidad volúmenes de control (método de volúmenes
finitos). Esto nuevamente debe ser elegido tal que, para cada uno de casos, se obtenga el
óptimo entre exactitud y tiempo de cálculo.
Para resolver el cálculo CFD, se utiliza el programa FLUENT 6.3 de la empresa ANSYS Inc.
Parámetros seleccionados
Cantidad de elementos de la malla de modelo base: aprox.

4.200.000
Cantidad de elementos de la malla de las variantes: aprox.

1.550.000
Modelo turbulento: RNG k    Modell
Esquema de la discrecionalidad: First Order Upwind
3. Resultado

27
3.1. Cálculo CFD para la estación AGUA BRUTA variante a)
De acuerdo al modelo base se definieron los siguientes parámetros
Ancho final del canal de entrada B = 7.120 m
Nivel del canal de entrada h = 0.74 m
Nivel de la cámara de la bomba h = -0.2 m
Nivel mínimo del agua h = 1.694 (1.7) m
Cantidad de bombas en operación 2
Modelo I
La transición entre el canal de entrada (no modelado) y la antecámara es mediante un

escalón. El arrastre de los sedimentos es mejorado con un plano inclinado de 10°.
Plano inclinado10°
Fig. 3: Modelo I con 2 bombas en operación
Para el análisis de las velocidades se consideraron las siguientes niveles o planos:
1. Nivel a 0.1 m del fondo de la cámara para describir el flujo cercano al fondo.
2. Corte en la cámara en el nivel h=0.74 m para describir el flujo base.
3. Flujo en la superficie de la cámara para identificar zonas de vórtices.
Los resultados se muestran en la Fig. 4-9.

27
La distribución de velocidades cercana al fondo de la cámara (Fig.4) muestra zonas de baja

velocidad a la izquierda y derecha de la cámara principal. Estas zonas son utilizadas para
permitir un rebalse desde el lado de descarga hacia la antecámara.
Fig. 4: Distribución de velocidad en el fondo de la cámara – 2 bombas en operación - nivel

mínimo del agua
En la zona de ingreso, la distribución de velocidades en el fondo hace notar el flujo inverso

debido al escalón (Fig.5). Las regiones más alejadas del ingreso ya no son alcanzadas por los
flujos del canal. Por lo cual aquí también aparecen zonas de baja velocidad. El flujo en la
zona de los rectificadores de flujo y en la cámara de las bombas está influenciado por la
succión de las bombas.

27
Fig5: Flujo inverso en la zona del escalón de entrada.
Fig. 6: Distribución de velocidad a la altura del fondo del canal – 2 bombas en operación
El flujo a la altura del fondo del canal se ve cualitativamente similar. La influencia del flujo
inverso en el escalón de entrada no se aprecia (Fig. 6).

27
Las velocidades más altas se observan en la superficie (Fig. 7).
Fig. 7: Distribución de velocidad en la superficie – 2 bombas en operación

27
Fig. 8: Líneas de corriente en la superficie – 2 bombas en operación
Las líneas de corriente mostradas en la Fig. 8 indican solo en zonas de baja velocidad la
tendencia a formar vórtices. Debido al bajo contenido de energía, estas zonas no son
relevantes. El carácter asimétrico del flujo posterior a los rectificadores, dado la utilización
de cámara cubierta, no tiene importancia en la operación de l a bomba. Las líneas de
corriente en la cámara Fig. 9 y la distribución de velocidades en la sección transversal de la
entrada de la bomba Fig. 10 tienen características simétricas.
Fig. 9.: Líneas de corriente en la zona de la cámara.

27
Fig. 10.: Distribución de velocidad en la entrada de la bomba
Las líneas de corriente calculadas en la zona cercana a la cámara muestran un flujo regular.
3.2. Cálculo CFD para la estación AGUA BRUTA variante b)
De acuerdo al modelo base se definieron los siguientes parámetros
Ancho final del canal de entrada B = 7.120 m
Nivel del canal de entrada h = 0.74 m
Nivel de la cámara de la bomba h = -0.2 m
Nivel máximo del agua h = 2.74 (2.7) m
Cantidad de bombas en operación 8
Caudal 8 x 2.53 = 20.24 m3/s
Modelo II
Con las 8 bombas en operación y nivel máximo de agua, en gran parte de la estación no hay
zonas de baja velocidad. Solo las últimas bombas 8 y 9 no son alcanzadas por corriente del
canal con lo que los niveles de velocidad disminuyen fuertemente. En la zona de los

27
rectificadores y frente a las bombas la distribución de velocidad se influencia por el efecto

de succión de las bombas y no por el flujo de entrada a la estación.
Fig. 11: Distribución de velocidades en el fondo – 8 bombas en operación

27
Fig. 12: Distribución de velocidad a la altura del fondo del canal – nivel máximo del agua – 2
bombas en operación
Fig. 13: Distribución de velocidad en la superficie – nivel máximo del agua – 8 bombas en
operación
La distribución de velocidades a la altura del fondo del canal y en la superficie muestran

cualitativamente características similares a la variante a). Las velocidades más altas se
observan lógicamente en la superficie.

27
Fig. 14.: Líneas de corriente desde la entrada – 8 bombas en operación
Las líneas de flujo superficiales mostradas en la Fig. 14 muestran solo en la zona cámaras
sin flujo una tendencia a formar vórtices. Debido al bajo nivel de velocidad, estas zonas no
son relevantes.
El flujo hacia la cámara de la bomba muestra un carácter regular.
4. Resumen
Los resultados se pueden resumir como sigue:
Al operar con 2 bombas y como era esperable, en la zona de la antecámara aparecen

regiones de media/baja velocidad, en torno a 0,4 m/s. En el flujo directo hacia las bombas,
debido a la cámara cubierta, el flujo de la antecámara prácticamente no tiene influencia
alguna. La forma de la cámara posibilita un bajo nivel de agua debido a un reducido
requerimiento de sumergencia. Bloqueos en el ingreso son compensados por la
aceleración.
En estas regiones, aún con la operación de las 8 bombas y a carga completa de la estación,
no habrá flujo y por lo tanto deberían limpiadas con un flujo inverso forzado a través del
rebalse del lado de descarga. Para el caudal máximo de la estación el flujo no muestra nada
anormal.
La zona muerta producto del escalón de entrada podría ser evitada mediante la ejecución
de una rampa de 10° al final del canal.
27
Debido a la cámara cubierta, en las dos condiciones de operación, no se observan

condiciones indeseadas para las bombas.
Springer
KSB Aktiengesellschaft
Frankenthal, Alemania

27
CAMARA DE BOMBEO DE AGUA A EMISOR

CBD 2
Índice
1. Introducción ......................................................................¡Error! Marcador no definido.18
2. Datos técnicos................................................................................................................. 19
2.1 dibujos disponibles .....................................................¡Error! Marcador no definido.19
2.2 Bombas .....................................................................¡Error! Marcador no definido.19
2.3 Base de cálculo........................................................................................................ 20
3. Resultado...........................................................................¡Error! Marcador no definido.21
4. Resumen......................................................................................................................... 25

27
1. Introducción.
La ESTACION DE BOMBEO DE AGUA A EMISARIO de PTAR Taboada se utiliza como estación

de bombeo de salida.
Esta estación consta de 8 bombas Amacan PA4 1400-870 instaladas en forma simétrica.
Para el estudio de las condiciones de flujo se utilizaron calculo CFD.
En la estación se utilizo un diseño estandarizado en ingeniería de aguas servidas.
El ancho de cámara utilizado B = 2800 mm corresponde al ancho estándar de KSB.
La ESTACION DE BOMBEO DE AGUA A EMISARIO está diseñada en forma simétrica. Por esta
razón se debe modelar solo una mitad de la estación. Las diferencias en los flujos de
entrada son consecuencia de la curva del canal de entrada. Para la operación de la cámara
de la bomba, dada la separación hidráulica de la cámara y la antecámara, estas diferencias
son de una relevancia despreciable. La Fig. 1 muestra el modelo utilizado.
Fig.1: Modelo-CFD de la ESTACION DE BOMBEO DE AGUA A EMISARIO

27
2. Datos técnicos
2.1 Dibujos disponibles
Planos de la ESTACION DE BOMBEO DE AGUA A EMISARIO – con línea azul se indica el límite
del cálculo.
2.2 Bombas
8 Amacan PA4 1400-870/H29012UTG2
Cantidad de bombas: 8
Punto de operación
Caudal / Capacidad: 2900 l/s

Altura de elevación (man.): 7,00 m
27
Velocidad: 585 rpm

Nivel de agua en la antecámara 4.7 m
2.3 Base de cálculo
Los cálculos realizados son resumidos en el concepto CFD (Computational Fluid Dynamics).
Bajo esto se entiende la determinación de flujos, intercambio de calor y materia, reacciones
químicas y fenómenos similares mediante soluciones numéricas de ecuaciones
matemáticas. Estas se obtienen de ecuaciones de, conservación, continuidad, movimiento,
movimiento rotatorio y energía.
En particular son:
Ecuación de continuidad
 
 ui   0
t xi
 Ecuación de la conservación del movimiento
 T
ui    ui u j   ij  f i
t xi x j
 Ecuación de la conservación de la energía
u

e   ui e  Tij j  hi  q
t xi xi xi

27
Donde  es la densidad del fluido, u i la componente de velocidad, Tij la componente de

tensión, f i la componente de fuerza por unidad de masa, e la energía interna específica,
h el vector del flujo de calor por unidad de área y q la fuente de calor.
Las tres ecuaciones se resumen como las ecuaciones de Navier-Stokes y que con la elección
de condiciones de borde adecuadas se pueden acoplar en un sistema no lineal y parcial de
ecuaciones diferenciales que pueden ser resueltas numéricamente. Para flujos
incompresibles, isotermales en un volumen se obtienen consecuentemente cuatro
ecuaciones con cuatro incógnitas (tres componentes de velocidad y la presión) que pueden
ser determinadas.
Para el cálculo de flujos turbulentos se utilizan los llamados modelos turbulentos. En los
cuales para el estudio se utiliza el modelo k    Modell de viscosidad. La viscosidad se
calcula de:
k²  u iu  i
t  C  , 
  x j x j
Donde C  es una constante empírica, k la energía cinética y  la tasa de disipación.
Para poder utilizar las descritas ecuaciones de Reynolds – Navier-Stokes en flujos, estos
deben estar divididos en una finita cantidad volúmenes de control (método de volúmenes
finitos). Esto nuevamente debe ser elegido tal que, para cada uno de casos, se obtenga el
óptimo entre exactitud y tiempo de cálculo.
Para resolver el cálculo CFD, se utiliza el programa FLUENT 6.3 de la empresa ANSYS Inc.
Parámetros seleccionados
Cantidad de elementos de la malla de modelo base: aprox. 1.850.000
Modelo turbulento: RNG k    Modell
Esquema de la discrecionalidad: First Order Upwind

27
3. Resultados
3.1. Cálculo CFD para la ESTACION DE BOMBEO DE AGUA A EMISARIO
La Fig. 2 muestra la característica del flujo. La alta velocidad al nivel del suelo de aprox. 2
m/s corresponde al diseño de cámara utilizado. Si se conserva el ancho de 2.8 m entonces
se debería aumentar la profundidad de la abertura de los actuales 0.8 m a 0.9 m.
Fig. 2: Distribución de velocidad a una altura de 0.1 m del fondo de la cámara
En las Fig. 3 y 4 se analizaron las líneas de corriente desde la abertura de entrada a la

cámara. Ambas figuras tienen la misma calidad pero diferentes largos.
De la Fig.4 se deduce que parte del flujo cercano al fondo pasa de largo por la bomba
aumentando la turbulencia en la cámara. Si se desea evitar esto, se debería utilizar una
corta cámara cubierta. En la Fig.5 se muestran análisis de las líneas de corriente desde la
27
entrada de la bomba hacia el flujo de entrada. También aquí se puede reconocer que una
parte del flujo circulante de la cámara es tomado por la bomba. Debido a la pequeña
proporción, esta característica no representa un gran problema para el flujo de entrada a la
bomba. El flujo se muestra mayormente estable e independiente del tiempo.
Fig. 3: Líneas de corriente desde la abertura de la antecámara

27
Fig. 4: Líneas de corriente desde la abertura de la antecámara – largo mayor
Fig. 5: Líneas de corriente desde las bombas

27
Zuströmrichtung
Fig. 6: Distribución de velocidad a nivel del impulsor
La distribución de velocidad representada en la Fig.6 considera también el nivel del

impulsor. Sin embargo se puede reconocer la dirección del flujo de alimentación. La
distribución de velocidad obtenida tiene la característica esperada.
4. Resumen
El diseño de la ESTACION DE BOMBEO DE AGUA A EMISARIO considera una configuración

estándar. Las cámaras de las bombas están separadas mediante una antecámara de los
efectos indeseados del flujo de entrada.
El ancho de la cámara de B = 2800 mm corresponde al valor documentado para la Amacan

PA4 1400-870. La distancia entre la bomba y la abertura de la antecámara, es decir la razón
A/B = 11.5/2.8 es suficientemente grande.
Las zonas de altas velocidades presentes cerca del fondo son completamente admitidas por
las bombas. A través de esto se revela una turbulencia reducida en la cámara.
Sería recomendable rellenar las esquinas de la cámara con el fin de evitar zonas de aguas
quietas.

27
Springer
KSB Aktiengesellschaft
Frankenthal, Alemania

27

Modelamiento Hidráulico Rev 1 101218

Cargado por

Copyright:

Formatos disponibles

Modelamiento Hidráulico Rev 1 101218

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Modelamiento Hidráulico Rev 1 101218

Cargado por

Copyright:

Formatos disponibles

PROYECTO CONSTRUCTIVO DE LA PLANTA

DE TRATAMIENTO DE AGUAS RESIDUALES

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 1 de

CAMARA DE BOMBEO DE AGUA CRUDA

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 2 de

El riesgo de la aparición de vórtices es eliminado mediante el uso de la cámara cubierta de

Fig.1: Modelo completo de la estación AGUA BRUTA

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 3 de

En el cálculo con menos de 5 bombas el modelo es completo considerando las bombas en

a) Cálculo con 2 bombas en operación (Modelo I) b) Cálculo con 8 bombas en operación

Fig. 2: Modelos estación AGUA CRUDA

Mientras que en la variante a) se evaluaron los flujos en la antecámara y en la cámara

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 4 de

2.1 Dibujos disponibles

En la planta de la estación de AGUA BRUTA, en azul se indica el límite de la zona de

El nivel de piso de la sentina en la zona de las bombas se fijo en h=-0.2m.

9 AmacanP A4 1200 870/19012UTG1

2.3 Base de cálculo

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 5 de

 Ecuación de la conservación del movimiento

 Ecuación de la conservación de la energía

Donde  es la densidad del fluido, u i la componente de velocidad, Tij la componente de

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 6 de

Donde C  es una constante empírica, k la energía cinética y  la tasa de disipación.

Cantidad de elementos de la malla de modelo base: aprox.

Cantidad de elementos de la malla de las variantes: aprox.

Modelo turbulento: RNG k    Modell

Esquema de la discrecionalidad: First Order Upwind

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 7 de

3.1. Cálculo CFD para la estación AGUA BRUTA variante a)

De acuerdo al modelo base se definieron los siguientes parámetros

Ancho final del canal de entrada B = 7.120 m

Nivel del canal de entrada h = 0.74 m

Nivel de la cámara de la bomba h = -0.2 m

Nivel mínimo del agua h = 1.694 (1.7) m

Cantidad de bombas en operación 2

La transición entre el canal de entrada (no modelado) y la antecámara es mediante un

Fig. 3: Modelo I con 2 bombas en operación

Para el análisis de las velocidades se consideraron las siguientes niveles o planos:

Los resultados se muestran en la Fig. 4-9.

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 8 de

La distribución de velocidades cercana al fondo de la cámara (Fig.4) muestra zonas de baja

Fig. 4: Distribución de velocidad en el fondo de la cámara – 2 bombas en operación - nivel

En la zona de ingreso, la distribución de velocidades en el fondo hace notar el flujo inverso

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 9 de

Fig5: Flujo inverso en la zona del escalón de entrada.

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 10 de

Las velocidades más altas se observan en la superficie (Fig. 7).

Fig. 7: Distribución de velocidad en la superficie – 2 bombas en operación

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 11 de

Fig. 8: Líneas de corriente en la superficie – 2 bombas en operación

Fig. 9.: Líneas de corriente en la zona de la cámara.

PROYECTO CONSTRUCTIVO DE LA PLANTA DE TRATAMIENTO DE AGUAS RESIDUALES TABOADA - PTAR TABOADAPágina 12 de

Fig. 10.: Distribución de velocidad en la entrada de la bomba

3.2. Cálculo CFD para la estación AGUA BRUTA variante b)

De acuerdo al modelo base se definieron los siguientes parámetros

Ancho final del canal de entrada B = 7.120 m

Nivel del canal de entrada h = 0.74 m

Nivel de la cámara de la bomba h = -0.2 m