Modelo Numerico de Las Caracteristiscas Del Flujo Traducido
Modelo Numerico de Las Caracteristiscas Del Flujo Traducido
Modelo Numerico de Las Caracteristiscas Del Flujo Traducido
Resultados en Ingeniería
a
Departamento de Ingeniería Civil y Ambiental, Universiti Teknologi PETRONAS, Seri Iskandar, Malasia b
Instituto de Energía Sostenible (ISE) Universiti Tenaga Nasional Selangor, Malasia
C
Centro de Investigación de Ingeniería de Ríos y Drenaje Urbano Universiti Sains Malaysia Pulau Pinang, Malasia
Palabras clave: Los hallazgos de la investigación han demostrado que la estimación precisa del rendimiento y los perfiles de la superficie del agua en las
Vertederos de Crump
estructuras hidráulicas es engorrosa y lleva mucho tiempo. Además, investigar experimentalmente las variables hidráulicas requería más
Flujo-3D
dinero para realizar las pruebas de campo y encontrar el equipo necesario. Además, las mediciones de laboratorio de los parámetros
Flujo sobre joroba
hidráulicos no son precisas en comparación con los resultados de la simulación numérica.
Simulación numérica
Por lo tanto, la simulación numérica de problemas hidráulicos como el flujo sobre la joroba utilizando la ciencia de la dinámica de fluidos
computacional (CFD) es una excelente alternativa para las pruebas de laboratorio. En este estudio, las propiedades de flujo de un vertedero
triangular de cresta afilada de dimensión 50 30 7 cm3 se modelaron numéricamente utilizando Flow 3D y se validaron mediante experimentos
de laboratorio. Los resultados mostraron que el software modeló los parámetros hidráulicos como la elevación de la superficie del flujo, el
número de Froudes, la velocidad y la profundidad del agua a través de la joroba triangular con un mayor grado de rendimiento y precisión
(variaciones mínimas entre 4% y 5%). Por lo tanto, se recomendó Flow 3D como una herramienta útil para predecir e investigar el
comportamiento del flujo aguas arriba y aguas abajo de los vertederos, así como en diferentes estructuras hidráulicas. Además, el modelado
similar con otros modelos de turbulencia dentro del software puede considerarse una gran contribución científica al campo y, por lo tanto, es
muy recomendable para futuros estudios.
Introducción por lo tanto, requiere ser resuelto utilizando técnicas numéricas precisas y simplificadas
que han sido validadas con resultados experimentales [6,7].
En canales abiertos, las estructuras hidráulicas más comunes utilizadas para las Hoy en día, la resolución de problemas de ingeniería mediante el enfoque de
mediciones de descarga de flujo son los vertederos debido a su precisión, simplicidad y simulación numérica se ha incrementado especialmente con la gran mejora en la capacidad
facilidad de diseño y construcción [1,2]. Sus principales funciones incluyen la gestión de y habilidad de los poderes de procesamiento de la computadora [8]. Sin embargo, se
los niveles de agua aguas arriba y aguas abajo, así como la medición del caudal y la requiere la validación de la solución numérica experimentalmente para garantizar que los
estabilización de los canales [3]. Además, tienen diferentes aplicaciones de ingeniería en resultados del modelado sean lo suficientemente precisos. En términos de simulación
el laboratorio y en el campo, especialmente en seguridad, por ejemplo, en la descarga de numérica de fluidos, el método de volumen finito y el método de elementos finitos se
agua sobrante de ríos y presas durante las épocas de inundaciones. Los vertederos se utilizan para resolver las ecuaciones de gobierno de fluidos, que son las ecuaciones de
pueden clasificar en diferentes tipos según su geometría y diseño. Estos incluyen continuidad y de NavierStokes, además de los modelos de turbulencia [9–12].
vertederos de cresta ancha, de cresta afilada, de cresta estrecha y en forma de conopial. La dinámica de fluidos computacional (CFD) se ha informado entre las técnicas de
Sin embargo, según la forma de la abertura, los vertederos se pueden clasificar en 3 tipos simulación numérica utilizadas para resolver problemas de flujo de fluidos [13].
principales, a saber, vertederos rectangulares, triangulares y trapezoidales [4]. Los En este campo, se han desarrollado muchos códigos fuente abiertos y comerciales para
vertederos triangulares de cresta corta o vertederos crump se encuentran entre los tipos simulaciones numéricas. Entre los más populares se encuentran Fluent, Flow-3D y Star
especiales de cresta ancha con una pendiente en la cara posterior de 1:5 y una pendiente CCMþ, mientras que el código open-FOAM es un ejemplo de código CFD de código
en la cara frontal de 1:2 [5]. En estos tipos de vertederos, la estimación precisa del abierto [14]. El software Flow-3D utiliza el método de volumen finito (FVM) para resolver
rendimiento y los perfiles de la superficie del agua a veces es engorrosa y requiere mucho las ecuaciones gobernantes y hay diferentes modelos turbulentos disponibles, incluido el
tiempo. Los problemas complejos son modelo k-ÿ de dos ecuaciones, renormalizado
* Autor correspondiente.
Dirección de correo electrónico: [email protected] (EHH Al-Qadami).
https://doi.org/10.1016/j.rineng.2019.100052 Recibido
el 2 de octubre de 2019; Recibido en forma revisada el 16 de octubre de 2019; Aceptado el 17 de octubre de
2019 2590-1230/© 2019 Los autores. Publicado por Elsevier BV Este es un artículo de acceso abierto bajo la licencia CC BY-NC-ND (http://creativecommons.org/licenses/by nc-nd/4.0/).
Machine Translated by Google
modelo de grupo (RNG), modelo kw de dos ecuaciones y modelo de energía turbulenta de una
ecuación [15–17]. Sin embargo, los modelos k-ÿ y (RNG) son las ecuaciones más avanzadas
comúnmente utilizadas para la ingeniería de fluidos [18–20]. La mayoría de los estudios previos
emplearon el uso de métodos experimentales y algunas veces con expresiones teóricas Fig. 3. Geometría y dimensiones de la joroba.
Teoría y ecuaciones de gobierno La ecuación de Navier-Stokes se puede dividir en tres ecuaciones para describir el
movimiento del flujo en tres sistemas de coordenadas como en las Ecuaciones (4)–(6).
Recientemente, las soluciones numéricas se han convertido en un medio para resolver
problemas complicados que son difíciles o costosos de lograr en el laboratorio. En este estudio,
Du ÿ2 u ÿ2 tu ÿ2 tu þ
se utilizó el software comercial de CFD, Flow 3D, que aplica las ecuaciones de movimiento de ÿp
ÿ ¼ ÿgx þÿ þ ÿz2 (4)
Dt ÿx ÿx2 ÿy2
fluidos para resolver las ecuaciones diferenciales de segundo orden, transitorias y no lineales
para describir el movimiento del fluido. Las ecuaciones gobernantes involucradas incluyen las
Ecuaciones de continuidad (1) y (2) y las Ecuaciones de Navier-Stokes promediadas por Dv ÿ2 v ÿ2 ÿ2 v
ÿp
ÿ ¼ ÿgy þÿ þ vþ (5)
Reynolds (3)–(6) [21]. Dt ÿy ÿx2 ÿy2 ÿz2
Dw ÿ2 w ÿ2 w ÿ2 w
ÿp
ÿÿ ÿ ¼ ÿgz þÿ + ÿy2 þ ÿz2 (6)
þ ÿ:ðÿvÞ ¼ 0 (1) Dt ÿz ÿx2
ÿt
En Flow-3D, cada celda en el dominio de fluido tiene una fracción de volumen de agua (F)
donde,
que oscila entre 0 y 1. Donde 1 representa las celdas que están completamente ocupadas con
agua, mientras que 0 representa las celdas que están completamente ocupadas con aire. Los
valores entre 1 y 0 representan la superficie libre entre el aire y el agua, como se muestra en la
Fig. 3. Flow-3D define la elevación de la superficie libre mediante el uso de la función de volumen
de fluido (VOF) en la Ecuación (7).
ÿF 1 ÿÿÿ FAxu
þ ðFAxuÞ þ R ðFAzwÞ þFAyv
ÿ þ ÿx ¼ FDIF þ FSOR (7)
ÿt FV ÿy ÿz X
Tabla 1
2
Machine Translated by Google
Fig. 4. Mallado y geometría después de aplicar el solucionador FAVOR del interior de Flow-3D.
w son componentes de velocidad en las direcciones x, y y z. aguas abajo de la joroba. Se llevaron a cabo cuatro casos de experimentos con diferentes
descargas de flujo Q y pendientes de lecho So , como se muestra en la Tabla 1.
Materiales y métodos El medidor de corriente utilizado para las mediciones de velocidad registró el número de
rotaciones (p) durante 30 s y el número de rotaciones de la hélice por segundo (n) se calculó
A. Pruebas de laboratorio utilizando n ¼ p/30. Las velocidades se calcularon usando las Ecuaciones (8) y (9) a continuación:
tanques de almacenamiento ubicados debajo del canal. Para reducir la turbulencia del flujo en
la entrada, se colocaron unos deflectores en la entrada (aguas arriba) del canal. Durante los
experimentos, se utilizó un calibrador puntual para medir la profundidad del agua y un medidor B. Configuración numérica
de flujo (hidrómetro SEBA) para medir la velocidad del flujo en diferentes lugares, aguas arriba
y En este estudio, se utilizó el código comercializado por CFD (Flow-3D) para simulaciones
numéricas. El software utilizó el método de volumen finito.
3
Machine Translated by Google
Se ubicaron 4 sondas de historia a lo largo del cauce, 2 sondas colocadas antes de la joroba y el
(FAVOR) se utilizó para definir los efectos del mallado en la geometría. Se definieron bloques de malla
resto colocadas después de la joroba. Estas sondas registraron datos hidráulicos, incluida la velocidad
hexaédrica para las simulaciones y el número total de bloques de malla fue de 1.920.000 celdas. La
del flujo, la profundidad del flujo, el número de Reynold, etc.
Fig. 4 muestra la malla y la geometría después de aplicar el solucionador FAVOR. Las condiciones de
Se simularon cuatro casos con diferentes caudales de caudal y taludes de lecho con un tiempo de
contorno se definieron con la entrada aguas arriba definida usando la velocidad del flujo y la
simulación promedio de 6 h cada uno. Se ha realizado un estudio de independencia de malla para
profundidad del flujo, mientras que aguas abajo se definió como una salida con presión 0.
encontrar el tamaño de malla adecuado y preciso para este caso.
Fig. 8. Perfil de agua numérico vs experimental caso 2. Fig. 10. Perfil de agua numérico vs experimental caso 4.
4
Machine Translated by Google
Fig. 11. a- Contorno de variación de la velocidad del caudal b- Contorno de variación del número de fraudes c- Contorno de variación de la profundidad del caudal d- Contorno de variación de la presión.
5
Machine Translated by Google
Fig. 13. Variación del número de Froud aguas arriba y aguas abajo.
Fig. 14. Variación de la velocidad media en profundidad aguas arriba y aguas abajo.
características. Se puede observar que la descarga más pequeña (30 m3 /h) tiene la menor Para acceder a la precisión del software Flow 3D, la comparación entre los resultados
profundidad de agua. numéricos y experimentales, los perfiles de la superficie del agua para todos los casos
probados se extrajeron de la simulación numérica y
A. Validación de resultados experimentales con trazado como se muestra en las Figs. 7–10. A partir de los gráficos, un fuerte acuerdo puede
6
Machine Translated by Google
se observa entre los resultados numéricos y experimentales como las variaciones oscilan otros modelos de turbulencia dentro de Flow 3D pueden considerarse una gran
entre el 4% y el 5%. Esta variación es similar a la de contribución científica al campo y, por lo tanto, son muy recomendables para futuros
el estudio de [22] quienes observaron variación entre 1 y 3.5%. estudios.
Además, la velocidad promedio de profundidad y las variables del número de Froud
se extrajeron de la simulación numérica y se compararon con los resultados de los
experimentos como se muestra en las Tablas 3 y 4. Se puede notar una buena Conflicto de intereses
concordancia entre ambos resultados en la ubicación aguas arriba (EE. UU.) y aguas abajo (DS).
Los autores declaran no tener conflicto de intereses.
B. Visualización de contornos de flujo
Reconocimiento
Para comprender la física de la velocidad, el número de Froude, la profundidad del
flujo y las variaciones de presión en el flujo con joroba, los gráficos de contorno se Los autores agradecen a Universiti Teknologi PETRONAS, Malasia por proporcionar
muestran en la Fig. 11. La Fig. 11a muestra la variación de la velocidad del flujo de aguas el entorno propicio y los equipos de laboratorio hidráulicos que se utilizarán para esta
arriba a aguas abajo. Se puede observar claramente que el flujo tiene mayor velocidad investigación y a iRMC Bold2025, número de subvención RJO10436494, Universiti
aguas abajo que aguas arriba, esto se debe al efecto de la joroba que generó más presión Tenaga Nasional, Malasia.
de cabeza en el lado de aguas arriba y menos en el lado de aguas abajo, como se
muestra en la Fig. 11d. Referencias
Además, la Fig. 11b muestra las variaciones del número de Froude desde la etapa
subcrítica hasta la supercrítica aguas abajo con el valor más alto en la parte superior del [1] RL Francia, Diseño de humedales: principios y prácticas para arquitectos paisajistas y planificadores
extremo derecho de la joroba. Esto se debe a la geometría de la joroba y la naturaleza del uso de la tierra, WW Norton and Company, 2003.
[2] A. Ferrari, simulación SPH de flujo en superficie libre sobre un vertedero de cresta afilada, Adv.
de los puntos de transición. La Fig. 11c muestra la variación de la profundidad del flujo
Recurso de agua . 33 (3) (2010) 270–276.
entre aguas arriba y aguas abajo de la joroba. Además, se puede observar una vista 3D [3] R. Singh, D. Manivannan, T. Satyanarayana, coeficiente de descarga de rectangular
más clara de la velocidad del flujo en la Fig. 12 con la velocidad más alta aguas abajo y vertederos laterales, J. Irrig. Fuga. Ing. 120 (4) (1994) 814–819.
la más baja aguas arriba debido a la mayor presión creada por la presencia de la joroba. [4] SM Borghei, MR Jalili, M. Ghodsian, Coeficiente de descarga para vertedero lateral de cresta afilada
en flujo subcrítico, J. Hydraul. Ing. 125 (10) (1999) 1051–1056.
[5] A. Al-Shukur, M. Al-jumaily, Z. Shaker, Experimental research of flow char
Además, a partir de los resultados numéricos obtenidos con Flow 3D, las variaciones características sobre el vertedero crump con diferentes condiciones, Saudi J. Eng. Tecnología 2
(10) (2017) 373–379.
aguas arriba y aguas abajo en tres parámetros, incluidos el número de Froude, la
[6] KS Erduran, G. Seckin, S. Kocaman, S. Atabay, Modelado numérico 3D del flujo alrededor del
velocidad del flujo y la profundidad del flujo, se muestran en las Figs. 13-14 para los
cruce de puentes sesgados, Eng. aplicación computar Mecánica de fluidos 6 (3) (2012)
cuatro casos de descargas. De las figuras, generalmente se puede observar que el flujo 475–489.
aguas arriba es de naturaleza subcrítica con valores del número de Froude menores a [7] ASI Abdurrasheed, KW Yusof, HB Takaijudin, A. Ab, BS Iskandar, Efectos del remanso en la
1.0 mientras que aguas abajo los flujos son supercríticos con valores del número de evaluación del rendimiento hidráulico de los módulos rainsmart en sistemas de drenaje
sostenibles, en: Conferencia Internacional sobre Recursos Hídricos, Langkawi, Malasia, noviembre
Froude mayores a 1.0 como en la Fig. 13. Además, se encontró que la velocidad del flujo de 2018 , págs. 27 y 28.
generalmente aumenta a medida que el flujo se mueve de aguas arriba a aguas abajo en [8] PG Chanel, JC Doering, Evaluación del modelado de aliviaderos usando dinámica de fluidos
computacional , Can. J. Civ. Ing. 35 (12) (2008) 1481–1485.
los cuatro casos, como en la Fig. 14. Por otro lado, se encontró que la profundidad del
[9] S. Dehdar-Behbahani, A. Parsaie, Numerical modeling of ow pattern in dam
flujo disminuye de aguas arriba a aguas abajo en todos los casos . evidente en la Fig. 15. Muro guía de aliviaderos. Estudio de caso: represa balaroud, Irán, Alejandría Ing. J. 55 (1)
(2016) 467–473.
[10] A. Parsaie, AH Haghiabi, A. Moradinejad, modelado CFD del patrón de flujo en aliviaderos
canal de aproximación, Sustain. Recurso de agua. Administrar 1 (3) (2015) 245–251.
[11] HK Versteeg, W. Malalasekera, Una introducción a la dinámica de UID computacional:
Conclusiones el Método del Volumen Infinito, Pearson Education, 2007.
[12] ASI Abdurrasheed, K.,W. Yusof, HB Takaijudin, AA Ghani, MM Muhammad, AT Sholagberu,
Avances y desafíos en la tecnología de módulos de drenaje subterráneo y BIOECODS: una
Este estudio modeló las características del flujo sobre una joroba triangular de cresta revisión, en: MATEC Web of Conferences, 203, EDP Sciences, 2018, 07005.
afilada usando simulación numérica y validado con resultados experimentales. A partir de
los resultados y las discusiones, se extraen las siguientes conclusiones: El software Flow [13] SY Kumcu, Investigación del modelado de flujo sobre aliviadero y comparación entre datos
experimentales y análisis CFD, KSCE J. Civil Eng. 21 (3) (2017) 994–1003.
3D modeló las características del flujo a través de una joroba triangular con un mayor
[14] M. Darw, F. Moukalled, M. Luca, Método de volumen finito en fluido computacional
grado de rendimiento y precisión. Esto se debe a que se observó una fuerte concordancia Dinámica: una introducción avanzada con OpenFOAM y Matlab, Springer, 2015.
entre los resultados numéricos y experimentales con variaciones mínimas que oscilan [15] H. Zahabi, M. Torabi, E. Alamatian, M. Bahiraei, M. Goodarzi, Efectos de la geometría y las
características hidráulicas de los embalses poco profundos sobre el atrapamiento de sedimentos,
entre el 4% y el 5%. En términos de propiedades de flujo, se puede observar que los
Water 10 (12) (2018) 17–25.
flujos eran de naturaleza subcrítica y supercrítica aguas arriba y aguas abajo, [16] G. Li, X. Li, J. Ning, Y. Deng, Simulación numérica y aplicación de ingeniería de un balde en forma
respectivamente. Esto se atribuye a la presencia de la joroba a lo largo de la trayectoria de cola de milano, Agua 11 (2) (2019) 242.
[17] DC Lo, JS Liou, S. Chang, Actuaciones hidrodinámicas de los flujos de aire y agua en
del flujo y puede ser muy aplicable en el diseño de estructuras hidráulicas, especialmente
sumideros con y sin paletas generadoras de turbulencia para sistemas de drenaje de edificios,
en la estabilización de canales. Además, se encontró que las velocidades de flujo Water 7 (2) (2015) 679–696.
generalmente aumentan a medida que el flujo se mueve de aguas arriba a aguas abajo [18] T. Cebeci, Modelos de turbulencia y su aplicación: métodos numéricos eficientes con programas
informáticos, Springer Science and Business Media, 2003.
en los cuatro casos (como se ve en las vistas 2D y 3D), mientras que se encontró que las
[19] B. Mohammadi, O. Pironneau, Análisis del modelo de turbulencia k-epsilon, 1993.
profundidades de flujo disminuyen de aguas arriba a aguas abajo en todos los casos. En [20] S. Patankar, Numerical Heat Transfer and Uid Ow, CRC press, 1980.
general, Flow 3D se puede recomendar como una herramienta útil para predecir el [21] L. Choufu, S. Abbasi, H. Pourshahbaz, P. Taghvaei, S. Tfwala, Investigación del patrón de flujo,
comportamiento de los vertederos en las mediciones de descarga, la estabilización del erosión y sedimentación alrededor de espigones variados bajo diferentes condiciones hidráulicas
y geométricas: un estudio numérico, Agua 11 (2 ) (2019) 235.
flujo aguas arriba y aguas abajo. Además, realizar un modelado similar con [22] S. Haun, NRB Olsen, R. Feurich, Numerical modeling of ow over trapezoidal wide -crested weir,
Eng. aplicación computar Mecánica de fluidos 5 (3) (2011) 397–405.