UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA CIVIL

DEPARTAMENTO ACADÉMICO DE HIDRAÚLICA E HIDROLOGÍA

LABORATORIO N°3

“SIMULACIÓN DE FLUJOS EN CONTRACCIONES Y

EXPANSIONES”
MECÁNICA DE FLUIDOS II (HH224G)

ESTUDIANTE CÓDIGO

ESPINAL CERRON PAUL JHOBBY 20192014J

DOCENTE:

Ing. Julio Martín Kuroiwa Zevallos

Ing. Belinda Noemi Guillén Vidal

FECHA: 18/06/2025

LIMA-PERÚ
 ÍNDICE

1. INTRODUCCIÓN ................................................................................................................ 3
2. OBJETIVOS ......................................................................................................................... 3
2.1. Objetivo principal........................................................................................................... 3
2.2. Objetivos secundarios ..................................................................................................... 3
3. MARCO TEÓRICO.............................................................................................................. 4
3.1. Ecuación general de flujo variado ................................................................................... 4
4. PROCEDIMIENTO .............................................................................................................. 5
4.1. Modelamiento hidráulico con HEC RAS ........................................................................ 6
4.2. Transiciones del modelado ............................................................................................. 8
4.2.1. Contracción o angostamiento ................................................................................... 8
5. RESULTADOS Y DISCUSIONES ....................................................................................... 9
6. CONCLUSIONES .............................................................................................................. 11

SIMULACIÓN DE FLUJOS EN CONTRACCIONES Y EXPANSIONES

 Resumen
Las transiciones en canales son cambios en la forma geométrica de los mismos, en este informe
se estudiarán los cambios de ancho en una sección rectangular. Una expansión puede servir para
desfogue o evacuación de fluidos, mientras que una reducción puede servir como obra de
encausamiento entre un canal. En el presente informe de laboratorio se busca identificar el
comportamiento de un flujo de 8.4 m3/s cuando se presenta un cambio geométrico de expansión
o reducción en un canal. Durante el desarrollo del laboratorio, se realizó un modelado con el
software HEC RAS para la obtención de las gráficas y resultados del análisis del flujo. Se planteó
un modelado para un canal de 112 m, de contracción, analizando secciones cada 8 m, se
obtuvieron 14 secciones de análisis. Con base en lo anterior, se hizo un análisis de los resultados:
tirantes y números de Froude. Finalmente se corroboró la influencia de la transición con las
mejores características hidráulicas y geométricas en el comportamiento del fluido,
principalmente, con el tirante del flujo.

1. INTRODUCCIÓN
Las transiciones hidráulicas son elementos ampliamente utilizados en canales abiertos para
conectar dos secciones transversales distintas, ya sea por su forma, tamaño o pendiente
longitudinal. Estas transiciones permiten que los cambios en la velocidad del flujo y en el tirante
se produzcan de manera gradual.
Por ejemplo, una transición puede servir para enlazar un tramo de canal con sección rectangular
a otro con sección trapezoidal, o bien para unir dos tramos rectangulares con diferentes anchos.
Estas situaciones pueden darse tanto en flujos subcríticos (lentos) como supercríticos (rápidos).
El presente informe de laboratorio tiene como objetivo realizar un análisis experimental, tanto
cualitativo como cuantitativo, del comportamiento del flujo en transiciones que implican
expansiones y contracciones de canal. Para ello, se ha empleado el software HEC-RAS versión
4.1, desarrollado por el Cuerpo de Ingenieros del Ejército de los Estados Unidos (USACE), con
el propósito de modelar hidráulicamente un canal que presenta un ensanchamiento y un
angostamiento.

2. OBJETIVOS
2.1. Objetivo principal
Efectuar el modelado de canales con transiciones (expansión y contracción) utilizando el
programa HEC RAS, para validar la precisión de sus resultados y determinar el comportamiento
de los tirantes en las transiciones.
2.2. Objetivos secundarios
 Calcular los parámetros hidráulicos y geométricos que se producen en las transiciones
en canales.
 Validar los resultados obtenidos con el modelo numérico HEC RAS al contrastar los
resultados obtenidos del diseño hidráulico con los conceptos teóricos.
3. MARCO TEÓRICO
3.1. Ecuación general de flujo variado
Partiremos de la ecuación de energía total, para obtener una expresión de flujo variado que
tenga en consideración la variación del ancho en un canal rectangular.

Donde le símbolo q(x) nos indica que q (caudal unitario) es función de x, es decir, el caudal
unitario dependerá de la posición de la sección a analizar en nuestro informe, entonces:

Además, el término dH/dx es el cambio de energía en relación con la distancia longitudinal o a

la fricción de la pendiente. Se define como:

Mientras que el término dz/dx es el cambio de elevación del fondo del canal con respecto a la
distancia o a la pendiente del fondo. Se define como:

Cuando 𝑄=𝑞(𝑥)𝑏 es un valor constante como en nuestra simulación, tendremos:

Reemplazando, la ecuación se simplifica de la siguiente forma, la llamaremos ecuación general:

Esta ecuación se puede subdividir en dos casos, dependiendo de la variación del ancho tenemos:
3.1.1. CASO 1: Canal con ancho constante
En este caso, el término db/dx es igual a cero debido a que no varía el ancho de la sección
transversal del canal. Además, considerando muy pequeña la pendiente del canal (dH/dx≈0),
esto solo nos servirá para estimar el comportamiento de los tirantes, mas no la obtención del
valor exacto del mismo.
La ecuación se reduce a:
En este momento puede considerarse una serie de subcasos.
a. Si dz/dx > 0 y F < 1, entonces (1 - F2) > 0 y dy/dx debe ser menor que cero. Por lo
tanto, el tirante del flujo se decrementa al incrementarse x.
b. Si dz/dx > 0 y F > 1, entonces (1 - F2) < 0 y dy/dx debe ser mayor que cero. Por lo
tanto, bajo estas condiciones el tirante del flujo se incrementa al incrementarse x.
c. Si dz/dx < 0 y F < 1, entonces (1 - F2) > 0 y dy/dx debe ser mayor que cero. Por lo
tanto, bajo estas condiciones, el tirante del flujo se incrementa al incrementarse x.
d. Si dz/dx < 0 y F > 1, entonces (1 - F2) < 0 y dy/dx debe ser menor que cero. Bajo estas
condiciones, el tirante del flujo se decrementa al incrementarse x.
3.1.2. CASO 2: Canal con ancho variable
En este caso, el gasto por unidad de ancho no es constante, además, aproximando que la
pendiente del canal es muy pequeña (dH/dx≈0) y que las pérdidas son despreciables (dz/dx≈0),
esto solo nos servirá para estimar el comportamiento de los tirantes, mas no la obtención del
valor exacto del mismo: La ecuación general se reduce a:

De esta última ecuación, podemos considerar 4 subcasos:

a) Si db/dx > 0 y F < 1, entonces (1 - F2) > 0 y dy/dx debe ser mayor que cero. En este
caso el tirante del flujo se incrementa al incrementarse x.
b) Si db/dx > 0 y F > 1, entonces (1 - F2) < 0 y dy/dx debe ser menor que cero. En este
caso el tirante del flujo disminuye al incrementarse x.
c) Si db/dx < 0 y F < 1, entonces (1 - F2) > 0 y dy/dx debe ser menor que cero. Bajo estas
condiciones, el tirante del flujo disminuye al incrementarse x.
d) Si db/dx < 0 y F > 1, entonces (1 - F2) < 0 y dy/dx debe ser mayor que cero. Bajo estas
condiciones, el tirante del flujo se incrementa al incrementarse x.

4. PROCEDIMIENTO
Se realizó simulaciones de flujo en canales abiertos con el programa HEC RAS 4.1. Los datos
son de acuerdo con el apellido paterno y materno.
Datos:
Apellido paterno: Espinal
Tabla 1: Datos correspondiente al apellido paterno

Apellido paterno Q (m3/s) W (m) H (m) n

D-K 8.4 3.2 3.2 0.012

Apellido materno: Cerron

Tabla 2: Datos correspondientes al apellido materno

Apellido materno So
A-C 0.0005
4.1. Modelamiento hidráulico con HEC RAS
Para el modelamiento del canal con el programa HEC-RAS se ha seguido el siguiente
procedimiento:
a) Definido el trazo del canal, sobre el perfil longitudinal se traza la rasante del canal, siguiendo
el procedimiento para el dimensionamiento del canal antes indicado.
b) Teniendo las cotas de la rasante del canal se define en el programa HEC-RAS la sección de
cada tramo de canal.

Figura 1: Trazo del canal HEC RAS

c) Cada sección en el programa se define como una estación, que requiere los siguientes datos:
las cotas en las secciones, la longitud a la que se encuentra la siguiente estación aguas abajo, el
valor de Manning para la sección y los límites del flujo en la sección.

Figura 2: Geometría del canal de la sección 1

d) Definida las estaciones se ingresa los datos hidráulicos (Steady Flow Data) del flujo
permanente, y agregamos la condición de borde. La condición de borde es que en la sección
aguas abajo (DownStream) la pendiente de la energía es igual a la del canal, es decir, empieza
como un flujo uniforme.
 Figura 3: Ingreso del caudal

Figura 4: Condiciones de borde

e) Definido las características e ingresado todos los tramos del canal en el programa se procede
a correr (Run and Compute) la información para obtener el comportamiento del flujo en el canal
de caudal 8.4 m3/s.

Figura 5: Ejecución del modelado en HEC RAS

4.2. Transiciones del modelado
4.2.1. Contracción o angostamiento
Primero se simulará una contracción en la que un flujo pasará a través de un canal rectangular de
ancho nominal de 3.2 m, que transporta un caudal 8.4 m3/s. El primer tramo aguas abajo tiene
una longitud igual a 32 m, el segundo tramo (desde el extremo ubicado aguas abajo) es una
transición en la que el ancho pasa de 3.2 m a 2.4 m y tiene una longitud igual a 16 m. El tercer
tramo es un tramo recto cuyo ancho es 2.4 m y tiene una longitud igual a 16 m Aguas arriba de
este tramo existe un tramo de transición en que el ancho pasa de 3.2m a 2.4 m. Aguas arriba de
este tramo existe un tramo recto cuya longitud es igual 32 m, y cuyo ancho es el ancho nominal.
Considerando que las secciones están espaciadas cada 8 metros y que la altura del canal es 3.2m.

Figura 6: Esquema del modelo de contracción

5. RESULTADOS Y DISCUSIONES
5.1. Contracción y angostamiento
De la ecuación simplificada de canal con sección de ancho variable (CASO II), teniendo en
referencia el subcaso 3 (dz/dx < 0), se puede observar que a partir de los resultados del
modelamiento el régimen del flujo es subcrítico (Fr < 1) en todo el recorrido del canal, entonces
el tirante en la zona de contracción (Sección 80 a 64) disminuya desde y80 =1.7108 (tirante
máximo) tal y como se observa en la tabla 3.
Tabla 3: Resumen de perfil de flujo del modelado de contracción o ensanchamiento
E.G.
Q Total Nivel de E.G. Área Ancho
X (m) Z (m) y (m) Slope V (m/s) Froude
(m3/s) agua (m) Elev (m) (m2) (m)
(m/m)
112 8.4 100.06 101.7648 1.7087 101.8851 0.000438 1.5363 5.4676 3.1998 0.3752
104 8.4 100.05 101.7613 1.7092 101.8815 0.000438 1.5359 5.469 3.1998 0.375
96 8.4 100.05 101.7578 1.71 101.8779 0.000438 1.5352 5.4715 3.1998 0.3748
88 8.4 100.04 101.7543 1.7104 101.8743 0.000437 1.5348 5.473 3.1998 0.3746
80 8.4 100.04 101.7508 1.7108 101.8707 0.000437 1.5344 5.4743 3.1998 0.3745
72 8.4 100.04 101.6952 1.6592 101.8617 0.00068 1.8082 4.6456 2.7999 0.4481
64 8.4 100.03 101.5872 1.5552 101.8452 0.001226 2.2506 3.7324 2.4 0.5761
56 8.4 100.03 101.5751 1.547 101.8358 0.001243 2.2624 3.7128 2.4 0.5806
48 8.4 100.02 101.5626 1.5386 101.8263 0.001261 2.2749 3.6925 2.4 0.5854
40 8.4 100.02 101.6129 1.5928 101.7937 0.000756 1.8836 4.4597 2.7999 0.4764
32 8.4 100.02 101.6416 1.6257 101.7744 0.0005 1.6148 5.2019 3.1998 0.4043
24 8.4 100.01 101.6376 1.6257 101.7704 0.0005 1.6148 5.2018 3.1998 0.4043
16 8.4 100.01 101.6335 1.6256 101.7664 0.0005 1.6149 5.2016 3.1998 0.4043
8 8.4 100 101.6295 1.6256 101.7624 0.0005 1.6149 5.2015 3.1998 0.4043
0 8.4 100 101.6255 1.6255 101.7584 0.0005 1.615 5.2014 3.1998 0.4043

5.1.1. Perfil de flujo en contracción

De la figura 7, se puede apreciar el perfil de flujo de canal con contracción en la sección
intermedia.

Figura 7: Perfil flujo en contracción

5.1.2. Tirante en contracción mediante fórmula general
Las estaciones desde 80m a 64m (Secciones 64, 72 y 80) pertenecen a la zona de contracción
(angostamiento). El ancho se reduce de 3.2 m a 2.4 m (db/dx<0)
Fórmula general

Cálculo de tirante en la sección 9 (𝒚𝟕𝟐)

Reemplazando los valores obtenidos de la sección 10 (y80 ) en el HEC RAS en la fórmula, para
calcular el tirante en la sección 9 (y72 )
(𝑦72 − 𝑦80 ) 𝑦72 (𝑏72 − 𝑏80 )
(1 − 𝐹𝑟72 2 ) − 𝐹𝑟72 2 ( ) − 𝑆𝑜 = −𝑆𝑓72
8 𝑏72 8
(𝑦72 − 1.7108) 𝑦72 (2.8 − 3.2)
(1 − 0.44812 ) − 0.44812 ( ) − 0.0005 = −0.00068
8 2.8 8
𝑦72 (𝑡𝑒ó𝑟𝑖𝑐𝑜 ) = 1.650 𝑚
Entonces el error obtenido del valor teórico con respecto al HEC RAS es de:
|1.650 − 1.6592|
𝑒𝑟𝑟𝑜𝑟 = × 100%
1.6592
𝑒𝑟𝑟𝑜𝑟 = 0.554%
Valor aceptable

Cálculo de tirante en la sección 8 (𝒚𝟔𝟒)

Reemplazando los valores obtenidos de la sección 9 (y72 ) en el HEC RAS en la fórmula, para
calcular el tirante en la sección 8 (y64 )
(𝑦64 − 𝑦72 ) 𝑦64 (𝑏64 − 𝑏72 )
(1 − 𝐹𝑟64 2 ) − 𝐹𝑟64 2 ( ) − 𝑆𝑜 = −𝑆𝑓64
8 𝑏64 8
(𝑦64 − 1.6592) 𝑦64 (2.4 − 2.8)
(1 − 0.57612 ) − 0.57612 ( ) − 0.0005 = −0.001226
8 2.4 8
𝑦64 (𝑡𝑒ó𝑟𝑖𝑐𝑜 ) = 1.524 𝑚
Entonces el error obtenido del valor teórico con respecto al HEC RAS es de:
|1.524 − 1.5552|
𝑒𝑟𝑟𝑜𝑟 = × 100%
1.5552
𝑒𝑟𝑟𝑜𝑟 = 2.006%
Valor aceptable
6. CONCLUSIONES
 Se confirmó que en la zona de contracción (secciones 80m a 64m), el tirante
disminuye progresivamente desde 1.7108 m hasta 1.5552 m, cumpliendo con
la teoría para flujos subcríticos (Fr < 1) donde db/dx < 0, validando el subcaso 3
de la ecuación general.
 Durante todo el recorrido del canal de 112 m, el flujo se mantuvo en régimen
subcrítico (números de Froude entre 0.3745 y 0.5854)
 Los errores obtenidos al comparar los valores teóricos con HEC-RAS fueron:
Sección 9 (72m): 0.554% de error y en la Sección 8 (64m): 2.006% de error.
Estos valores confirman la validez de la ecuación general de flujo variado para
el análisis de transiciones.