Trabajo # 3

Simulación 3D de un Vertedero en Open Foam

Obras Hidráulicas.

Autores:
Jose Manuel Duque Duque
Yerik Sánchez Sánchez
Manuela Rendón Salazar

Profesor:
Juan David Montoya Velilla

Fecha de entrega: 26/06/2023

Universidad Nacional de Colombia

Facultad de Minas
Departamento de Ingeniería Civil
Semestre 2023-1
Medellín
Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

Índice

1. Introducción 2

1.1. Características del vertedero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

2. Mallado 3

3. Parámetros Open Foam 4

3.1. Definición de Parámetros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3.1.1. Consideración importante: Número de Courant [C] . . . . . . . . . . . . . . . 5

3.2. InterFoam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

3.2.1. PIMPLE Foam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

3.2.2. Efectos de la presión hidrostática. . . . . . . . . . . . . . . . . . . . . . . . . 6

4. Simulación y Resultados 6

4.1. Simulación completa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

4.2. Simulación reducida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

5. Análisis y Conclusiones 7

Universidad Nacional de Colombia 1 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

1. Introducción
Un vertedero es una estructura u obra hidráulica cuya finalidad principal es controlar el paso de un
fluido, principalmente agua, en escurrimientos superficiales, a veces permitiendo controlar el paso
de dicho fluido por medio de compuertas o dejando un flujo libre a partir de determinada altura,
siendo bastante útil para sistemas de riego o como disipador de energía. En este informe se realizará
una simulación 3D de un vertedero donde se explica el funcionamiento y los conceptos teóricos que
conlleva su diseño. También se presentan los aspectos generales necesarios para correr una simulación
utilizando Open Foam y otros Softwares como AutoCAD y SALOME que son complementarios para
este trabajo.

1.1. Características del vertedero

Se optó por diseñar un canal simple y recto donde se puedan presenciar los fenómenos de flujo
subcrítico y supercrítico, remanso y resalto hidráulico. A continuación se muestran las dimensiones
del canal y el vertedero diseñadas en AutoCAD.

Figura 1: Dimensiones canal y vertedero.

Figura 2: Dimensiones canal.

Universidad Nacional de Colombia 2 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

El canal consta de una compuerta de entrada por donde ingresará el fluido, en este caso y de nuestra
elección será agua, el agua luego se encuentra con el vertedero diseñado, y al superar la altura de
este descenderá y chocará con el escalón final del canal, donde no hay salida del agua. Esto es lo
que se planea simular, para conocer sus velocidades y el tiempo que tardaría el agua en sobrepasar
el vertedero.

2. Mallado
Existen 2 tipos de mallado:

Malla estructurada: Las celdas o elementos que discretizan el dominio son regularmente
organizados y están estructurados en una forma geométrica predefinida, como rectángulos o
hexágonos en 2D, o prismas y hexaedros en 3D. Recomendable para canales rectos.

Malla no estructurada: Las celdas o elementos de la malla no siguen una estructura

geométrica regular y pueden tener formas arbitrarias, como triángulos en 2D o tetraedros en
3D. Las celdas vecinas no están necesariamente alineadas entre sí y pueden tener tamaños y
formas variables. Recomendable para canales con pendientes.

Para el presente análisis se utiliza una malla no estructurada, se diseña con el Software SALOME
importando el diseño del canal y el vertedero desde AutoCAD, el mallado fue el siguiente:

Figura 3: Mallado del Canal y Vertedero.

Universidad Nacional de Colombia 3 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

Figura 4: Mallado del Canal y Vertedero.

3. Parámetros Open Foam

Para comenzar una simulación en Opem Foam se requieren de 3 parámetros iniciales:

Definir las condiciones de frontera.

• Presión
FoamFile version 2.0; format ascii; class volScalarField; object p_rgh; // * * * * * * *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField atmosphere type totalPressure; p0 uniform 0;
".*" type fixedFluxPressure; value uniform 0;
• Velocidad
FoamFile version 2.0; format ascii; class volVectorField; object U; // * * * * * * * * *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (1 0 0);
boundaryField inlet type flowRateInletVelocity; volumetricFlowRate constant 2;
walls type fixedValue; value uniform (0 0 0);
atmosphere type pressureInletOutletVelocity; value uniform (0 0 0);
outlet type inletOutlet; inletValue uniform (0 0 0); value $internalField;

Definir los valores constantes.

• Gravedad
FoamFile version 2.0; format ascii; class uniformDimensionedVectorField; location çons-

Universidad Nacional de Colombia 4 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

tant"; object g; // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
//
dimensions [0 1 -2 0 0 0 0]; value (0 -9.81 0);
• Propiedades de transporte
FoamFile version 2.0; format ascii; class dictionary; location çonstant"; object transport-
Properties; // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
//
phases (water air);
water transportModel Newtonian; nu 1e-06; rho 1000;
air transportModel Newtonian; nu 1.48e-05; rho 1;
sigma 0.07;
Definir las propiedades de solución.

3.1. Definición de Parámetros

Para las condiciones de frontera definimos unas variables: el volumen de agua que entrará en la
simulación, la viscosidad cinemática, el algoritmo para correr ecuaciones de flujo turbulento, la
densidad y la velocidad de flujo.

1. Velocidad = 2 m
s

2. Volumen de agua: Especificamos la fracción volumétria, el 100 % del agua entrará por la única
entrada que tiene el canal.
3. Densidad: presión atmosférica 0 y bajo el agua se calculará con la simulación.
2 2
4. Viscosidad = 0.003 ms y Uniformidad = 0.5 ms

En las constantes se encuentra la malla, la aceleración de la gravedad (9.81 sm2 ), el tipo de flujo
(Turbulento) y las fases:

Agua: µ = 1e−6 ; ρ = 1
Aire: µ = 1.48e−5 ; ρ = 1

Parámetros de corrida: Duración de 64 segundos, paso = 0.01, intervalo = 0.2 y definimos un número
de Courant = 60.

3.1.1. Consideración importante: Número de Courant [C]

Uo · ∆t
C=
∆x

Universidad Nacional de Colombia 5 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

C = N úmerodeCourant
U = V elocidaddelf lujo
∆t = T amañopasotemporal
∆x = T amañodelpasoespacial

El número de Courant para una simulación de un proyecto debe ser menor o igual a 1 para asegurar
resultados confiables, esto implica diseñar una malla mucho más detallada creando más elementos,
al igual que tener un paso de simulación menor, cosa que generaría que la simulación tarde semanas
o incluso más tiempo que no es necesario invertir para fines académicos.

3.2. InterFoam

El solver Interfoam se encarga de resolver las ecuaciones matemáticas que describen el flujo del agua
en el vertedero, para obtener soluciones aproximadas del flujo.

En esta modelación se trabajaron los siguientes métodos:

3.2.1. PIMPLE Foam

Este método combina los métodos PISO (Presión-Implicita Separada Operacionalmente) y SIMPLE
(Semi-Implicito para las Ecuaciones de Navier-Stokes).

El algoritmo PIMPLE se utiliza principalmente para problemas de flujo transitorio en los que las
variaciones temporales son significativas, pero también se tienen en cuenta las variaciones espaciales
y las relaciones entre la presión y la velocidad del fluido.

3.2.2. Efectos de la presión hidrostática.

Estabilidad de la interfaz: Ejerce una fuerza hacia abajo en el aguqa y hacia arriba en el aire.
A medida que el agua se acumula en el vertedero, la presión hidrostática aumenta con la
profundidad. Esto se refleja en la distribución de la presión en el dominio y es importante
para evaluar la carga hidrostática en las estructuras y superficies de control.
Los gradientes de presión hidrostática generan fuerzas que aceleran o desaceleran el flujo de
agua a través del vertedero. Estos efectos son cruciales para analizar el caudal y la descarga
de agua en el vertedero.

4. Simulación y Resultados
Se corre el programa en C++ y se revisa la simulación en la extensión de ParaView y los resultados
fueron los siguientes:

Universidad Nacional de Colombia 6 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

4.1. Simulación completa

https://drive.google.com/file/d/1sjI5f8_058s5vGbUsO5bjidPE2FbJlzG/view?usp=sharing

Tiempo de simulación: 10.687 segundos

Tamaño malla: 0.2 m y 385.000 elementos

4.2. Simulación reducida

https://drive.google.com/file/d/1BR-6OqVte_q-vuLd6ju_GrZcCPk7zVqU/view?usp=sharing

Tiempo de simulación: 3.437 segundos

Tamaño malla: 0.25 m y 154.000 elementos

5. Análisis y Conclusiones
El flujo inicialmente se encuentra en estado subcrítico mientras se llena el volumen, una vez
se desborda por el vertedero el flujo se vuelve supercrítico, la altura de flujo disminuye y la
velocidad aumenta para acumularse en el otro tramo del canal.

El modelo representa de manera acertada las transiciones esperadas en una obra de control
hidráulico, en el cual se espera que se presenten transiciones del régimen del flujo y por tanto
disipación de energía.

El mallado se mostró entre valores adecuados para realizar la simulación, este se refinó en las
paredes como puede verse en 4 ya que al definir la velocidad en estas igual a cero, se requiere
que los volúmenes sean pequeños para que la pantalla de agua inmóvil no obstruya el resto de
la simulación.

El mallado también se refinó en el vertedero, ya que al ser su geometría circular, se requieren

volúmenes más pequeños para lograr un contorno adecuado.

Podemos notar la gran diferencia en el tiempo de corrido y el resultado de la simulación

solamente al refinarla y hacerla 5 cm más pequeña, muestra de la sensibilidad y capacidad de
estos softwares que ahora son vitales para el diseño y construcción de obras hidráulicas.

La simulación mostró exitosamente el comportamiento del agua al correr por un vertedero un

deja percibir los fenómenos típicos que esta estructura genera en un canal.

En la simulación por fines académicos se opto por definir un numero de Courant extrema-
damente alto, esto se hizo pues no consideramos que se necesita una resolución espacial y
temporal extremadamente alta en la simulación numérica que se presenta ya que de esta no
se extraerá información que requiera alto grado de precisión además que observamos que los

Universidad Nacional de Colombia 7 Obras hidráulicas

Universidad Nacional de Colombia
Facultad de minas Obras Hidráulicas

datos son consistentes con la teoría en cuanto a los regímenes de flujo.

De ser necesario realizar una simulación con un mayor nivel de detalle, se hace indispensable
entonces realizar un refinamiento de la malla y ajustar el time step para que estos me brinden
mayor resolución espacial y temporal además que se haría indispensable validar los datos con
respecto a las ecuaciones que se encuentren en la literatura o con investigación de campo. No
obstante queda por sentado que el CFD es una muy buena herramienta para el análisis de
problemas hidráulicos.

Universidad Nacional de Colombia 8 Obras hidráulicas