Traducido del inglés al español - www.onlinedoctranslator.

com

www.nature.com/scientificreports

ABIERTO Evaluación experimental

de la eficiencia de disipación de energía de la
sección de flujo de vórtice de los ejes de caída

Mohammad Mahmoudi‑Rad1*&Mohammad Najafzadeh2

En las redes de recogida y drenaje de aguas residuales urbanas, se utilizan estructuras de vórtice para transferir
fluido entre dos conductos con diferencias de nivel significativas. Durante el eje de caída, además de evitar que el
fluido caiga debido a la formación de un flujo de vórtice, una cantidad significativa de la energía del fluido se
disipa debido a la fricción de la pared del eje vertical. En el presente estudio, mediante la construcción de un
modelo físico a escala 1:10 fabricado en plexiglás, se ha investigado la eficiencia de disipación de energía en el eje
vertical. De esta manera, la realización del análisis dimensional indica que el número de Froude de flujo (fr.) y la
relación entre la altura total de caída y el diámetro del eje (L⁄D) son parámetros que afectan la eficiencia de la
disipación de energía del flujo en el eje vertical (ηs). Esta investigación considera cuatro niveles defr.factor (1,77,
2,01, 2,18 y 2,32) y tres niveles deL⁄Dfactores (10, 13 y 16). Además, cuatro repeticiones para 12 combinaciones
posibles nos permiten realizar 48 experimentos y el método factorial completo. Los resultados demostraron que
la eficiencia de disipación de energía en los cambios de eje vertical varía de 10,80 a 62,29%. Además,ηsLos valores
disminuyen con el aumento defr.mientras que la eficiencia aumenta al aumentarL⁄Drelación. Además, el análisis
de regresión arrojó una ecuación polinómica de segundo orden que es función defr.yl⁄Dpara estimar con
precisión la eficiencia de disipación de energía del flujo en el eje vertical.

En los sistemas de alcantarillado y drenaje urbanos, la transferencia de fluidos desde el nivel superior al nivel inferior se realiza
ocasionalmente a través de estructuras de descenso.1. Los dos tipos más comunes de estas estructuras son los pozos de acceso y
los pozos de descenso.2,3. Las arquetas de bajada se utilizan para alcanzar el objetivo de desperdiciar energía y reducir la velocidad
del flujo en diferencias de nivel inferiores a 7 m.4mientras que los pozos de caída se utilizan en desniveles de más de 7 a 10 m5. Los
ejes de caída se dividen en dos grupos según el tipo de flujo formado: ejes de caída de vórtice y ejes de caída de inmersión.6,7.
Mantener un patrón de flujo constante y desperdiciar más energía en diferentes descargas hace que los ejes de caída de vórtice
sean superiores a los ejes de caída con ejes de caída hundidos.8,9. La estructura del vórtice consta de tres partes principales: la
estructura de entrada, el eje vertical y la estructura del disipador de energía. La cierta geometría de la estructura de entrada que
se forma en varias configuraciones, como la entrada en espiral.10,11, entrada tangencial8y desplácese hacia abajo12, es capaz de
convertir el flujo de aproximación en flujo de vórtice13. El flujo de vórtice se adhiere a la pared vertical del pozo para formar un
núcleo de aire en el centro de la columna y se mueve hacia abajo de forma anular. Además, la fricción de la pared es el factor más
importante en la disipación de energía. Esta disipación de energía en el eje vertical depende principalmente de la fricción de la
pared y de la turbulencia del flujo de vórtice. Por el contrario, las turbulencias y los procesos de arrastre de aire en el fondo del
pozo de caída no influyen en la disipación de energía en el pozo de caída.5. En la estructura de salida, además de la desaireación, la
energía restante del flujo se disipa en gran medida.14.
Las estructuras de caída de vórtice se emplean básicamente en diferentes campos. De esta manera, se introducen los ejemplos
prácticos más importantes como (i) aguas residuales municipales del área de Milwaukee con descarga de hasta 90 m3/s en una
altura de caída de 80 m15, (ii) la descarga de 140 m3/s y la altura de 170 m, Curban, Italia14, (III) central eléctrica de Shapai en China
con un caudal de 200 m3/s y la altura de caída de unos 100 mdieciséis, y (IV) transmisión de descarga de 1400 m3/s con una altura de
caída de 190 m hasta el túnel de desvío en las centrales eléctricas de Xiaowan en China17.
Generalmente, se puede esperar una fracción considerable de la disipación de energía del flujo en los ejes de caída del vórtice. Los
resultados de la investigación experimental descubrieron que la disipación de energía del flujo probablemente se debe a la relación entre la
altura y el diámetro (L/D) del eje. En los estudios anteriores, las fracciones de 85% y 90% de disipación de energía del flujo han

1Departamento de Ingeniería Civil, Complejo de Educación Superior de Bam, PO Box 76615314, Bam, Irán.2Departamento de
Ingeniería Hidráulica, Facultad de Ingeniería Civil y Topografía, Universidad de Posgrado de Tecnología Avanzada,
PO Box 76315116, Kerman, Irán.*correo electrónico: m.m6456@yahoo.com ; m.mahmoudirad@bam.ac.ir

Informes Científicos| (2023) 13:1679 |https://doi.org/10.1038/s41598-023-28762-2 1

Vol.:(0123456789)
 www.nature.com/scientificreports/

estimado para el eje de caída conl=50D14yl=100D15, respectivamente. Estos hallazgos se han establecido basándose en el
supuesto de que los pozos de caída con una gran relación altura-diámetro pueden permitir que el flujo alcance la
velocidad terminal. Por el contrario, para ejes de caída con una relación altura-diámetro relativamente pequeña, se han
reportado proporciones de disipación de energía del 62% y 34%.l=9D18y por la mitad de la longitud de un eje de caída con
l=14D19, respectivamente. Crispino et al.20indicó que la eficiencia de la pérdida de energía dependía principalmente de los
efectos del flujo y la turbulencia que tenía lugar en la sección de la cámara. Propusieron una relación derivada de un
modelo teórico sencillo para aproximar el coeficiente de pérdida de energía. En el caso de los sistemas de eje de caída-
túnel, el estudio experimental realizado por Chan y Chiu21indicó que el número de Froude del eje de caída y la relación
entre el diámetro del eje de caída y el diámetro del túnel tuvieron profundos impactos en la formación de diversas
estructuras de flujo.
El diseño de experimentos (DoE) implica métodos que examinan los efectos específicos sobre los resultados (respuestas)
mediante la realización de cambios específicos en uno o más factores. En estos métodos, para lograr la respuesta, se prueban los
factores simultáneamente y también se consideran las interacciones entre ellos. Refiriéndose a22,23, considerando el diseño de
experimentos (es decir, factorial completo), el propósito es mostrar primero la posibilidad de investigar los principales factores
efectivos sobre la respuesta del modelo físico. Para mayor explicación, se puede decir que después de seleccionar los distintos
niveles para cada factor considerando las repeticiones (al menos dos repeticiones para calcular el error) para todas las
combinaciones posibles, los resultados obtenidos se analizan utilizando el ANOVA y elpag-Índice estadístico de valor. En los
últimos años, se han realizado esfuerzos formidables para investigar la disipación de energía del flujo y el mecanismo de arrastre
de aire en las estructuras de los vórtices.22,23. A partir de estas investigaciones, se puede demostrar que la usabilidad exitosa del
DOE en el análisis de la investigación experimental demostró valores óptimos de los parámetros adimensionales de diseño al
construir estructuras de vórtice en las aplicaciones prácticas. Además, se estudió la investigación más reciente en la que se
estudiaron los efectos de la cámara de disipación sobre las pérdidas de energía en la energía del vórtice.24. Mahmoudi-Rad y
Najafzadeh24encontró que los valores óptimos de los parámetros efectivos (es decir,Fr yl/D) arrojó 2,32 y 13,901, respectivamente;
para mantener la eficiencia de la pérdida de energía del flujo en su nivel más alto. Obviamente, una gran cantidad de energía de
flujo en las estructuras de vórtice se disipa en la sección del eje vertical. Las pérdidas de energía en el eje vertical pueden deberse a
las propiedades geométricas de los ejes y a la velocidad del flujo que se aproxima aguas arriba a la sección de entrada de las
estructuras de vórtice. Bajo el auspicio de los estudios más relevantes, sólo un factor efectivo (i.mi., número de Froude que se
refiere a la condición hidráulica) se ha investigado en varios niveles14,15,18,19mientrasl/DLa proporción se mantuvo constante
durante las investigaciones. Las principales diferencias entre el presente estudio y otras investigaciones previas sobre la pérdida
de energía en los fustes verticales es que se consideraría el efecto simultáneo de dos factores, introducidos como condiciones
hidráulicas y geométricas, en muchos niveles sobre el comportamiento de la estructura. De esta manera, la ecuación de regresión
dada para calcular la pérdida de energía del flujo dependía del número de Froude, mientras que los parámetros geométricos de
las estructuras han jugado un papel destacado en la estimación de las pérdidas. Sin embargo, existe una ferviente necesidad de
proporcionar un modelo matemático para predecir la respuesta de la estructura del eje a las variaciones en las condiciones
hidráulicas y geométricas. En el presente estudio, el mayor inconveniente de las investigaciones anteriores se ha resuelto
considerando el diseño de ensayos factoriales completos para el eje vertical.
La organización de la investigación de este estudio se presenta de la siguiente manera, (i) se introducen factores efectivos para
estimar la eficiencia de disipación de energía del flujo en el eje vertical de la estructura de vórtice mediante análisis dimensional, (ii
) los efectos simultáneos de estos factores sobre laηsen el eje vertical se investigan debido al uso limitado del diseño experimental
en la literatura anterior, (III) los experimentos se diseñan y analizan mediante el método factorial completo para evaluar la
eficiencia de disipación de energía del flujo en el eje vertical de la estructura de vórtice, (IV) se proporciona una ecuación de
regresión robusta para describir cuán efectivos son los factores principales, y (v) los resultados del presente estudio se comparan
con la literatura relevante en términos de cantidad y calidad.

Eficiencia de disipación de energía de flujo en el eje vertical.

Investigaciones anteriores demostraron queηsEl valor en el eje vertical se calcula como19,25:

s= (1-hyo+1/Hi)×100% (1)
dónde,hiyhi+1son cabeza de energía total enz/D=10.3 (sección "Introducción") yz/D=4.4 (sección "Resultados y discusiones
") en la estructura de caída del vórtice (ver Fig.1), respectivamente. Según investigaciones experimentales anteriores19,24,
la carga de energía total en cada sección del eje de caída se calcula utilizando las siguientes ecuaciones:

V2 V2t r
PAG( )
=z +z + + (2)
2gramo 2gramo ρg

∫r [ ]
V2 1 1 1
(r) = ρtdr = ρC2 − (3)
r 2 R2(1 -t)2 r2
R-b

( ) 1 ( )
q 3 4 D mi
=Vtr= gramo porque3b − (4)
mi 2 2

q 4q
z= = (5)
A 2
πDt(2-t)

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 2

Vol:.(1234567890)
www.nature.com/scientificreports/

Figura 1.Definición de eje de caída de vórtice: (a) vista lateral, (b) vista en planta de entrada tangencial, (C) fotografía del
modelo físico.

4C
t= (6)
(2 -t)D
en el cual,z=altura desde la elevación cero (lecho de origen del canal de salida),Vz=velocidad vertical,Vt=velocidad
tangencial,PAG(r) = distribución de presión en cada sección del eje de caída,r=coordenada radial,R=radio del eje vertical,D
=diámetro del eje vertical,b=espesor del flujo de vórtice,t=espesor relativo del flujo (t=b/R), C=circulación constante,q=
descarga de diseño,mi=ancho de entrada en la unión del eje vertical,gramo=gravedad de la aceleración, yb= ángulo de la
pendiente inferior.
La medición de la velocidad dentro del eje se realiza reemplazando el espesor relativo del flujo de vórtice (t) en las
ecuaciones. (5) y (6), se calculan las componentes vertical y tangencial de la velocidad. Por lo tanto, en cada sección
transversal del eje vertical, los valores promedio del espesor relativo del flujo de vórtice se reemplazan en las Ecs. (5) y (6)
el instrumento mostrado en la Fig.4. Según Zhao et al.17investigación, una combinación de cabezal de presión (referido a
la presión de la pared donder=R=D/2) y altura de velocidad tangencial obtenida:

V2t P(r) 2C2

+ = (7)
2gramo ρg gD2(1 -t)2

En este estudio, se utilizó la mano derecha de la ecuación reciente en lugar de la carga de presión más la carga de velocidad.

Análisis dimensional
Durante la última década, los resultados de investigaciones experimentales han descubierto que las variables independientes que afectan ηs
Los valores en el eje vertical se introducen como19,20,25–27:

s= ϕ(Q, B, l, e, β, L, D, f , g, ρ,µ, σ) (8)

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 3

Vol.:(0123456789)
 www.nature.com/scientificreports/

Figura 2.Piezas prefabricadas de vidrio para cambio de altura de hueco.

donde = símbolo funcional, = densidad del fluido, = viscosidad dinámica, = tensión superficial. Según la Fig.1, otras
variables sonB=ancho del canal de acceso,yo=longitud de entrada tangencial,l=caída de altura total yF=coeficiente de
fricción. Cabe señalar que la profundidad del flujo (h) se ha medido mediante piezómetros instalados en el fondo del
canal de entrada y, además, la velocidad media del flujo (V) al comienzo de la entrada tangencial se calcula mediante
variables conocidas deqyB. Considerandoq,Dy teorema como de Buckingham:
parámetros repetidos y la ecuación. (9) se reorganiza usando

( )
QV2hρV B le L
s= ϕ β, f , , ,√,,,, (9)
νh σ gh DDDD
√
dónde,r= Q/h,r= V/gh,y =V2hρ/σson el número radial de Reynolds, el número de Froude y el número de Weber
número de flujo de aproximación respectivamente. En el presente estudio, el número de Reynolds radial varía de 213.000 a
360.000 y adicionalmente el número de Weber está entre 861 y 4066. Según Mulligan et al.28investigación, los efectos de la tensión
superficial y la viscosidad para el flujo de vórtice se ignoraron debido al hecho de queW.yRreran mayores que 120 y 1000. Por lo
tanto, la ecuación. (10) se reescribe de la siguiente manera:
( )
B le L
s= ϕ β, f , Fr, , , , (10)
DDDDD

A través de este estudio,b,F,B/D,l/D, ymi/Dtienen valores constantes que son 29,7°, 0,02, 1,125, 2,18 y 0,25, respectivamente.
Así, la ecuación. (10) se puede reescribir como,
( )
l
s= ϕ Fr , (11)
D

Diseño del modelo físico

El modelo físico de la estructura del vórtice de la red municipal de alcantarillado de Teherán con una escala de 1:10 está hecho de plexiglás
transparente. Cifra1c muestra una representación esquemática de la configuración experimental. El modelo físico consta de un canal de
aproximación rectangular, una entrada tangencial, un eje vertical, una cámara de disipación y un túnel de salida rectangular. La estabilidad
del flujo de vórtice, la circulación adecuada del aire y la prevención de fluctuaciones del flujo dependen directamente de la determinación del
diámetro apropiado para el eje de caída. Por lo tanto, jainista8propuso la relaciónD=k× [q2/gramo]0,2 d
para determinar el diámetro del eje de caída. En una relación reciente,kes el factor de seguridad ydes la descarga máxima de
diseño. El rangok=1–1.25 conduce al diseño económico de la estructura de vórtice.8,10. Para el presente modelo de investigación
con un caudal máximo de diseño igual a 19,4 l/s, el valor deD=0,16 m dieron como resultadok=1.22. Esta investigación utilizal/D
ratios de 10, 13 y 16, y adicionalmente dos piezas prefabricadas conl/D=3, como se muestra en la Fig.2. Además, para ajustar la
descarga del flujo y, en consecuencia, el número de Froude del flujo de aproximación (fr.) en cuatro niveles (1,77, 2,01, 2,18 y 2,32),
se utilizó un caudalímetro electromagnético (MFC 300, Irán Farasanj Abzar, Teherán, Irán) con una precisión del 1%. En esta
investigación, el núcleo de aire formado en la entrada del eje vertical se evalúa mediante una regla de cuatro patas de diseño
propio, como se muestra en la Fig.3. Utilizando varillas con un diámetro de 2 mm y

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 4

Vol:.(1234567890)
www.nature.com/scientificreports/

Figura 3.Dispositivo fabricado para medir el espesor del flujo del vórtice.

creando intervalos de 5 mm con diferentes colores en cada uno de ellos, el error máximo en la medición del valor
dimétrico del núcleo de aire se limita a ± 1 mm.

Método factorial completo

El diseño de experimentos (DoE) se utiliza generalmente para investigar el efecto de los principales factores y sus interacciones sobre el
fenómeno estudiado.21,23,29. Independientemente de cualquiera de los factores en el diseño de las pruebas, es imposible mencionar con
firmeza cuán efectivos son estos parámetros adimensionales. Según la literatura, las variaciones defr. yl/DSe han investigado los factores que
influyen en la pérdida de energía del flujo en el eje vertical.14,15,18,19. Esta investigación tuvo como objetivo estudiar los efectos de los factores
principales sobre la respuesta mediante el uso de un diseño de prueba factorial completo. Estos dos factores se encuentran entre los factores
que se han evaluado en las investigaciones más recientes relacionadas únicamente con la estructura del vórtice. Por otro lado, dado que el
sellado de la estructura del vórtice es el problema más desafiante al cambiar las condiciones geométricas de un modelo de laboratorio, la
investigación de varios factores geométricos al mismo tiempo será mucho más difícil. Además, con la adición de una o dos pruebas
factoriales más (teniendo en cuenta al menos tres niveles para estos factores), el número de pruebas aumentó de 48 experimentos a 144 y
432, respectivamente. No se puede negar el hecho de que un aumento en el número de pruebas afecta de manera eficiente los costos de
realización de la investigación.
En este estudio, se ha utilizado un diseño factorial completo para investigar el efecto defr.yl/Dfactores enηsvalores.
Considerando 4 y 3 niveles parafr.yl/D, un diseño factorial completo incluye 12 (4 × 3) combinaciones posibles de estos
factores. Aplicando 4 repeticiones para cada combinación posible se realizan 48 (12 × 4) experimentos y adicionalmente,
Tabla1Da detalles del diseño de los experimentos.

Resultados y discusiones
Estudio paramétrico.Cifra3proporciona a los lectores información sobre la variación significativa de la eficiencia de
disipación de energía del flujo en el eje vertical (ηs) versus el número de Froude (fr.) para varios niveles del/D. A simple
vista, queda muy claro que los cambios deηscontrafr.valores en todos los valores del/Dhan atravesado tendencias a la
baja. Para todos los valores defr.,ηsaumentó con el aumento de los valores del/D. Por ejemplo, parafr.=1,77,ηsLos valores
aumentaron de aproximadamente el 32% enl/D= 10 a aproximadamente 62% enl/D= 16. Además, cuando se trata de un
cierto valor del/D,ηsLos valores se desplomaron con un aumento enfr.número. Como ejemplo, la Fig.3demostró que, porl
/D= 13,ηslos valores disminuyeron de alrededor del 51% enfr.=1,77 a aproximadamente el 25% enfr.=2.32. Esta
investigación sigue el 2Dakota del Norteregresión polinómica de orden para la descripción de variaciones deηsversusl/Den
varios niveles de fr.números. Como ejemplo, paral/D= 10, los resultados del análisis estadístico indicaron que la
expresión de segundo orden proporcionó la predicción más precisa deηsvalores (coeficiente de determinación [R2] =
0,9975) en Material complementario (ver Fig. S1) en comparación con otros tipos de ecuaciones de regresión: lineal (R2=
0,932), potencia (R2= 0,8122), logarítmico (R2= 0,9112), y exponencial (R2= 0,8407). Todos los rendimientos de las
ecuaciones de regresión para varios niveles del/Dfueron presentados en la tabla2.

Patrón de flujo en el eje vertical.La correcta formación del flujo de vórtice a lo largo del eje vertical está influenciada por el
núcleo de aire que se forma al principio. Según Yu y Lee27estudio, la relación del área del núcleo de aire (λ) se aproxima porλ=d2/D
2en el cualdes igual al diámetro del núcleo de aire en la entrada del eje vertical. Como se ve en la Fig.4, el flujo de vórtice en el eje
vertical permanecerá estable si la relación λ es mayor que 0,25. El núcleo de aire medido en el

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 5

Vol.:(0123456789)
 www.nature.com/scientificreports/

fr.
L/D 1,77 2.01 2.18 2.32
33,89 28.21 21.38 11.06

34.04 28,73 21.51 10.80

10
33,79 28.60 21.46 11.06

34.04 28.60 21.14 12.01

50.34 44,73 37,45 22,84

50,37 44,74 37.24 23.45

13
50,48 45.16 37,48 22,98

50,51 45.23 37,47 23.30

62.10 57,84 51.02 36,99

62,28 58.30 51,69 37,53

dieciséis
62,29 58.20 51,78 37,60

62.24 58.17 51,79 37,81

Tabla 1.Diseño factorial completo conηs(%) valores relacionados con cada prueba.

l/D Tipo de regresión ηs R2

Polinomio − 65,14fr.2+ 225,66 francos–161,5 0.9975

Exponencial 1087.5e−1,885fr. 0.8407

10 Lineal − 40.042fr.+106,66 0.932

logarítmico − 80,23Ln(fr.) + 81.737 0.9112

Fuerza 331.02fr.−3,755 0.8122

Polinomio − 98.856fr.2+ 356,12 francos–270,44 0.9924

Exponencial 565.86e-1.312fr. 0.8183

13 Lineal − 47.115fr.+136,51 0.8885

logarítmico − 94,12Ln(fr.) + 106,99 0.8634

Fuerza 246fr.−2.61 0.7887

Polinomio − 103,66fr.2+ 380,55fr.–286,85 0.9885

Exponencial 300e-0,852fr. 0.8059

dieciséis Lineal − 42,28fr.+139,88 0.8524

logarítmico − 84,28Ln(fr.) + 113,24 0.8246

Fuerza 174,83fr.−1.694 0,7755

Tabla 2.Resultados de varios ajustes de curvas para comprender las variaciones.ηsde contraP.

comienzo del eje vertical para diferentes valores de lafr.El número se muestra en la Fig.4. Se puede observar que la relación del
área del núcleo de aire disminuye con un aumento en lafr.valores numéricos (o descarga de flujo). Paraq= 1,7qd= 33 l/s, la
perturbación de la circulación del aire hace que la superficie libre de flujo se vuelva horizontal en la entrada de la estructura y el
flujo de vórtice desaparezca (Fig.5). Cifra6a – c muestra los cambios en el espesor relativo del flujo (t), velocidad vertical (Vz), y
velocidad tangencial (Vt) a lo largo del eje vertical para descarga de diseño y diferentesl/Dniveles, respectivamente. Cifra6a ilustra
que, para diferentes niveles del/Drelación, los valores medios del espesor relativo del flujo en niveles más altos (z/D> 5,9)
disminuyen rápidamente y se observa un ligero aumento en niveles más bajos. Por ejemplo, enl/D= 13, los valores medios de
espesor relativo disminuyen de 0,24 pulgadasz/D= 10,3 a 0,16 pulgadasz/D= 5,9 mientras que paraz/D= 4,4, el valor medio del
espesor relativo aumentó ligeramente a 0,17. Cifra6b ilustró que los valores de velocidad vertical tenían tendencias ascendentes
en elz/D> 5,9, mientras que se observaron tendencias opuestas en la elevación más baja (z/D= 4,4–5,9). Por ejemplo, paral/D= 16,
los valores de velocidad vertical aumentan de 2,08 m/s enz/D= 13,3 a 3,21 m/s enz/D= 5,9 y luego disminuyó a 2,86 m/s enz/D= 4,4.
La reducción de los valores de la velocidad vertical en los niveles inferiores del eje vertical provoca un impacto más débil del flujo
de gotas dentro de la cámara de disipación. Este débil golpe hace que el flujo de energía residual en la cámara de disipación se
desperdicie. Cifra6c demuestra que, para diferentes niveles del/Dfactor, la velocidad tangencial en el nivel superior (z/D> 5,9)
tuvieron tendencias a la baja en elz/D> 5,9 y mientras que en niveles inferiores se produce un ligero aumento. Por ejemplo, enl/D=
10, los valores de velocidad tangencial disminuyen de 1,18 m/s en z/D= 7,4 a 1,15 m/s enz/D= 5,9 y luego aumenta a
aproximadamente 1,16 m/s enz/D= 4,4. Dado que la velocidad tangencial para diferentes valores de lal/DEl factor varía entre 1,13
m/s y 1,20 m/s, por lo que a partir de la proximidad de los valores obtenidos, se puede concluir el patrón de flujo del vórtice de
estabilidad en el eje vertical. Según el patrón de flujo presentado en la Fig.6a-c, el nivel dez/D= 5,9 es la sección transversal del
cambio de tendencia del flujo del vórtice en el eje vertical. Mesa3presenta los valores de eficiencia de disipación de energía del
flujo en diferentes partes del eje vertical (distancias entre las secciones 1-2, 2-3, 3-4 y 4-5 en la Fig.1). Como se puede observar,
parafr.≥ 2,01, la eficiencia de pérdida de energía entre secciones "Eficiencia de disipación de energía de flujo en el eje vertical." & "
Diseño del modelo físico" es mínimo y para fr.=1.77 la eficiencia mínima de pérdida de energía se produce entre tramos "Diseño
del modelo físico" y "factorial completo

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 6

Vol:.(1234567890)
www.nature.com/scientificreports/

fr.=1,77 fr.=2.01
λ = 0,61
λ = 0,48

fr.=2.18 fr.=2.32
λ = 0,38 λ = 0,31
d

Figura 4.Medición del núcleo de aire al inicio del eje vertical.

Figura 5.Perturbación observada en la formación del flujo de vórtice paraq> 33 l/s.

método". Por otro lado, para todos los valores defr.número (fr.=1,77–2,32), la mayor eficiencia de disipación de energía
del flujo se produce en las partes inferiores del eje vertical (intervalo entre secciones "Método factorial completo" & "
Resultados y discusiones"). Según la tabla3, las eficiencias máxima y mínima de disipación de energía del flujo entre
diferentes niveles del eje vertical son 26,43% (fr.=1,77) y 1,88% (fr.=2.32), respectivamente.

Modelo matemático propuesto deηs.El análisis de varianza es uno de los métodos estadísticos utilizados para
analizar datos de laboratorio y expresar la relación entre ellos. Con este método, el modelo matemático que rige los
resultados se puede extraer generando una regresión de datos y examinando los errores. También el

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 7

Vol.:(0123456789)
 www.nature.com/scientificreports/

14 14
L/D = 16

12 L/D = 13 12
L/D = 10

10 10

z/D 8

z/D
8

6 6
L/D = 16
4 4
(a) L/D = 13
L/D = 10
(b)
2 2
0.1 0,15 0,2 0,25 0.3 1.5 2 2.5 3 3.5
t=b/R Vz(EM)

14
L/D = 16

12 L/D = 13
L/D = 10

10
z/D

4
(C)
2
1,12 1,14 1,16 1,18 1,2 1,22
Vt(EM)

Figura 6.Cambios en los componentes del flujo de vórtice a través del eje vertical para varios valores del/D: (a) espesor relativo
del flujo de vórtice, (b) velocidad del flujo vertical, y (C) velocidad del flujo tangencial.

1−2(%) 2−3(%) 3−4(%)

fr. s s s s4−5(%)
1,77 12.21 13.75 11.00 26.43

2.01 10.78 10.16 12.45 21.56

2.18 8.45 7.03 9.62 18.63

2.32 4,99 1,88 5.62 12.64

Tabla 3.Eficiencia de disipación de energía entre ambas secciones transversales del eje vertical.

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 8

Vol:.(1234567890)
www.nature.com/scientificreports/

Fuente Estimación del coeficiente (i) Suma de cuadrados Grado de libertad Cuadrado medio valor F pag-valor Significativo

Modelo – 10.730,91 5 2146.18 1418.16 < 0,0001 *

Interceptar − 328.3180 – 1 – – –
X1:fr. 325.6207 3912.81 1 3912.81 2585.51 < 0,0001 *
X2:L/D 8.2099 6465.25 1 6465.25 4272.11 < 0,0001 *
X1X2 − 0,3739 1,69 1 1,69 1.12 0.2963

X21 − 89.2151 426,65 1 426,65 281,92 < 0,0001 *

X2 − 0,1028 9.13 1 9.13 6.03 0.0183 *
Residual – 63,56 42 1,51 – –

Tabla 4.ANOVA del modelo cuadrático completo [Ec. (12)].

Fuente Estimación del coeficiente (i) Suma de cuadrados Grado de libertad Cuadrado medio valor F pag-valor Significativo

Modelo – 10.729,22 4 2682.31 1767.54 < 0,0001 *

Interceptar − 318.2574 – 1 – – –
X1:fr. 320.7605 3912.81 1 3912.81 2578.40 < 0,0001 *
X2:L/D 7.4360 6535.67 1 6535.67 4306.78 < 0,0001 *
X21 − 89.2151 426,65 1 426,65 281.15 < 0,0001 *
X2 − 0,1028 9.13 1 9.13 6.01 0.0183 *
Residual – 65,25 43 1,52 – –
Modelo de indicadores estadísticos.

2 0.994 COV (%) 3.21

2 0.993 AP 132.936
adj.
2 0,992
pred

Tabla 5.ANOVA de la ecuación. (13) con sus indicadores estadísticos.

Se medirá la interacción de los factores de laboratorio en este método.29. Según las variaciones en los valores de respuesta en la
Fig. S1, el modelo matemático inicial [Ec. (12)] para la eficiencia de disipación de energía del flujo en el eje vertical (ηs) fue
seleccionado como polinomio de segundo orden:

s=c0+C1(Fr) +C2(L/D) + c3(Fr)(L/D) + c4(Fr)2+C5(A/D)2 (12)

donde C0a c5Los coeficientes de nota se obtienen utilizando el método de mínimos cuadrados. La adecuación del modelo
cuadrático ha sido evaluada porR2y los demás parámetros estadísticos se adquirieron del análisis de varianza (ANOVA). Se obtuvo
la significancia estadística del modelo y sus términos.pag-valor al 5% de nivel de significancia. Valores de los coeficientes de
regresión (C0aC5) y resultados del ANOVA para el modelo inicial [Ec. (13)] se muestran en la tabla4. En consecuencia, ANOVA
demuestra que, para elηsparámetro, los efectos lineales y cuadráticos de los factores son muy significativos (pag-valor < 0,05). El
grandepag-valores para el efecto de interacción de los factores (pag-valor > 0,05) presentado en la tabla4demostrar que los
efectos defr.yl/DLos factores son solo factores que afectan la pérdida de energía en el eje vertical. Además, dado que sólo los
términos conpag-el valor < 0,05 permanece en el modelo, el efecto de interacción defr.yl/D Los factores se eliminan del modelo
inicial y luegoC3= 0 se considera. Después de eliminar el efecto de interacción de fr.yl/Dfactores del modelo inicial, el modelo
cuadrático reducido final se obtuvo volviendo a analizar el ANOVA de la siguiente manera:

s= -318.2574 + 320.7605(Fr) +7,4360(L/D) −89.2151(Fr)2− 0,1028(L/D)2 (13)

Los resultados del análisis de varianza del modelo cuadrático reducido con algunos indicadores estadísticos se
presentan en la Tabla5. Los siguientes resultados se pueden obtener del análisis ANOVA:

i. El modelo final conpag-valores < 0,0001 indica que el modelo es estadísticamente significativo con un nivel de confianza del
99,99%.
ii. Según la tabla5,R2= 0,994 indica que el 99,4% de los cambios de respuesta se pueden describir.
2 El valor siempre es igual o menor queR2. De esta manera, se puede decir que el
III. En los modelos reducidos, el adj.
El modelo reducido tiene un número razonable de términos matemáticos.30. De la mesa5, la calidad del modelo reducido
se muestra en un número razonable de términos matemáticos con
2
adj.=0.993y r2= 0,994.
IV. El índice estadístico de2predmuestra la capacidad del modelo para predecir un conjunto de datos nuevos. Según la tabla5,
los altos valores de este índice ( 2
pred=0,992)mostrar la capacidad del modelo final [Ec. (13)] a este respecto.
v. El coeficiente de varianza (COV) es la relación entre la desviación estándar y el valor medio de la respuesta
medida (como porcentaje). El COV es el factor de reproducibilidad de un modelo matemático. En el otro

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 9

Vol.:(0123456789)
 www.nature.com/scientificreports/

70
Línea de mejor ajuste

60 Ec. (13) [R2=0,994]

+ 5% de línea de error
50
- Línea de error del 5%

Predichoηs(%)
40

30

20

10

0
0 10 20 30 40 50 60 70
Medidoηs(%)

Figura 7. El medidoηsversus los predichos por la ecuación. (13).

Por otro lado, dado que los modelos con COV<10% son modelos reproducibles, el modelo final con COV=3,21% también se
beneficia de esta capacidad.
vi. La precisión adecuada (AP) representa la relación entre la diferencia en el valor de respuesta predicho del
modelo con el valor promedio del error predictivo. AP ≥ 4 muestra el desempeño apropiado del modelo
matemático en la predicción de valores de respuesta. De la mesa5,PA =132.936indica la capacidad del
modelo final [Ec. (13)].

Además de los ítems mencionados, el desempeño gráfico del modelo final [Ec. (13)] se muestra en la Fig.7 proporcionando un
diagrama de dispersión de la medidaηsdatos versus los datos pronosticados. Este gráfico indica una concordancia adecuada entre
los resultados observados.ηsvalores y estimados por la Ec. (13).

Efecto de los principales factores sobreηs.De acuerdo apag-índice estadístico de valor presentado en la tabla5, los efectos de
fr. yL/DenηsLas variables se muestran en el material complementario (ver Fig. S2). De la Fig. S2a, variaciones deηsversus números
de Froudefr.tuvo una tendencia cuadrática a la baja paral/D= 13. La figura S2a ilustra que el efecto del término de segundo orden
de los números de Froudefr.sobre elηsvariable es considerable. Esta conclusión se puede obtener de ANOVA (Tabla5). Elpag
-valores defr.2término es menor que 0.0001, lo que demuestra que elfr.2El término es significativo en el modelo final. Aunque
según la tabla5, el efecto de los términos de primer y segundo orden del/DEl factor en los valores de respuesta es
estadísticamente significativo. La figura S2b muestra que los cambios en los valores de respuesta están más influenciados por el
término de primer orden del/Dfactor. Además del elemento mencionado, la Fig. S2b muestra que el valor deηsaumenta con un
aumento enl/Dvalor.

Desempeño cualitativo del modelo matemático.En esta sección, el desempeño cualitativo

del modelo final [Ec. (13)] al predecirηsvalores para diferentes niveles defr.yl/DSe investigan los factores. Cifra8a–d ilustra la
tendencia de las variaciones en los valores medidos y predichos deηsversus ell/Dfactor para fr.=1,77, 2,01, 2,18 y 2,32. Más
concretamente, la Fig.8ilustra que las variaciones deηscontral/Den niveles dispares defr.El número tuvo una tendencia ascendente.
Por ejemplo, a nivel defr.=1.77, figura.8un demostrado queηsLos valores aumentaron suavemente, pasando de aproximadamente
el 33% enl/D= 10 a aproximadamente 62% enl/D= 16. Al mismo nivel del/Dvalores, la eficiencia de disipación de energía en el eje
disminuyó con un aumento en los valores de Fr. Como ejemplo notable, paral/D= 13,ηslos valores aumentaron del 50% enfr.=1,77 a
aproximadamente el 45% enfr.=2.01 como se indica en la Fig.8a, b, respectivamente. Entonces,ηsse desplomó a aproximadamente
el 38% y el 20% enfr.=2.18 y 2.32, como se ve en la Fig.8c, d, respectivamente. En este estudio, el polinomio de segundo orden se
ajustó a los datos experimentales de la Fig.8.
Cifra9indicó la robustez de la ecuación. (13) para describir variaciones deηsvalores versusfr.número para todos los niveles del/
D. Como se ve en la Fig.9, cambios enηsvalores en contrafr.Los valores siguen una tendencia a la baja. Por ejemplo, el desempeño
de la Ec. (13) indicó que, porl/D=16, lo previstoηslos valores disminuyeron de aproximadamente el 62% en fr.=1,77 a alrededor del
38% enfr.=2.32. Además, al mismo nivel defr.número,ηslos valores aumentaron con un aumento enl/D. Como ejemplo, parafr.=
2,32,ηslos valores aumentan desde aproximadamente el 12% enl/D=10 a aproximadamente 38% en l/D=dieciséis.

Según las Figs.8y9, la diferencia absoluta máxima deηsvalores (diferencia máxima entre los valores medidos y
predichos) con una cantidad insignificante del 2,21% muestra que la ecuación. (13) tiene un desempeño
aceptable en la predicciónηsvalores.

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 10

Vol:.(1234567890)
www.nature.com/scientificreports/

70 70

60 (a) 60 (b)

50 50

ηs(%)

ηs(%)
40 40

30 30
Fr = 1,77 (medido)
20 Fr = 1,77 [Predicho por la ecuación. 20 Fr = 2,01 (medido)
Fr = 2,01 [Predicho por la ecuación.
(13)] Poli. (Fr = 1,77 (Medido))
(13)] Poli. (Fr = 2,01 (medido))
10 10
7 10 13 dieciséis 19 7 10 13 dieciséis 19
L/D L/D

70 50
Fr = 2,18 (medido) Fr = 2,32 (medido)
Fr = 2,18 [Predicho por la ecuación. Fr = 2,32 [Predicho por la ecuación.
60 (13)] Poli. (Fr = 2,32 (Medido))
(13)] Poli. (Fr = 2,18 (Medido))
40
50
30
ηs(%)

40 ηs(%)
20
30

20 (C) 10
(d)
10 0
3 6 9 12 15 18 7 10 13 dieciséis 19
L/D L/D
Figura 8.variaciones deηsvalores versusl/Dfactor para varios valores defr.factor: (a)fr.=1.77, (b)fr.=2.01, (c) fr.=
2.18, y (d)fr.=2.32.

Comparación de los resultados actuales con la literatura.En esta sección, los resultados de los experimentos.
Las investigaciones se han comparado con la literatura relevante. El presente estudio empleó una amplia gama de descarga de
flujo, número de Froude yl/Dproporción, en comparación con las investigaciones anteriores. Como se ve en la tabla6, muchos
estudios (por ejemplo, Vischer y Hager14, Jain y Kennedy15) realizó experimentos con los valores sin cambios del/Dratios y, en
consecuencia, no se extrajo inevitablemente ningún modelo matemático debido al rango límite de parámetros efectivos. Esta
investigación proporcionó una ecuación empírica que incluyel/Dyfr.. Además,l/Drelación utilizada por Zhao et al.19La investigación
estuvo en el rango de la presente investigación.l/D(10-16). Por cierto, la eficiencia energética obtenida a partir de los resultados de
Zhao et al.19experimentos fue del 34%, que estaba en el rango del presente estudio (10,80–62,29%). Como mérito importante en
comparación con la literatura, los resultados actuales se pueden generalizar para rangos adicionales de factores efectivos (es
decir,fr.yL/D) que no se observan en este estudio. De esta manera, la interacción existente entre los factores efectivos y la
eficiencia de la pérdida de energía se puede obtener sin la realización de experimentos cuando los expertos aplicaron factores
efectivos cuyos límites están dentro del rango de los parámetros de prueba actuales.

Consentimiento para publicar.Todos los autores otorgan a la Editorial el permiso de los autores para publicar el trabajo de investigación.

Conclusión
Este estudio investigó el modelo experimental de la estructura de vórtices de aguas residuales urbanas. Se ha estudiado la pérdida de
energía del flujo en el eje vertical para proporcionar un modelo matemático con el fin de investigar varios efectos de los parámetros
contribuyentes (fr.yl⁄D) sobre el comportamiento del flujo en el eje vertical de la estructura de vórtice. De esta manera, los principales
hallazgos de esta investigación se pueden resumir de la siguiente manera:

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 11

Vol.:(0123456789)
 www.nature.com/scientificreports/

60
40
(b)
(a)
35 50

30
40
25

ηs(%)
ηs(%)

20 30

15
L/D = 10 (Medido) 20 L/D = 13 (Medido)
10 L/D = 10 [Predicho por la ecuación. L/D = 13 [Predicho por la ecuación.
(13)] Poli. (L/D = 10 (Medido)) (13)] Poli. (L/D = 13 (Medido))
5 10
1.5 1.7 1.9 2.1 2.3 2.5 1.5 1.7 1.9 2.1 2.3 2.5
fr. fr.

70
(C)
60

50
ηs(%)

40

30 L/D = 16 (Medido)
L/D = 16 [Predicho por la ecuación.
(13)] Poli. (L/D = 16 (Medido))
20
1.5 1.7 1.9 2.1 2.3 2.5
fr.

Figura 9.variaciones deηsvalores versusfr.factor para varios valores del/Dfactor: (a)l/D= 10, (b)l/D= 13, y (C
)l/D= 16.

Autores q(metro3/s) l/D ηS(%)

Vischer y Hager14 140 50 85
Jain y Kennedy15 90 100 90
Jeanpierre y Lachal18 88 9 62
Zhao y cols.19 0.051 14 34*
Este estudio 0,0097–0,0271 10-16 10,80–62,29

Tabla 6. Comparación de los resultados actuales con la literatura relevante. *Para la mitad de la longitud de un eje de caída con
l=14D.

• Paraq≥ 1,7qdEn este caso, el núcleo de aire formado al principio del eje vertical se destruyó y, aunque se perturba la
circulación del aire, la superficie libre del flujo en el canal de entrada es casi horizontal y se perturba la formación de
un flujo de vórtice.
• Los valores de espesor relativo del flujo de vórtice a lo largo del eje vertical para todosl/DLos niveles de los factores comienzan con una
tendencia decreciente relativamente significativa y alcanzan un valor casi constante en elevaciones bajas.

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 12

Vol:.(1234567890)
www.nature.com/scientificreports/

• Reducir los valores de velocidad vertical para todos los niveles del/DEl factor en niveles cercanos a la cámara de
disipación hace que el flujo de gotas ingrese a esta sección con menos impacto y mejora la eficiencia de disipación de
energía en esta sección.
• Valores muy cercanos de velocidades tangenciales en diferentes niveles del eje vertical para todosl/DLos niveles de factor indicaron que
el flujo del vórtice se mantuvo estable a lo largo del eje vertical.
• El análisis de varianza mostró que la interacción defr.yl/Dfactores sobre los valores de la eficiencia de disipación de
energía del flujo (ηs) no fue estadísticamente significativo. Por lo tanto, se puede decir que la tendencia en las
variaciones de los valores delηsversus un factor (ofr.yl/D) en diferentes niveles de otro factor se mantuvo sin cambios.
• Los resultados mostraron que los efectos de los términos de primer y segundo orden del/DLa relación entre los valores de respuesta fue
significativa. También se demostró que el efecto del término de primer orden sobre los cambios en los valores de respuesta era más
importante y la tendencia de los cambios no tenía una tendencia determinada.
• Variaciones de los valores de respuesta versus elfr.El factor se vio afectado principalmente por el término de segundo orden de este
factor y la tendencia de cambios en los valores de respuesta siguió una forma cuadrática.
• La eficiencia de disipación de energía del flujo en el eje vertical (ηs) se desplomó cuandofr.tuvo una tendencia creciente. A
diferencia de,ηsyl/Dla relación tenía dependencia inversa.
• La eficiencia de disipación de energía del flujo en el eje vertical (ηs) se obtuvo entre 10,80 y 62,29%.
• Se ha propuesto una ecuación polinómica cuadrática para comprender el significado físico de laηs
valores en contrafr.yl/Djunto conR2=0,994. La ecuación propuesta fue confiable parafr.=1,77–2,32 y l/D
=10–16.
• El nivel dez/D= 5,9 se introdujo como una sección transversal en la que se observó la tendencia de las variaciones en las propiedades del
flujo del vórtice (es decir, velocidad tangencial, velocidad vertical y espesor relativo del flujo) en el eje vertical.

Disponibilidad de datos
Los datos no están disponibles públicamente debido a restricciones, como que contienen información que podría comprometer la privacidad de los
participantes de la investigación. Póngase en contacto con el autor correspondiente para solicitar datos.

Recibido: 19 de julio de 2022; Aceptado: 24 de enero de 2023

Referencias
1. Jain, SC, Kennedy, JF (1983). Estructuras de caída de flujo de vórtice para el sistema de almacenamiento en línea del Distrito Metropolitano de Alcantarillado de
Milwaukee.IIDH Rep. No. 264, Univ. de Iowa, Iowa City, Iowa.
2. Padulano, R. & Del Giudice, G. Tuberías de drenaje y rebosadero verticales: revisión de la literatura y nuevos datos experimentales.J. Irrig. Drenar. Ing.
144(6), 04018010 (2018).
3. Fernandes, J. & Jónatas, R. Caracterización experimental del flujo en un eje de caída de vórtice en espiral.Ciencia del agua. Tecnología.80(2), 274–281
(2019).
4. Granata, F. Cascadas Dropshaft en sistemas de drenaje urbano.Ciencia del agua. Tecnología.73(9), 2052-2059 (2016).
5. Del Giudice, G., Gisonni, C. y Rasulo, G. Diseño de una entrada de vórtice de espiral para flujo de aproximación supercrítico.J. Hidráulico. Ing.136(10),
837–841 (2010).
6. Ma, Y., Zhu, DZ y Rajaratnam, N. Arrastre de aire en un eje de caída de flujo alto y hundido.J. Hidráulico. Ing.142(10), 04016038 (2016).
7. Liu, ZPet al.Investigación experimental y numérica del flujo en un aliviadero de eje de caída de vórtice recientemente desarrollado.J. Hidráulico. Ing.
144(5), 04018044 (2018).
8. Jain, SC )1984(. Entrada de vórtice tangencial.J. Hidráulico. Ing., 110)12(, 1693–1699.
9. Pfister, M., Crispino, G., Fuchsmann, T., Ribi, JM y Gisonni, C. Múltiples ramas de flujo de entrada en un eje de caída de vórtice de tipo supercrítico.
J. Hidráulico. Ing.144(11), 05018008 (2018).
10. Hager, WH Entrada de caída Vortex para flujo de aproximación supercrítico.J. Hidráulico. Ing.116(8), 1048-1054 (1990).
11. Crispino, G., Pfister, M. y Gisonni, C. Aspectos de diseño hidráulico para flujo supercrítico en ejes de caída de vórtice.Agua Urbana J.dieciséis(3), 225–
234 (2019).
12. Chan, SN, Guo, J. & Lee, JHW Modelado físico y numérico del flujo giratorio en una entrada de vórtice de espiral.J. Hidroambiente. Res.40,
64–76 (2022).
13. Jain, SC, Ettema, R. (1987). Tomas de flujo de vórtice.Manual de diseño de estructuras hidráulicas del IAHR, vol. 1, AA Balkema, Róterdam,
Países Bajos.
14. Vischer, DL, Hager, WH (1995). Gotas de vórtice.Disipadores de energía: manual de diseño de estructuras hidráulicas., núm. 9, cap. 9, AA
Balkema, Rotterdam, Países Bajos, 167–181.
15. Jain, SC, Kennedy, JF (1984). Caída de flujo de vórtice.Proc., 1984 Int. Síntoma. sobre Hidrología Urbana, Hidráulica y Control de Sedimentos, Univ. de
Kentucky, Lexington, Ky., 115–120.
16. Dong, X., Gao, J. (1995). Informe sobre el estudio modelo de modernización de un túnel de desvío en un aliviadero de eje de descenso de vórtice en Shapai
PowerStation.Representante de investigación del IWHR., Instituto Chino de Investigación sobre Recursos Hídricos y Energía Hidroeléctrica, China (en chino).
17. Zhao, CH, Sun, SK y Liu, ZP Estudio óptimo sobre la profundidad del pozo tranquilizador para un túnel de liberación de inundaciones con eje de flujo de rotación.El
poder del agua2001(5), 30–33 (2001) ((en chino)).
18. Jeanpierre, D. & Lachal, A. Dissipation d'énergie dans un puits a vortex.Houille Blanche21(7), 825–831 (1966).
19. Zhao, CH, Zhu, DZ, Sun, SK y Liu, ZP Estudio experimental del flujo en un eje de caída de vórtice.J. Hidráulico. Ing.132(1), 61–68 (2006).

20. Crispino, C., Contestabile, P., Vicinanza, D. y Gisonni, C. Disipación de carga de energía y presiones de flujo en ejes de caída de vórtice.Agua 13(2), 165
(2021).
21. Chan, SN y Chiu, CH Regímenes de flujo de un sistema de túnel de eje de caída sobrecargado.J. Irrig. Drenar. Ing.147(10), 04021045
(2021).
22. Singh, TS, Rajak, U., Samuel, OD, Chaurasiya, PK, Natarajan, K., Verma, TN, Nashine, P. (2021). Optimización de los parámetros de rendimiento y
emisiones de un motor diésel de inyección directa alimentado con microalga Spirulina (L.): metodología de superficie de respuesta y enfoque del
método factorial completo.Combustible,285, 119103.
23. Mahmoudi-Rad, M. y Khanjani, MJ Disipación de energía del flujo en la estructura del vórtice: investigación experimental.J. Sistema de tuberías. Ing.
Practica.10(4), 04019027 (2019).
24. Mahmoudi-Rad, M. & Najafzadeh, M. Papel de la cámara de disipación en la pérdida de energía de estructuras de vórtice: evaluación experimental.
Fluir. Medidas. Instrumento.88, 102232 (2022).

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 13

Vol.:(0123456789)
 www.nature.com/scientificreports/

25. Hager, WH (2010).Hidráulica de aguas residuales: teoría y práctica. Springer, Nueva York.
26. Toda, K. & Inoue, K. Diseño hidráulico de estructuras de toma de sistemas de túneles subterráneos ubicados profundamente.Ciencia del agua.
Tecnología. 39(9), 137-144 (1999).
27. Yu, D. & Lee, J. Hidráulica de toma de vórtice tangencial para drenaje urbano.J. Hidráulico. Ing.135(3), 164-174 (2009).
28. Mulligan, S., Casserly, J. y Sherlock, R. Efectos de la geometría sobre fuertes vórtices de superficie libre en flujos de aproximación subcríticos.J.
Hidráulico. Ing.142(11), 04016051 (2016).
29. Montgomery, DC (2017).Diseño y análisis de experimentos, 6ª edición, Wiley.
30. Sangsefidi, Y., Mehraein, M., Ghodsian, M. & Motalebizadeh, MR Evaluación y análisis del flujo sobre vertederos en arco utilizando metodologías
tradicionales y de superficie de respuesta.J. Hidráulico. Ing.143(11), 04017048 (2017).

Contribuciones de autor
MM-R.; Preparar trabajos experimentales; análisis dimensional y redacción de resultados de estudios experimentales; MINNESOTA;
Escritura: preparación del borrador original. Escritura: revisión y edición.

Conflicto de intereses
Los autores declaran no tener conflictos de intereses.

Información adicional
Información suplementariaLa versión en línea contiene material complementario disponible enhttps://doi.org/10.1038/
s41598-023-28762-2.

Correspondenciay las solicitudes de materiales deben dirigirse a MM-R. Información

de reimpresiones y permisos.está disponible enwww.nature.com/reprints.

nota del editorSpringer Nature se mantiene neutral con respecto a reclamos jurisdiccionales en mapas publicados y
afiliaciones institucionales.

Acceso abiertoEste artículo está bajo una licencia internacional Creative Commons Attribution 4.0,
que permite su uso, intercambio, adaptación, distribución y reproducción en cualquier medio o
formato, siempre y cuando dé el crédito apropiado a los autores originales y a la fuente, proporcione un enlace a la
licencia Creative Commons e indique si se realizaron cambios. Las imágenes u otro material de terceros en este artículo
están incluidos en la licencia Creative Commons del artículo, a menos que se indique lo contrario en una línea de crédito
al material. Si el material no está incluido en la licencia Creative Commons del artículo y su uso previsto no está permitido
por la normativa legal o excede el uso permitido, deberá obtener permiso directamente del titular de los derechos de
autor. Para ver una copia de esta licencia, visitehttp://creativecommons.org/licenses/by/4.0/.

© El autor (es) 2023

Informes Científicos| (2023) 13:1679 | https://doi.org/10.1038/s41598-023-28762-2 14

Vol:.(1234567890)