Universidad Nacional Mayor de San Marcos

Universidad del Perú. Decana de América

Facultad de Ciencias Físicas
Escuela Profesional de Ingeniería Mecánica de Fluidos

Simulación numérica de la dinámica del flujo

secundario helicoidal y del transporte de sedimentos en
el modelo hidráulico del río Madre de Dios – tramo La
Pastora

TESIS

Para optar el Título Profesional de Ingeniero Mecánico de Fluidos

AUTOR
Magaly CUSIPUMA AYUQUE

ASESOR
Dr. Ing. Carlos Augusto Antonio CARBONEL HUAMAN

Lima, Perú

2020
Reconocimiento - No Comercial - Compartir Igual - Sin restricciones adicionales

https://creativecommons.org/licenses/by-nc-sa/4.0/
Usted puede distribuir, remezclar, retocar, y crear a partir del documento original de modo no
comercial, siempre y cuando se dé crédito al autor del documento y se licencien las nuevas
creaciones bajo las mismas condiciones. No se permite aplicar términos legales o medidas
tecnológicas que restrinjan legalmente a otros a hacer cualquier cosa que permita esta licencia.
Referencia bibliográfica

Cusipuma, M. (2020). Simulación numérica de la dinámica del flujo secundario

helicoidal y del transporte de sedimentos en el modelo hidráulico del río Madre de
Dios – tramo La Pastora. Tesis para optar el título de Ingeniero Mecánico de Fluidos.
Escuela Profesional de de Ingeniería Mecánica de Fluidos, Facultad de Ciencias
Físicas, Universidad Nacional Mayor de San Marcos, Lima, Perú.
 Hoja de metadatos complementarios

Código ORCID del autor 0000-0002-9855-1016

DNI o pasaporte del autor 70312459

Código ORCID del asesor 0000-0001-9465-8118

DNI o pasaporte del asesor 08772534

Grupo de investigación “—“

Agencia financiadora Financiamiento propio.

Lima, Escuela de Ingeniería Mecánica de

Ubicación geográfica donde se Fluidos.
desarrolló la investigación 12°03'24.4"S 77°05'13.7"W

Año o rango de años en que

2017 - 2019
se realizó la investigación

Ingeniería civil
http://purl.org/pe-repo/ocde/ford#2.01.01

Ingeniería estructural y municipal

Disciplinas OCDE http://purl.org/pe-repo/ocde/ford#2.01.04

Ingeniería del Petróleo, Energía, Combustibles

http://purl.org/pe-repo/ocde/ford#2.07.03
 Universidad Nacional Mayor de San Marcos
(Universidad del PERÚ, Decana de América)

FACULTAD DE CIENCIAS FÍSICAS

ESCUELA PROFESIONAL DE INGENIERÍA MECÁNICA DE FLUIDOS

ACTA DE SUSTENTACIÓN DE TESIS EN LA MODALIDAD VIRTUAL

PARA OPTAR EL TITULO PROFESIONAL DE INGENIERO MECANICO DE
FLUIDOS

Siendo las 11:00 a.m del viernes 4 de diciembre de 2020, en la Sala de Sesión Virtual de la
Facultad de Ciencias Físicas, bajo la Presidencia del Ing. GUIDO AMÉRICO ROZAS
OLIVERA (Presidente), MSc. Ing. DOUGLAS DONAL SARANGO JULCA
(Miembro), Dr. Ing. JORGE LUIS GASTELO VILLANUEVA (Miembro) y el Dr. Ing.
CARLOS AUGUSTO ANTONIO CARBONEL HUAMÁN (Asesor), se dio inicio a la
Sesión Pública Virtual de Sustentación de Tesis para la obtención del Título Profesional de
Ingeniero Mecánico de Fluidos, mediante la Modalidad M1 – Por Sustentación y
aprobación de Tesis Profesional, de la Bachiller:

MAGALY CUSIPUMA AYUQUE

El Presidente del Jurado en primer lugar dio lectura al documento formal que designa al
Jurado y aprueba el acto de sustentación de la Tesis, así como del Resumen del Expediente
e invitó a la citada Bachiller a realizar la presentación y exposición de su Tesis titulada:
“SIMULACIÓN NUMÉRICA DE LA DINÁMICA DEL FLUJO SECUNDARIO
HELICOIDAL Y DEL TRANSPORTE DE SEDIMENTO EN EL MODELO
HIDRÁULICO DEL RÍO MADRE DE DIOS – TRAMO LA PASTORA”.

Concluida la exposición de la candidata y luego de las preguntas de rigor por parte del
Jurado, el Presidente del Jurado con el apoyo de la Unidad de Informática, invitó a la
Bachiller y al público a abandonar momentáneamente la Sala de Sesión Virtual, para dar
paso a la deliberación y calificación por parte del Jurado.

Al término de la deliberación el Presidente de Jurado invitó a la candidata y al público en

general a pasar a la Sala de Sesión Virtual, para dar lectura al Acta de Calificación, en cuyo
contenido se ha registrado la nota obtenida por la bachiller, la misma que ha sido de:

17 Diecisiete
..................................... .....................................
(NÚMERO) (LETRAS)

Habiendo concluido la Sustentación de la Tesis, el Presidente del Jurado remitirá el Acta al

Vicedecano Académico de la Facultad de Ciencias Físicas, quien la suscribirá dando fe de
la realización del acto virtual de manera remota.

Ciudad Universitaria, Pabellón de Ciencias Físicas, Ca. Germán Amézaga N° 375. Lima, 1
Teléfono: 619-7000 anexo 3806 E-mail: [email protected]
 Universidad Nacional Mayor de San Marcos
(Universidad del PERÚ, Decana de América)

FACULTAD DE CIENCIAS FÍSICAS

ESCUELA PROFESIONAL DE INGENIERÍA MECÁNICA DE FLUIDOS

El Acta de Calificación será remitida al Director de la Escuela Profesional de Ingeniería

Mecánica de Fluidos, quien a su vez solicitará la aprobación del Consejo de Facultad para
el otorgamiento del Título Profesional y su elevación al Consejo Universitario para conferir
el Título Profesional de Ingeniero Mecánico de Fluidos a la Bachiller MAGALY
CUSIPUMA AYUQUE.

Siendo las 13:15 horas, se dio por concluido el acto académico, expidiéndose tres (03)
Actas originales de la Sustentación de la Tesis.

Ciudad Universitaria, 4 de diciembre de 2020

ING. GUIDO AMÉRICO ROZAS OLIVERA MSC. ING. DOUGLAS DONAL SARANGO JULCA
PRESIDENTE DEL JURADO MIEMBRO DEL JURADO

DR. ING. JORGE LUIS GASTELO VILLANUEVA DR. ING. CARLOS A. A. CARBONELHUAMÁN
MIEMBRO DEL JURADO ASESOR DEL JURADO

Firmado digitalmente por LOZANO

BARTRA Whualkuer Enrique FAU
20148092282 soft
Motivo: Soy el autor del documento
Fecha: 09.12.2020 10:24:56 -05:00

DR. WHUALKUER ENRIQUE LOZANO BARTRA

VICEDECANO ACADÉMICO FCF

Datos de la plataforma virtual institucional del acto de sustentación:

https://zoom.us/j/94419018709?pwd=b2l5U0FTWUFLVUJzaG91aHNNVTBkQT09
ID de reunión: 944 1901 8709
Código de acceso: 866353
Grabación archivada en el siguiente enlace:
https://www.facebook.com/FCFUNMSM

Ciudad Universitaria, Pabellón de Ciencias Físicas, Ca. Germán Amézaga N° 375. Lima, 1
Teléfono: 619-7000 anexo 3806 E-mail: [email protected]
Dedicatoria

Agradecer a mis padres por su apoyo incondicional a lo largo de toda mi formación profesional,
no me alcanzara una vida para agradecerles su dedicación y motivación de superarme.
Virgilio e Isabel.
Agradecimientos

Agradecer al "Laboratorío Nacional de Hidráulica - Alfonso Alcedán la Cruz" por

la oportunidad brindada para concluir con mi formación profesional mediante el desarrollo
del trabajo de la presente tesis. Además, agradecer al Ing. Luis Fernando Castro Inga, ac-
tual Investigador del Laboratorio Nacional de Hidráulica-Universidad Nacional de Ingeniería,
quien me motivo para iniciarme en el ámbito de la investigación.

La fase experimental de esta investigación fue financiada mediante aportes del Contrato No
358-PNICP-PIAP-2014, firmado entre la UNI y el Programa Nacional de Innovación para la
Competitividad y Productividad del Ministerio de la Producción, INNOVATE PERU.

En la fase analítica se empleó el cluster de alto rendimiento Tipón, que fue financiado me-
diante el CONVENIO No 257-INNOVATEPERU-EC-2016, “Adquisición de un servidor de
alto rendimiento para investigar la dinámica de los flujos geofísicos y ambientales”, firmado
entre la UNI e INNOVATE PERU.

Agradezco al Dr. Julio Martín Kuroiwa Zevallos, actual Director del Laboratorio Nacional
de Hidráulica por su apoyo y confianza brindada para redactar el árticulo: ’SIMULACIÓN
NUMÉRICA DE LA HIDRODINÁMICA Y TRANSPORTE DE SEDIMENTOS EN EL
MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS - TRAMO LA PASTORA’, que fue
enviado al Congreso Latinoamerico de Hidráulica de Buenos Aires, Septiembre del
2018.
 iii

Igualmente, también a los Dres: Pablo Tassi, Riadh Ata y Pat Prodanovic del Laboratoire
National d’Hydraulique et Environnement LNHE EDF–R&D (Francia), quienes
me permitieron tener una pasantia y por el apoyo en el desarrollo de las simulaciones nu-
méricas. En especial a Pat Prodanovic quien fue un amigo que me brindo su apoyo en mi
estancia en el LNHE, asimismo a Adrien Bourgoin, Sofiane Martel.

Asimismo, agradecer a la Universidad Nacional Mayor de San Marcos y a la Escuela de

Ingeniería Mecánica de Fluidos, donde inicie mi formación profesional. Agradecer al Ing.
Bernardino Salvador quien me guio con el desarrollo de la presente tesis.
Resumen

El río Madre de Dios, en el SE peruano, lleva una alta carga de sedimentos en suspen-
sión y la tasa de migración lateral del cauce en el tramo estudiado es aproximadamente 6
m/año, generando fenómenos de socavación lateral en su margen derecha. El proyecto N 0
358-PNICP-PIAP-2014, auspiciado por INNOVATE, (2014) indica que, en el tramo curvo
cerca de la ciudad de Puerto Maldonado, cuya longitud es aproximadamente 1540 metros, se
está generando el colapso de la ribera dejando desprotegido el talud. Este fenómeno afecta
a los inmuebles, sistemas de abastecimiento de agua y desagüe de la ciudad de Puerto Mal-
donado, y, de no haberse tomado medidas correctivas, hubiera puesto en riesgo más de 100
metros de la Carretera Interoceánica. La erosión lateral ha formado bancos de gran altura
con taludes muy empinados.La verticalidad de los taludes hace que el banco sea más inesta-
ble ante una rápida variación del caudal. Para analizar el problema en 2014 se construyó el
modelo hidráulico en el Laboratorio Nacional de Hidráulica – LNH.
Se construyó un modelo hidráulico de lecho móvil que representa el tramo La Pastora del
río Madre de Dios a escala 1/60 en el Laboratorio Nacional de Hidráulica - LNH. El modelo
hidráulico tiene una longitud de 50 m, un ancho medio, B, de 6.9 m, un radio de curvatura,
R, de 13.8 m y un tirante máximo, y, de 0.25 m. Por lo tanto, la relación radio sobre ancho,
R/B, es 2. El lecho móvil está compuesto de arena cuyo diámetro medio, d50 es 0.4 mm. Antes
del inicio de las pruebas se ejecutó el levantamiento topográfico mediante escaneo láser 3D,
procedimiento que se repitió al día siguiente de ejecutadas las pruebas para poder comparar
la topografía del lecho antes y después de las pruebas. Cuatro bombas sumergibles permiten
probar caudales de hasta 550 L/s que son transportados al modelo a través de tuberías de
PVC de 356 mm (14”).
En la presente investigación se realizó la simulación numérica bidimensional del modelo hi-
 v

dráulico del Río Madre de Dios (Tramo la Pastora) para analizar cualitativamente la dinámi-
ca del flujo secundario y del transporte de sedimentos empleando el módulo TELEMAC-2D
(hidráulico) y SISYPHE (transporte de sedimentos) del sistema de código abierto Open
TELEMAC-MASCARET, desarrollado por el Laboratoire National d’Hydraulique et Envi-
ronnement LNHE EDF–RD (Francia).

En la configuración del modelo numérico se construyó una malla triangular no estructurada

con 12345 elementos y 4578 nodos, con un tamaño medio del elemento de 0.3 m. Las condi-
ciones de contorno que se impusieron son: Q = 0.490 m/s en el límite aguas arriba y el tirante
inicial aguas abajo. Como condición inicial de profundidad de agua cero. La rugosidad que
se prescribió fue de Manning 0.020. El esquema numérico que usa el programa es elementos
finitos, la metodología de turbulencia que se uso es RANS (Reynolds Average Navier Stoke)
y modelo de turbulencia k-epsilon. Además, dos parámetros fueron analizados para usar la
corrección de flujo secundario en el modelo bidimensional que son el coeficiente de produc-
ción de flujo secundario (As) que tiene un valor de 7.071 y el coeficiente de disipación de
flujo secundario (Ads) donde su valor es 0.5, ambos parámetros fueron calibrados. Al usar
ambos parámetros los resultados revelan que representa mejor el comportamiento del flujo
en el tramo curvo que si no se usará.
Para analizar el transporte de sedimentos y evolución del lecho móvil se realizó una simulación
de forma acoplada con el modulo hidrodinámico (TELEMAC-2D) y transporte de sedimentos
(SISYPHE) donde el cálculo de las variable de ambos módulos se obtiene en cada iteración.
En el módulo morfodinamico el transporte de sedimentos fue puramente advectivo y difusivo
también se usó la corrección de pendiente del lecho y corrección de flujo secundario.
Abstract

The Madre de Dios River, in the Peruvian SE, carries a high load of suspended sediments and
the lateral migration rate of the riverbed in the studied section is approximately 6 m year-1,
generating lateral undermining phenomena on its right bank. The INDECI report (2006)
indicates that, in the curved section near the city of Puerto Maldonado, whose length is
approximately 1540 meters, the collapse of the riverbank is being generated leaving the slope
unprotected. This phenomenon affects the buildings, water supply and drainage systems of
the city of Puerto Maldonado, and, had no corrective measures been taken, would have put
more than 100 meters of the Interoceanic Highway at risk. Lateral erosion has formed high-rise
benches with very steep slopes. The verticality of the slopes makes the bank more unstable
in the face of rapid flow variation. To analyze this problem in 2014 built the hydraulic model
at National Laboratory Hydraulic-LNH. A mobile bed hydraulic model was constructed that
represents the La Pastora section of the Madre de Dios river at 1/60 scale at the National
Hydraulics Laboratory - LNH. The hydraulic model has a length of 50 m, an average width,
B, of 6.9 m, a radius of curvature, R, of 13.8 m and a maximum shoulder, and, of 0.25 m.
Therefore, the ratio of radius to width, R / B, is 2. The moving bed is composed of sand who-
se average diameter, d50 is 0.4 mm. The topography was reconsidered using steel rods welded
to a steel plate that was fixed to the bed. Before the start of the tests, the topographic survey
was carried out using 3D laser scanning, a procedure that was repeated the day after the tests
were executed in order to compare the topography of the bed before and after the tests. Four
submersible pumps make it possible to test flow rates of up to 550 L / s that are transported
to the model through 356 mm (14”) PVC pipes. The flows enter two tanks of rectangular
section and at the exit the flows can be measured by rectangular landfills. In each tank it was
possible to measure by measuring the water level using Neyrpic limnimeters. The experience
 vii

in the experimental model lasted 24 hours, where the velocity, depth and bathymetry that
were used to calibrate and validate the numerical model were collected. In the present inves-
tigation, the two-dimensional numerical simulation of the hydraulic model of the Madre de
Dios River (La Pastora Section) was performed to qualitatively analyze the dynamics of the
secondary flow and sediment transport using the TELEMAC-2D and SISYPHE module of
the open source system Open TELEMAC-MASCARET, developed by the Laboratoire Na-
tional d’Hydraulique et Environnement LNHE EDF – R D (France). In the configuration of
the numerical model, an unstructured triangular mesh was constructed with 12345 elements
and 4578 nodes, with an average element size of 0.3 m. The boundary conditions that were
imposed are: Q = 0.490 m / s at the upstream limit and the initial downstream tie. As initial
condition of zero water depth. The roughness that was prescribed was Manning 0.020. The
numerical scheme used by the program is finite elements, the turbulence methodology used
is RANS (Reynolds Average Navier Stoke) and the k-epsilon turbulence model. In addition,
two parameters were analyzed to use the secondary flow correction in the two-dimensional
model which are the secondary flow production coefficient (As) that has a value of 7.071 and
the dissipation coefficient of for secondary flow (Ads) where its value is 0.5 both parameters
were calibrated. When using both parameters, the results reveal that it best represents the
flow behavior in the curved section. To analyze sediment transport and evolution of the mo-
ving bed, it was carried out in conjunction with the hydrodynamic module (TELEMAC-2D)
and sediment transport (SISYPHE) where the calculation of the variables of both modules
is obtained at each iteration. In the morphodynamic module, sediment transport was purely
advective and diffusive. Bed slope correction and secondary flow correction were also used.
Índice general

,. Introducción 1
1.1. Descripción del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2. Formulación del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2.1. Problema general . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2.2. Problemas específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3.1. Objetivo general . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3.2. Objetivos específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4. Justificación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.5. Metodología . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.5.1. Acopio de Información . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5.2. Manejo del software Open TELEMAC . . . . . . . . . . . . . . . . . 10
1.5.3. Simulación numérica computacional . . . . . . . . . . . . . . . . . . . 10
1.5.4. Validación del modelo numérico computacional . . . . . . . . . . . . 11
1.6. Motivación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

,,. Marco teórico 12

2.1. Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.2. Metódos para describir el flujo de fluidos . . . . . . . . . . . . . . . . . . . . 14
2.3. Fuerzas en un elemento de fluido . . . . . . . . . . . . . . . . . . . . . . . . 16
2.4. Ecuaciones diferenciales de flujo de fluidos . . . . . . . . . . . . . . . . . . . 17
2.4.1. Ecuación de conservación de masa . . . . . . . . . . . . . . . . . . . . 17
2.4.2. Ecuación de conservación de momentum . . . . . . . . . . . . . . . . 18
 ÍNDICE GENERAL ix

2.4.3. Ecuación de transporte . . . . . . . . . . . . . . . . . . . . . . . . . . 21

2.4.4. Métodos de solución numérica de las ecuaciones del flujo turbulento . 22
2.4.4.1. Simulación numérica directa (DNS) . . . . . . . . . . . . . . 22
2.4.4.2. Simulación de grandes remolino (LES) . . . . . . . . . . . . 22
2.4.4.3. Metodologia (RANS) . . . . . . . . . . . . . . . . . . . . . . 23
2.5. Transporte de sedimentos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.5.1. Características físicas del sedimento . . . . . . . . . . . . . . . . . . . 26
2.5.2. Fuerzas que actúan en la partícula de sedimento . . . . . . . . . . . . 29
2.5.3. Papel de la rugosidad efectiva . . . . . . . . . . . . . . . . . . . . . . 31
2.5.4. Inicio del movimiento del sedimento . . . . . . . . . . . . . . . . . . . 32
2.5.5. Formas de transporte de sedimentos . . . . . . . . . . . . . . . . . . . 36
2.5.6. Fórmulas de transporte de sedimentos no cohesivo . . . . . . . . . . . 38
2.5.6.1. Meyer-Peter and Muller . . . . . . . . . . . . . . . . . . . . 38
2.5.6.2. Einstein y Brown . . . . . . . . . . . . . . . . . . . . . . . . 39
2.5.6.3. Engelund y Hansen . . . . . . . . . . . . . . . . . . . . . . . 41
2.5.6.4. Van Rijn . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42

,,,. Modelo hidráulico del río Madre de Dios 43

3.1. Similitud entre el prototipo y modelo hidráulico . . . . . . . . . . . . . . . 43
3.2. Descripción del modelo hidráulico del río Madre de Dios - Tramo la Pastora 46
3.2.1. Detalles de replanteo del lecho del río Madre de Dios en el modelo
hidráulico para el Tramo la Pastora . . . . . . . . . . . . . . . . . . . 47
3.2.2. Granulometría del material del lecho . . . . . . . . . . . . . . . . . . 47
3.2.3. Caudal y tasa de transporte . . . . . . . . . . . . . . . . . . . . . . . 50
3.3. Datos obtenidos durante los ensayos . . . . . . . . . . . . . . . . . . . . . . . 50
3.3.1. Batimetría del lecho . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.3.2. Medición de las velocidades . . . . . . . . . . . . . . . . . . . . . . . 53
3.3.2.1. Velocidad media . . . . . . . . . . . . . . . . . . . . . . . . 57
3.3.2.2. Tirantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
 ÍNDICE GENERAL x

,9. Sistema Open TELEMAC MASCARET 62

4.1. Descripción del software . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.2. Características del sistema open TELEMAC MASCARET . . . . . . . . . . 63
4.3. Entorno computacional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
4.4. Modulo hidrodinámico TELEMAC-2D . . . . . . . . . . . . . . . . . . . . . 65
4.5. Modulo de transporte de sedimentos SISYPHE . . . . . . . . . . . . . . . . . 68
4.5.1. Modificación de la magnitud y dirección de la carga del lecho . . . . . 70
4.5.1.1. Efecto de la pendiente del lecho . . . . . . . . . . . . . . . . 70
4.5.1.2. Corrección del esfuerzo cortante en la dirección del flujo se-
cundario . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.5.1.3. La correción del flujo secundarío y la pendiente transversal
del lecho debido a la gravedad . . . . . . . . . . . . . . . . . 74
4.5.1.4. Dirección del transporte de sedimentos . . . . . . . . . . . . 74
4.6. Pre-procesamiento con QGIS y Blue Kenue . . . . . . . . . . . . . . . . . . . 75
4.6.1. Características de la batimetría . . . . . . . . . . . . . . . . . . . . . 76
4.6.2. Modelo digital del terreno (DEM) . . . . . . . . . . . . . . . . . . . . 76
4.6.3. Características de la malla . . . . . . . . . . . . . . . . . . . . . . . . 76
4.6.3.1. Malla bidimensional . . . . . . . . . . . . . . . . . . . . . . 78
4.6.4. Condiciones de frontera . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.6.4.1. Condición de frontera solida . . . . . . . . . . . . . . . . . . 79
4.6.4.2. Condición de frontera liquida . . . . . . . . . . . . . . . . . 79
4.7. Técnicas numéricas usadas por el sistema Open TELEMAC MASCARET . . 80
4.7.1. Esquemas numéricos . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.7.2. Algoritmo de solución . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.8. Post-Procesamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83

9. Aplicación de TELEMAC-2D y SISYPHE en el modelo hidráulico del Río

Madre de Dios - Tramo la Pastora 84
5.1. Delimitación del dominio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
5.1.1. Generación de la malla . . . . . . . . . . . . . . . . . . . . . . . . . . 85
 ÍNDICE GENERAL xi

5.1.2. Generación de los archivos de entrada de TELEMAC-2D . . . . . . . 89

5.1.2.1. Generación del archivo de geometría .slf . . . . . . . . . . . 89
5.1.2.2. Generación del archivo de condiciones de borde .cli . . . . . 90
5.1.3. Archivos necesarios para ejecutar TELEMAC-2D . . . . . . . . . . . 92
5.1.3.1. Descripción de los parámetros del archivo .cas . . . . . . . . 92
5.1.4. Ejecución de TELEMAC-2D . . . . . . . . . . . . . . . . . . . . . . . 95
5.1.5. Configuración del módulo SISYPHE . . . . . . . . . . . . . . . . . . 96

9,. Análisis de resultados de la simulación numérica 99

6.1. Calibración del módulo hidrodinámico . . . . . . . . . . . . . . . . . . . . . 99
6.1.1. Comparación del resultado numérico y experimental de las velocidades 100
6.2. Calibración del módulo morfodinamico . . . . . . . . . . . . . . . . . . . . . 112

&RQFOXVLRQHV
5HFRPHQGDFLRQHV
Índice de figuras

1.1. Río Madre de Dios ubicado en la zona de estudio tramo ’La Pastora - Puerto
Maldonado’. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2.1. Componentes de esfuerzo viscoso que actúan en un volumen de control. . . . 17

2.2. Entrada y salida de masa de flujo en un volumen de control. . . . . . . . . . 18
2.3. Fuerzas superficiales que actúan en una partícula de fluido. . . . . . . . . . . 19
2.4. Fuerzas en la partícula de sedimento con lecho inclinado. . . . . . . . . . . . 29
2.5. Fuerzas en la partícula de sedimento . . . . . . . . . . . . . . . . . . . . . . 32
2.6. Fuerzas que actuan en la partícula de sedimento por Yang(1996) . . . . . . . 33
2.7. Parámetro de Shields con respecto al número de Reynolds de la partícula. . . 35
2.8. Variación de las dimensiones de la velocidad critica. . . . . . . . . . . . . . . 36

3.1. Modelo hidráulico del río Madre de Dios - Tramo la Pastora construido en el
Laboratorio Nacional de Hidráulica. . . . . . . . . . . . . . . . . . . . . . . . 44
3.2. Prototipo del río Madre de Dios - Tramo La Pastora. . . . . . . . . . . . . . 44
3.3. Curva granulométrica del modelo hidráulico del rio Madre de Dios. . . . . . 49
3.4. Configuración del lecho después de realizar el experimento, se observa el cambio
de forma del lecho tipo duna. . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.5. Batimetría del lecho después de la etapa experimental, capturado por un Es-
caner Laser, la flecha indica la dirección del flujo . . . . . . . . . . . . . . . . 51
3.6. Sección transversal S2, Ver Figura 3.5. . . . . . . . . . . . . . . . . . . . . . 52
3.7. Sección transversal S4, Ver Figura 3.5. . . . . . . . . . . . . . . . . . . . . . 52
3.8. Sección transversal S6, Ver Figura 3.5. . . . . . . . . . . . . . . . . . . . . . 53
 ÍNDICE DE FIGURAS xiii

3.9. Sección transversal S8,Ver Figura 3.5. . . . . . . . . . . . . . . . . . . . . . . 53

3.10. Siete secciones transversales en planta y distribución de las velocidades en seis
puntos medidos desde la margen derecha hacia la margen izquierda. . . . . . 54
3.11. Las velocidades se midieron con Velocímetro Acústico Doppler en las siete
secciones transversales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.12. Velocímetro Acústico Biestático está colocado a una distancia de 5cm por
encima del lecho. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.13. Registro y almacenamiento de las velocidades en una laptop para su posterior
procesamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.14. Perfil de velocidad de la Sec1 en el plano x-y. . . . . . . . . . . . . . . . . . . 58
3.15. Perfil de velocidad de la Sec2 en el plano x-y,ver Figura 3.10. . . . . . . . . . 58
3.16. Perfil de velocidad de la Sec3 en el plano x-y,ver Figura 3.10. . . . . . . . . . 59
3.17. Perfil de velocidad de la Sec4 en el plano x-y,ver Figura 3.10. . . . . . . . . . 59
3.18. Perfil de velocidad de la Sec5 en el plano x-y,ver Figura 3.10. . . . . . . . . . 60
3.19. Perfil de velocidad de la Sec6 en el plano x-y, Figura 3.10. . . . . . . . . . . 60
3.20. Perfil de velocidad de la Sec7 en el plano x-y, ver Figura 3.10. . . . . . . . . 61

4.1. Módulos del software Open Telemac Mascaret. . . . . . . . . . . . . . . . . . 63

4.2. Pendiente del lecho en el transporte del material de fondo y suspensión. . . . 71
4.3. Muestra el flujo secundario helicoidal generado en la vertical debido a la fuerza
centrífuga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.4. Representación de un elemento triangular 2D y número de nodos. . . . . . . 78

5.1. Dimensiones del domino físico. . . . . . . . . . . . . . . . . . . . . . . . . . . 85

5.2. Un total de 1239922 puntos vizualizado en el programa Klue Kenue. . . . . . 86
5.3. Se realizó el resample (mismo distanciamiento entre nodos) del dominio . . . 86
5.4. Seleccionar la opción de T3 Mesh Generator . . . . . . . . . . . . . . . . . . 87
5.5. Se selecciono la distancia entre nodos que es de 0.3 m . . . . . . . . . . . . . 87
5.6. Opción interpolacion (2D Interpolator) . . . . . . . . . . . . . . . . . . . . . 88
5.7. Asignarle elevación a la malla generada por medio de (Map Object) . . . . . 88
5.8. Archivo SELAFIN Object. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
 ÍNDICE DE FIGURAS xiv

5.9. Archivo Boundary Conditions (Conlim) aguas arriba. . . . . . . . . . . . . . 90

5.10. Archivo Boundary Conditions (Conlim) aguas abajo. . . . . . . . . . . . . . 91
5.11. Condiciones de frontera (Boundary Conditions) inflow Q [m3 /s] y outflow H
[m] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
5.12. Descripción del archivo .cas de TELEMAC2D. . . . . . . . . . . . . . . . . . 94
5.13. Uso de los coeficientes de flujo secundario y esquema numérico . . . . . . . . 95
5.14. Archivo de extensión .cas para el módulo morfodinamico. . . . . . . . . . . . 96
5.15. Características físicas de la partícula de sedimento. . . . . . . . . . . . . . . 97
5.16. Acoplamiento del módulo hidrodinámico (TELEMAC-2D) y transporte de se-
dimentos (SISYPHE). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
5.17. Algoritmo de solución para el módulo de transporte de sedimentos . . . . . . 98

6.1. Perfil de velocidad en el plano x-y, se muestra que al considerar la corrección

de flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . 101
6.2. Coeficiente de correlación entre los datos obtenidos de la simulación numérica
y datos experimentales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
6.3. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 102
6.4. Coeficiente de correlación entre los datos obtenidos de la simulación numérica
y datos experimentales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
6.5. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 104
6.6. Coeficiente de correlación entre los datos obtenidos de la simulacion numérica
y datos experimentales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
6.7. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 105
6.8. Coeficiente de correlación entre los datos obtenidos de la simulación numérica
y datos experimentales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
6.9. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 107
 ÍNDICE DE FIGURAS xv

6.10. Coeficiente de correlación entre los datos obtenidos de la simulación numérica

y datos experimentales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
6.11. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 108
6.12. Coeficiente de correlación entre los datos obtenidos de la simulación numérica
y datos experimentales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
6.13. Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de
flujo secundario se aproxima a los datos experimentales. . . . . . . . . . . . . 110
6.14. Coeficiente de correlación entre los datos obtenidos de la simulación numérica
y datos experimentales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
6.15. Se muestra las velocidades máximas en el tramo curvo, donde la velocidad está
en el rango de 0-0.5m/s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
6.16. Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S2) . . . . . . . . 113
6.17. Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S4) . . . . . . . . 113
6.18. Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S6) . . . . . . . . 114
6.19. Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Seccción S8) . . . . . . . . 114
Índice de tablas

1.1. Estadística del cambio de canal a lo largo de cuatro ríos en la Amazonía

peruana antes de 1992 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

3.1. Relación de escalas Modelo - Prototipo . . . . . . . . . . . . . . . . . . . . . 45

3.2. Porcentaje que pasa o es retenido por los tamices. . . . . . . . . . . . . . . . 48
3.3. Velocidades medias en seis puntos de cada sección transversal . . . . . . . . 57
3.4. Tirantes en seis puntos en las secciones 1,2 hasta 7 medidos desde la margen
izquierda hacia la margen derecha . . . . . . . . . . . . . . . . . . . . . . . . 61

6.1. Datos de la velocidad en la sección 1 [sec1] . . . . . . . . . . . . . . . . . . . 100

6.2. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
6.3. Datos de la velocidad en la sección 2 [sec2] . . . . . . . . . . . . . . . . . . . 102
6.4. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
6.5. Datos de la velocidad en la sección 3 [sec3] . . . . . . . . . . . . . . . . . . . 103
6.6. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
6.7. Datos de la velocidad en la sección 4 [sec4] . . . . . . . . . . . . . . . . . . . 105
6.8. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
6.9. Datos de la velocidad en la sección 5 [sec5] . . . . . . . . . . . . . . . . . . . 106
6.10. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
6.11. Datos de la velocidad en la sección 6 [sec6] . . . . . . . . . . . . . . . . . . . 108
6.12. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
6.13. Datos de la velocidad en la sección 7 [sec7] . . . . . . . . . . . . . . . . . . . 109
6.14. Parámetros Estadísticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
Capítulo ,

Introducción

Resumen:
En este capítulo se describe el problema que afecta a ríos meándricos y cómo la erosión
lateral afecta al río Madre de Dios en tramo la Pastora. Debido a la magnitud del problema
se construye un modelo hidráulico en el Laboratorio Nacional de Hidráulica - UNI.

1.1. Descripción del problema

Los ríos de la selva peruana presentan cauces con pendientes bajas, las aguas discurren en
un suelo aluvial y los caudales varían dependiendo de la época del año, y sus cursos son:
anastomosados, trenzados y meándricos, constituidos por curvas alternadas en forma de S.
Sioli et al.(1984), señalan que los ríos amazónicos a menudo se dividen en tres categorías:
en función a sus propiedades físicas, composición química de sus aguas y en la práctica de
acuerdo con el transporte de sedimentos en suspensión y fondo.
La tasa de erosión lateral (migración) es más rápido en los ríos de la Selva Baja. Una medida
de cómo los ríos han migrado es observando y comparando imágenes satelitales en el tiempo.
Por ejemplo Puhakka et al. (1992), utilizaron imágenes satélites Landsat MSS para observar
la migración lateral de algunos ríos de la selva peruana. Ellos observaron las estadísticas
de cambio del canal a lo largo de cuatro ríos (Amazonas, Marañón, Ucayali y Manu) en la
amazónica peruana, en un período de 4 años. Las observaciones señalan que el río Amazonas,
de tipo anastomosado, presenta tasa promedio de erosión lateral de 115 m/año, con un
 CAPÍTULO 1. INTRODUCCIÓN 2

máximo de erosión de 300 m/año y un total de área erosionada de 803.3 ha/año. El río
Ucayali presenta tasas máximos de erosión de 150 m/año ; y el río Manu, que tiene canales
meándrico, presenta una tasa promedio de erosión lateral de 12 m/año. La tasa máxima de
erosión del río Manu fue de 70 m/año y un área total erosionada de 0.93 ha/año ver Tabla
(1.1).

Tabla 1.1: Estadística del cambio de canal a lo largo de cuatro ríos en la Amazonía peruana
antes de 1992
.
Descripción Amazonas Marañon Ucayali Manu
Tipo de canal anastomosado anastomosado meándrico meándrico
Longitud de faja 110 90 110 70
meándrica (km)
Sinuosidad 1.5 1.2 1.9 2.1
Ancho del canal 4 4 4 14
min-max.(km)
Tasa de erosión lateral 115 110 78 12
promedio(m/año)
Tasa máxima de 300 200 150 70
erosión lateral(m/año)
Mayor unidad 1.5 0.4 0.5 0.03
erosinada (km2 /año)
Área total de 803.3 268.2 261.7 0.93
Erosión (ha/año)
Número de unidades 6.5 4 5.2 13.1
de bosque erosionada
Fuente:Estudio de Línea de Base – FINCyT.

La erosión lateral es un fenómeno natural que se presenta en casi todos los ríos tropicales;
sin embargo, en los últimos años estos procesos se han acelerado por causas antropogénicas,
produciendo destrucción de las riveras, cauces y planicies de inundación. Una causa ha sido
 CAPÍTULO 1. INTRODUCCIÓN 3

la “fiebre de oro” que desde la década de los 70 ha llevado a centenares de miles de mineros a
buscar oro en los ríos provocando serios impactos en los ecosistemas fluviales más complejos
y diversos de la selva peruana. La minería aurífera en suelos aluviales de la Amazonía se está
expandiendo rápidamente en los últimos años más que en cualquier época historica.(Indeci,
2006).

Figura 1.1: Río Madre de Dios ubicado en la zona de estudio tramo ’La Pastora - Puerto
Maldonado’.

Fuente: Google Earth(2018), elaboración propia.

En el cauce principal del río Madre de Dios hay una alta carga de sedimentos en suspensión
y la tasa de erosión lateral del canal es aproximadamente 6 m/año generando fenómenos de
socavación lateral en la margen derecha del río. En el tramo curvo cerca de la ciudad de
Puerto Maldonado, que tiene una longitud de más de 1540 m, se está generando el colapso de
la ribera dejando desprotegido el talud. La erosión lateral afectará a los inmuebles, sistema de
 CAPÍTULO 1. INTRODUCCIÓN 4

abastecimiento de agua y desagüe de la ciudad de Puerto Maldonado, y en el futuro pondrá

en riesgo más de 1 km de la carretera interoceánica. Así mismo ha dejado a los bancos con
taludes muy empinados, la verticalidad en los taludes hace que los bancos sean más inestables
ante el paso de cualquier incremento y disminución del caudal, generando la caída de grandes
masas de suelo debido a que el material es muy fácil de desprenderse, ya que la mayoría de
estos bancos presentan suelos finos como: arcillas, limos y arenas. Los bancos erosionados
presentan una pendiente inclinada,que en algunos tramos van desde los 560 hasta 900 .
Para proteger las riberas de los ríos existen estructuras de estabilización que mitigan la erosión
lateral de los canales aluviales. Los métodos de estabilización de ríos pueden ser clasificados
bajo dos enfoques:

Fortalecimiento de la ribera.

Reducción de las fuerzas.

De acuerdo a Julien (2002) ambos enfoques son necesarios para mitigar los efectos de la ero-
sión lateral en los ríos Amazónicos ya que sus riberas están compuestas de materiales que son
muy fácil de erosionarse, y como la geometría en planta del río es curvo, genera flujos alta-
mente tridimensionales que ayudan aún más a erosionar los bancos débiles. Es por esto que,
si se reduce las fuerzas hidrodinámicas, la erosión lateral podría ser mitigada. Las estructuras
de control de flujo son diseñadas para reducir las fuerzas hidrodinámicas que impactan en las
riberas mediante el control de la dirección del flujo, la velocidad o del tirante del agua que
fluye.
Conocer la morfología de un río meandrico es fundamental para plantear soluciones técnicas
adecuadas en el diseño y construcción de obras fluviales. El comportamiento morfológico se
puede estudiar en un modelo hidráulico, o implementando un modelo numérico hidrodiná-
mico y de transporte de sedimentos. Estos modelos, que casi siempre son complementarios,
ayudan a entender, reproducir y predecir los procesos físicos que ocurren en un río, y así
constituyen herramientas muy valiosas para el diseño hidráulico y estructural de las obras
fluviales.
El modelamiento hidráulico o numérico de transporte de sedimento, de por sí es difícil aún
para canales artificiales prismáticos y, con mayor razón, lo es para ríos. Si un rio presenta
 CAPÍTULO 1. INTRODUCCIÓN 5

curvatura, el modelamiento en el tramo curvo presenta mayor dificultad lo cual no es posible

predecir con un modelo numérico bidimensional. Un río meandrico por naturaleza presenta
tramos con curvatura. En consecuencia, el modelamiento de transporte de sedimento en tra-
mos curvos es un problema de mayor complejidad. En general, el transporte se caracteriza
por un alto nivel de sedimentación en el banco interior y una marcada erosión en el banco
exterior. Esta dinámica se traduce finalmente en un movimiento horizontal, o migración la-
teral, del cauce del río.
El problema de la migración de los ríos, se debe a la presencia de flujos secundarios en la
sección del tramo curvo. Las corrientes secundarias superficiales, y las que se encuentran
cerca del lecho, generan el fenómeno de erosión en un banco y sedimentación al otro lado, un
fenómeno de naturaleza tridimensional donde se presentan flujos helicoidales que empiezan
cuando el agua entra en el tramo con curvatura del cauce.
El estudio en el modelo hidráulico llevado a cabo en el Laboratorio Nacional de Hidráulica de
la Universidad Nacional de Ingeniería, para condiciones sin estructuras de protección ribereña
del tramo del río Madre de Dios se pudo observar que los efectos y la fuerza del flujo secunda-
rio varía considerablemente. Esto también se puede apreciar con una simple relación del radio
de curvatura local y el ancho del canal, haciendo que este tramo del río presente una alta
curvatura que induce flujo helicoidal o tridimensional con altas velocidades,alta turbulencia
y altos esfuerzos cortantes que producen erosión en el lecho del río; por otro lado, se hace
aún más complicado en los casos donde la batimetría o geometría del río es muy compleja.
El modelo hidráulico del tramo del río Madre de Dios permitió comprender los fenómenos
presentes en este río, pero es necesario utilizar un modelo numérico para complementar la
investigación, ya que los equipos empleados miden de forma puntual los datos. Un modelo
numérico puede explorar mejor la física del problema además su costo es menor comparado
con un modelo físico y las simulaciones numéricas de alta resolución podrían ayudar entender
dichos fenómenos.
La turbulencia en el tramo de estudio es usualmente bastante fuerte debido a los límites irre-
gulares, aunque no podamos observar visualmente la turbulencia interna que está generando
altos esfuerzos de corte cerca del lecho y en la superficie libre. Por otro lado, la relación del
transporte de sedimentos de fondo con el transporte de sedimentos en suspensión es de 1/10,
 CAPÍTULO 1. INTRODUCCIÓN 6

este valor podría deberse a una alta turbulencia que mantiene a las partículas en suspensión.
Visualizar y cuantificar la turbulencia es casi imposible ya que es un fenómeno muy comple-
jo, que está acompañado de vórtices tridimensionales que en su mayoría descienden hacia el
fondo del río y que varían en el tiempo.
En general, en cuanto a estudios en proyectos hidráulicos, un estudio numérico es una herra-
mienta que nos puede ayudar a complementar los estudios experimentales llevados a cabo en
laboratorio. Los modelos hidráulicos son muy caros, demandan mucho tiempo en la construc-
ción, calibración y toma de datos, por otro lado, los estudios de campo tampoco son faciles
de llevar acabo, sobre todo no se puede tomar datos debido a la complejidad, variación en el
tiempo y los diversos parámetros que resultan difícil de medir. En este contexto la simulación
numérica es una herramienta que presenta ventajas.

Meanderización de ríos
Los ríos meándricos se caracterizan por una sucesión de bucles curvos alternados, cada
bucle es un tramo del río entre dos puntos de inflexión; y un meandro consiste en un
par de bucles en dirección opuesta.
Numerosas hipótesis se han planteado para explicar la razón de la meanderización,
siendo éste un fenómeno en el cual la vorticidad juega un papel muy importante. Al
observar el fenómeno desde el punto de vista de la circulación de flujos secundarios en el
tramo curvo, la meanderización es generada a través de un cambio en la orientación del
río aguas abajo. Aquí las líneas de corriente cerca de la superficie son deflectadas hacia
el banco exterior, mientras las que están cerca del lecho son desviadas al banco interno.
La velocidad cerca al lecho, el esfuerzo cortante tangencial del lecho y el arrastre de las
partículas es generalmente direccionada hacia el banco interno produciéndose erosión
cerca al banco externo y sedimentación en el banco interno.
La migración de un río meándrico es el resultado de la erosión del banco externo com-
binado con la sedimentación cerca del banco interno.La erosión en el pie del banco
externo desplaza el thalweg hacia fuera del río curvo y causa empinamiento del banco
externo. El incremento del empinamiento del material del banco externo es una de las
causales de falla.
Las fallas en los bancos pueden ser caracterizados por la forma de la superficie de falla
 CAPÍTULO 1. INTRODUCCIÓN 7

o por el modo de falla.Por la superficie de la falla puede ser rotacional o plana.

Las fallas rotacionales son perjudiciales por que ocasionan pérdidas de tierra y por lo
general son menos comunes a lo largo de las riberas (Simon, 1989); en este tipo de falla
la tensión de cizallamiento aumenta más rápidamente con la profundidad que realiza
la resistencia al corte (Terzaghi et al., 1996).
Las fallas planas generalmente resultan de la erosión en la parte inferior de la ribera o el
desarrollo de una grieta de tensión casi vertical en la parte superior del banco (Lohnes,
1968).

Modos de falla típicos en los ríos aluviales

Las fallas en los bancos se han reportado con mayor frecuencia durante el periodo de
crecida. Existen tres modos de falla típicos en los ríos aluviales:
a) Con material granular no cohesivo: la pérdida de granos al pie del banco exterior
induce el deslizamiento del material granular, tan pronto al ángulo de inclinación excede
el ángulo de reposo del material.
b) Con material cohesivo: la falla rotacional es típica y la presencia de la tensión de
corte podría acelerar los procesos de erosión del banco.
c) Con depósitos estratificados: el material no cohesivo suprayacente es movilizado,
dejando así el material cohesivo suprayacente sujeto a la tensión de grietas y fallas en
voladizo. En general, los bancos más erosivos son los que están compuesto de arena y
limo, mientras que los menos erosivos de arcilla y grava.

1.2. Formulación del problema

1.2.1. Problema general

¿Es posible simular numéricamente la dinámica del flujo secundario helicoidal y del transporte
de sedimentos en el tramo curvo de un cauce de lecho movil, tal como elmodelo hidráulico
del rio Madre de Dios?
 CAPÍTULO 1. INTRODUCCIÓN 8

1.2.2. Problemas específicos

¿Es posible simular numéricamente las características del flujo helicoidal que se forma
en la sección transversal del tramo curvo del modelo hidráulico del rio Madre de Dios?

¿De qué manera afecta el flujo helicoidal a las riberas del tramo curvo del cauce del
modelo hidráulico del rio Madre de Dios?

¿De qué manera afecta el flujo helicoidal a la dinámica del transporte de sedimentos en
el tramo curvo del modelo hidráulico del rio Madre de Dios?

¿Es posible simular numéricamente la configuración del proceso de erosión-sedimentación

en la sección transversal del tramo curvo del modelo hidráulico del rio Madre de Dios?

1.3. Objetivos

1.3.1. Objetivo general

Realizar una simulación numérica de la dinámica del flujo secundario helicoidal y del trans-
porte de sedimentos en un tramo curvo del modelo hidráulico del Río Madre de Dios, para
analizar la influencia del flujo secundario en la morfología del lecho.

1.3.2. Objetivos específicos

Simular numéricamente las características del flujo helicoidal que se forma en la sección
transversal del tramo curvo del modelo hidráulico del rio Madre de Dios.

Interpretar los efectos del flujo helicoidal sobre las riberas del tramo curvo del modelo
hidráulico del rio Madre de Dios.

Interpretar los efectos del flujo helicoidal en la dinámica del transporte de sedimentos
en un tramo curvo del modelo hidráulico del rio Madre de Dios.

Simular numéricamente la configuración del proceso de erosión y sedimentación en la

sección transversal del tramo curvo del modelo hidráulico del rio Madre de Dios.
 CAPÍTULO 1. INTRODUCCIÓN 9

1.4. Justificación
Debido al acelerado proceso de erosión lateral en ríos meándricos es necesario comprender la
dinámica del flujo y el transporte de sedimentos en este tipo de ríos, ya que en muchos de ellos
se proyectan estructuras hidráulicas los cuales se ven afectadas por el proceso de transporte
de sedimentos. En el caso específico del rio Madre de Dios, la migración lateral afecta a
los inmuebles aledaños, al sistema de abastecimiento de agua y desagüe y a la carretera
interoceánica que pasa por la ciudad de Puerto Maldonado.
Por otro lado, la simulación numérica constituye una ayuda para analizar el flujo helicoidal y el
proceso de transporte de sedimentos cerca de las riberas empinadas de ríos con tramos curvos.
Por lo tanto, con esta herramienta numérica se puede simular el fenómeno de migración lateral
que se presenta en ríos meándricos, estimar su magnitud y, de ser el caso, analizar alternativas
de mitigación.
La principal contribución de este trabajo es mostrar que en un río meándrico, se presenta
flujo secundario tridimensional. En situaciones de este tipo es importante usar la corrección
de flujo secundario en modelos numéricos bidimensionales con el objetivo de parametrizar las
variables en la vertical o usar un modelo numérico tridimensional.

1.5. Metodología
La simulación numérica de flujo secundario helicoidal y transporte de sedimentos, se ha
realizado resolviendo las ecuaciones de conservación de momentum lineal y de continuidad,
complementada con fórmulas empíricas y semi-empíricas del transporte de sedimentos en
ríos. Estas ecuaciones están implementadas en el software libre de código abierto denominado
Open TELEMAC MASCARET. Los resultados de la simulación numérica se ha constrastado
con la información obtenida en el modelo hidráulico del río Madre de Dios realizado en el
Laboratorio Nacional de Hidráulica de la Universidad Nacional de Ingeniería.
Durante el desarrollo de la tesis se ha seguido los siguientes pasos:
 CAPÍTULO 1. INTRODUCCIÓN 10

1.5.1. Acopio de Información

Revisión de bibliografía general sobre el tema, con énfasis sobre la dinámica de fluidos compu-
tacional. La revisión bibliográfica abarcó trabajos de tesis, artículos de investigación, mono-
grafías y otros estudios realizados sobre: hidrodinámica de canales curvos, transporte de
sedimentos en canales curvos, desarrollo y evolución de lechos móviles. Asimismo, abarcó
búsqueda del tema en páginas de internet y acopio de información relacionado con el softwa-
re de simulación, etc.
Acopio y sistematización de datos obtenidos en los ensayos del modelo hidráulico del río Ma-
dre de Dios, de escala lineal de 1:60, realizado en en el Laboratorio Nacional de Hidráulica
(LNH-UNI), donde se recolectaron datos de velocidad, tirante y se escaneó la topografía del
modelo hidráulico antes y después del ensayo experimental, estos datos permitieron validar
el modelo numérico usado.

1.5.2. Manejo del software Open TELEMAC

Se recopiló el manual de usuario y “demos” del software Open TELEMAC MASCARET. Se

analizó el método de discretización utilizado, tanto del dominio (enmallado), como de las
ecuaciones involucradas. Se identificó las ecuaciones discretizadas en el software, tanto del
módulo hidrodinámico como del módulo de transporte de sedimentos. Se realizó simulaciones
con los módulos TELEMAC-2D y SISYPHE para casos simples o hipotéticos, con la finalidad
de obtener pericia en el manejo de los módulos.
Se adecuó, las ecuaciones del software (Open TELEMAC MASCARET) para el estudio
del problema objetivo. Estos pasos comprendieron tanto al módulo hidrodinámico como al
módulo de transporte de sedimentos.

1.5.3. Simulación numérica computacional

Este paso consistió en:

Delimitar el dominio computacional en programa Blue Kenue.

Generación de la malla computacional no estructurada, con un espaciamiento de 0.3 m.

 CAPÍTULO 1. INTRODUCCIÓN 11

Asignar las condiciones de frontera aguas arriba (caudal) y aguas abajo (tirante),además
generar los archivos de extensión .cli y .slf, que son los archivos necesarios para la
simulación.

Descargar el archivo de extensión .cas de la pagina web www.opentelemac.org, que

tiene todos los parametros hidráulicos como: (Rugosidad, time de paso, etc).

Se realizar la simulación hidrodinámica con TELEMAC-2D para generar el archivo

hort-start.slf que contiene resultados de las variables hidrodinámicas,necesarios para
inicializar las variables y hacer la simulación de forma acoplada.

Simular, de forma acoplada con el modulo TELEMAC-2D y SISYPHE empleando los

archivos de condición de frontera (.cli), malla computacional (.slf), el archivo .cas y el
hort-start generado previamente. Ademas, realizar analisis de sensibilidad para diferen-
tes pasos de tiempos (time-steps) hasta obtener resultados coherentes.

Con los resultados obtenidos hacer el post-proceso para contrastar con los datos expe-
rimentales obtenidos del modelo hidráulico.(Ver capitulo 6)

1.5.4. Validación del modelo numérico computacional

Se realizó la validación del modelo numérico con los resultados experimentales obtenidos en
el modelo hidráulico del río Madre de Dios tramo la Pastora, (ver item 1.5.1).

1.6. Motivación
El interés por realizar el presente trabajo surgió cuando tuve la oportunidad de realizar
prácticas pre-profesionales en el Laboratorio Nacional de Hidráulica que se encuentra en la
Universidad Nacional de Ingeniería (UNI), donde participe durante los ensayos y en proce-
samiento de datos del modelo hidráulico del rio Madre de Dios.
Por otro lado, en el Laboratorio Nacional de Hidráulica se presentó la oportunidad de usar
el software Open TELEMAC MASCARET para resolver las ecuaciones de flujo de fluidos y
las ecuaciones empíricas del transporte de sedimentos.
Capítulo ,,

Marco teórico

Resumen:
En este capítulo se describe las ecuaciones que gobiernan el movimiento del fluido, transporte
de sedimento y sus propiedades.

2.1. Antecedentes
La simulación numérica de la hidrodinámica, transporte de sedimentos y la morfología de
lecho de ríos son temas relativamente reciente y va emparejado con el avance de la tecnología
de la computación y el desarrollo de métodos numéricos, sin estas herramientas sería casi
imposible abordar el tema.
Las ecuaciones que gobiernan a los citados temas son ecuaciones tridimensionales, no perma-
nentes y en derivadas parciales, por lo que su solución es vía la simulación numérica.
Existen varios software de simulación hidrodinámica y transporte de sedimentos, los más co-
nocidos en el mundo académico e industrial, que son libres y validados, son: open TELEMAC
MASCARET, y DELFT-3D (versión estudiantil).
A la fecha, Open TELEMAC MASCARET tiene la última version v8p0 y cuenta con la
version 2D, 3D y el modulo de transporte de sedimentos.
Open TELEMAC fue desarrollado por el Laboratorie Nacional de Hidráulica y Medio Am-
biente del departamento de la División de Investigación y Desarrollo de Electricidad de Fran-
cia, (Laboratoire National d’Hydraulique et Environnement - LNHE, Recherche et Dévelop-
 CAPÍTULO 2. MARCO TEÓRICO 13

pement de Électricité de - EDF - R & D) Francia.

El código TELEMAC-2D resuelve ecuaciones de flujo de superficie libre promediadas en la
profundidad (ecuaciones de Saint Venant). La principal aplicación de TELEMAC- 2D está
en hidráulica marítima e hidráulica fluvial y tiene en cuenta los siguientes fenómenos:

Turbulencia.

Flujo subcrítico y supercrítico.

Influencia de la temperatura horizontal y gradiente de la salinidad en la densidad.

Planicies de marea y llanura de inundación.

Tratamiento de singularidades: vertederos, diques y alcantarillas.

Transporte de sedimentos.

A continuación, se presenta los antecedentes internacionales más próximos en el estudio de

la dinámica del flujo secundario helicoidal y del transporte de sedimentos en el tramo curvo
de ríos, especialmente los que han empleado el software Open TELEMAC MASCARET:

Villaret et al. (2013), en un artículo titulado “Morphodynamic modeling using the

Telemac finite-element system” describen que la principal originalidad del sistema
de código abierto open TELEMAC MASCARET radica en la eficiencia y la flexibilidad
de los elementos finitos, además de la optimización de esquemas numéricos, con el
progreso en el rendimiento de los ordenadores.

Morell et al. (2014) elaboraron un estudio bajo el título “Flow pattern at a river
diffluence at the alluvial system of the Paraná River” , en la ciudad de Santa Fe,
(476 km al noroeste de Buenos Aires). El objetivo del estudio fue mostrar el patrón de
flujo desarrollado en la difluencia y plantearon una solución al problema de sedimen-
tación crónica visto en la entrada del canal de acceso. Para dicho estudio utilizaron el
software Open TELEMAC MASCARET.

Tassi (2014), en un artículo titulado "Flow and bed morphology in large-amplitude

meanders, comparison of two and three-dimensional numerical models", pre-
sentó la simulación bidimensional y tridimensional del flujo en canales meandricos de
 CAPÍTULO 2. MARCO TEÓRICO 14

gran amplitud empleando la metodología RANS (Reynolds-Averaged Navier Stokes),

basado en los estudios de Whiting y Dietrich. En la simulación demostraron que con el
modelo bidimensional (profundidad promediada) no es posible tener en cuenta efectos
tridimensionales.

Riesterer et al. (2016) presentaron un artículo títulado “Bed load transport mode-
ling of a secondary flow influenced curved channel with 2D and 3D numerical
models”. Los autores realizaron simulación numérica de un canal meándrico de labo-
ratorio, donde se analizaron la influencia del flujo secundario en el desarrollo del lecho
y en el transporte de carga de fondo del canal. Además, realizaron la comparación de
la resolución espacial de la simulación hidrodinámica, llevado a cabo con software open
TELEMAC MASCARET.

En el ámbito nacional no se ha encontrado trabajos relacionados al tema.

2.2. Metódos para describir el flujo de fluidos

Exites dos diferentes enfoques de descripción del movimiento de los fluidos: Euleriano y Lan-
gragiano.
El enfoque Euleriano analiza el movimiento del fluido mediante un elemento fluido infinite-
simal o volumen de control que es fijo en el espacio. Este es el enfoque más común ya que
el interés generalmente está en la variación de las propiedades físicas (por ejemplo, presión,
densidad, temperatura y velocidad) en un punto fijo en el espacio en lugar de los puntos en mo-
vimiento. En este enfoque las ecuaciones gobernantes se obtienen directamente del elemento
fluido infinitesimal o volumen de control mediante la aplicación de principios fundamentales
(por ejemplo: continuidad, Segunda Ley de Newton, Primera Ley de la Termodinámica) y
se expresa en ecuaciones diferenciales parciales. El enfoque Euleriano da lugar a la forma
conservativa de las ecuaciones gobernantes del fluido (Versteeg, 1995).
El enfoque Langragiano rastrea las propiedades físicas (por ejemplo: densidad, presión, tem-
peratura y velocidad) de una partícula fluida mientras se mueve a través del elemento de
fluido. En este enfoque las ecuaciones de flujo de fluidos son directamente obtenidos desde el
 CAPÍTULO 2. MARCO TEÓRICO 15

elemento fluido infinitesimal aplicando principios fundamentales de ecuaciones diferenciales

parciales.
El enfoque Langragiano da lugar a la "forma no-conservativa"de las ecuaciones gobernantes
del fluido (Versteeg, 1995).

Derivada material
La derivada material se usa para describir el cambio de las propiedades de una partícula
fluida. Cualquier propiedad φ de una partícula por unidad de masa es usualmente
una función de su posición y del tiempo t φ(x, y, z, t). La tasa de cambio de esta
propiedad φ de una partícula fluida con respecto al tiempo ( Dφ
Dt
), es llamado derivada
material (White, 1991), derivada total o sustancial (Anderson & Wendt, 1995), derivada
absoluta (Chanson, 1999).( Dφ
Dt
) puede ser presentado matemáticamente por el enfoque
Langragiano utilizando la expansión en la serie de Taylor sin considerar los términos de
orden superior, o como en el presente caso, empleando la regla de la cadena del cálculo
diferencial.
Dφ dφ ∂φ ∂φ dx ∂φ dy ∂φ dz
= = + + + (2.1)
Dt dt ∂t ∂x dt ∂y dt ∂z dt
dx dy dz
Dado que una partícula de fluido sigue la línea de corriente, dt
= u, dt
=vy dt
= w,
con lo que resulta:
Dφ ∂φ ∂φ
= + u.∇φ = + u.gradφ (2.2)
Dt ∂t ∂t
∂φ
En la ecuación (2.2) el término ∂t
es llamado "derivada local". El cual físicamente es
la tasa de cambio de alguna propiedad φ con respecto al tiempo en un punto fijo o
elemento fluido. El termino u.gradφ es la derivada advectiva. El cual físicamente es la
tasa de cambio de alguna propiedad φ debido al movimiento del elemento fluido de un
punto a otro en el campo de flujo donde las propiedades son espacialmente diferentes.
Dφ
Es importante notar que Dt
indica la tasa de cambio de una propiedad φ por unidad
de masa. Por lo tanto debe ser multiplicado por la densidad del fluido ρ para obtener
la tasa de cambio de alguna propiedad φ por unidad de volumen. La tasa de cambio de
la propiedad φ por unidad de volumen de un elemento fluido está dado por la ecuación:
∂φ Dφ
ρ =ρ − u.∇(ρφ) (2.3)
∂t Dt
 CAPÍTULO 2. MARCO TEÓRICO 16

En otras palabras, la derivada material establece que la tasa de incremento de cual-

quier propiedad φ de un elemento fluido por unidad de volumen es igual a la tasa de
incremento de cualquier propiedad φ de una partícula fluida (moviéndose a través del
elemento) por unidad de volumen, menos la tasa neta de φ fuera del elemento fluido
por unidad de volumen.
Las ecuaciones de momento y energía se derivan usando el concepto de derivada mate-
rial.

2.3. Fuerzas en un elemento de fluido

Hay dos tipos diferentes de fuerzas que actúan en el volumen de control (Anderson & Wendt,
1995).

Fuerzas volumétricas: actúan directamente sobre la masa volumétrica del elemento

fluido. Estas fuerzas actúan a una distancia de la superficie, por ejemplo la fuerza
gravitacional, centrifuga, coriolis y fuerzas electromagnéticas.

Fuerzas superficiales: actúan directamente sobre la superficie del elemento fluido.

Estas fuerzas se deben a la distribución de la presión impuesta por el fluido exterior que
rodea el elemento fluido, que actúa a través de los esfuerzos normales y la distribución
del esfuerzo viscoso impuesto por el fluido externo en la superficie por medio de la
fricción.

Por simplicidad la contribución debido a las fuerzas volumétricas son frecuentemente agru-
pados como "termino fuente", mientras que las contribuciones debido a las fuerzas de super-
ficiales son considerados separadamente.
De acuerdo con (Anderson & Wendt, 1995), el esfuerzo del elemento fluido puede ser definido
en términos de presión y nueve componentes de esfuerzos viscosos (normal y tangencial) que
actúan en el elemento fluido.
 CAPÍTULO 2. MARCO TEÓRICO 17

Figura 2.1: Componentes de esfuerzo viscoso que actúan en un volumen de control.

Fuente: (Versteeg, 1995)

Las componentes del esfuerzo viscoso τij dependen del gradiente de velocidad (por ejemplo
∂u
∂x
). De acuerdo con (Anderson & Wendt, 1995), el esfuerzo viscoso normal (tal como: τxx )
es mucho más pequeño que el esfuerzo viscoso tangencial, por lo tanto estos esfuerzos son
insignificantes en el análisis del flujo viscoso.

2.4. Ecuaciones diferenciales de flujo de fluidos

2.4.1. Ecuación de conservación de masa

La ecuación de conservación de masa deriva del principio de conservación de masa. La tasa

de incremento de masa en un elemento fluido es igual a la tasa neta de flujo de masa en el
elemento de fluido.
Si consideramos un elemento fluido infinitesimal de lados dx,dy y dz fijo en el espacio como
se muestra en la Figura (2.2).
 CAPÍTULO 2. MARCO TEÓRICO 18

Figura 2.2: Entrada y salida de masa de flujo en un volumen de control.

Fuente: (Versteeg, 1995)

∂ρ
La tasa de incremento de masa en el elemento fluido por unidad de volumen es ∂t
. La tasa
neta de flujo de masa por unidad de volumen a través de sus caras puede ser expresado
como: ∂(ρu)
∂x
+ ∂(ρv)
∂y
+ ∂(ρw)
∂z
, donde: u,v, y w es la velocidad a lo largo de la dirección x,y,z
respectivamente.
En un fluido compresible, la ecuación de conservación de masa no estacionario puede ser
escrito como: (Anderson & Wendt, 1995).

∂ρ ∂(ρu) ∂(ρv) ∂(ρw)

+ + + =0 (2.4)
∂t ∂x ∂y ∂z
∂ρ
Para un fluido incompresible, de densidad constante ρ, ∂t
= 0 la ecuación (2.4) puede ser
reescrito como:
∂(u) ∂(v) ∂(w)
+ + =0 (2.5)
∂x ∂y ∂z

2.4.2. Ecuación de conservación de momentum

La ecuación de momentum se basa en el principio de conservación de momentum lineal, o en

la segunda ley de Newton. "La tasa de incremento de momentum en el volumen de control
es igual a la suma de fuerzas que actúan en la partícula fluida".
 CAPÍTULO 2. MARCO TEÓRICO 19

Usando el concepto de la derivada material, la tasa de incremento en las direcciones x,y, e

z por unidad de volumen de una partícula fluida de densidad constante está dado por: ρ Du
Dt
,
ρ Dv
Dt
y ρ Dw
Dt
respectivamente.
En la Figura (2.3) se ilustra solo el componente de las fuerzas superficiales (usando la serie
de Taylor sin tomar en cuenta los términos de orden superior) debido a la presión p y al
esfuerzo cortante τxx , τyx y τzx .
Un diagrama similar puede ser construido para las componentes y e z de las fuerzas superfi-
ciales.

Figura 2.3: Fuerzas superficiales que actúan en una partícula de fluido.

Fuente: (Versteeg, 1995)

La fuerza volumétrica se describe mediante el vector SM agrupando el termino fuente de una

partícula fluida por unidad de volumen y el tiempo.SM x ,SM y y SM z son las componentes de
SM en las direcciones x,y e z respectivamente.
La resultante de las fuerzas superficiales en la dirección x por unidad de volumen de una
partícula fluida puede ser escrito como (Versteeg, 1995):

∂(−p + τxx ) ∂τyx ∂τzx

+ + =0 (2.6)
∂x ∂y ∂z

El efecto general de las fuerzas volumétricas en la dirección x de la ecuación de momentum

 CAPÍTULO 2. MARCO TEÓRICO 

se puede incluir en la ecuación (2.7) agregando el termino fuente SM x .

El componente x de la ecuación de momentum establece que la tasa de cambio en la dirección
x de la partícula fluida es igual a la fuerza total en la dirección x del elemento debido a las
fuerzas superficiales más la tasa de aumento en la dirección x debido al término fuente
Versteeg (1995).
Du ∂(−p + τxx ) ∂τyx τzx
ρ = + + + SM x (2.7)
Dt ∂t ∂y ∂z
Similarmente en la dirección y e z la ecuación de momento puede ser derivado como:

Dv ∂τxy ∂(−p + τyy ) τzy

ρ = + + + SM y (2.8)
Dt ∂x ∂y ∂z

Dw ∂τxz ∂τyz ∂(−p + τzz )

ρ = + + + SM z (2.9)
Dt ∂x ∂y ∂z
El signo asociado con la presión (p) es opuesto al esfuerzo viscoso normal porque la conven-
ción de signos asume que el esfuerzo normal es positivo. Si la gravedad es la única fuerza
volumétrica, SM x = 0 y SM y = 0, mientras que SM z = ρg donde g es el vector gravedad para
un sistema de coordenadas con eje z concidiendo con la dirección vertical.
La ecuación de momento contiene incógnitas de esfuerzo viscoso τji , para la simplificación,
del cálculo del esfuerzo es usualmente basado en un fluido newtoniano, donde el esfuerzo
viscoso es proporcional a la tasa de deformación. Otros supuestos son la incompresibilidad
e isotropía de un fluido newtoniano,(Schlichting et al., 1955) muestran que la tasa lineal de
deformación de un elemento de fluido tiene nueve componentes en tres dimensiones, seis de
los cuales son independientes.
La ecuación de Stokes expresa los esfuerzos viscosos en términos de tasa de deformación y
proporcionalmente constante.Los nueve componentes de esfuerzo viscoso para flujo incompre-
sible puede ser expresado en términos de la viscosidad dinámica µ y la tasa de deformación
(Versteeg, 1995).
∂u ∂v ∂w
τxx = 2µ τyy = 2µ τzz = 2µ (2.10)
∂x ∂y ∂z

∂u ∂v
τxy = τyx = µ( + ) (2.11)
∂x ∂x
 CAPÍTULO 2. MARCO TEÓRICO 21

∂u ∂w
τxz = τzx = µ( + ) (2.12)
∂z ∂x

∂v ∂w
τyz = τzy = µ( + ) (2.13)
∂z ∂y
Sustituyendo las ecuaciones anteriores en la ecuación de momentum se tiene una ecuación
genérica del flujo de fluido que puede ser escrito usando la notación de Newton (White, 1991)
y una forma general (Versteeg, 1995):

Du ∂ ∂ui ∂uj
ρ = −∇p + [µ( + )] + SM (2.14)
Dt ∂xi ∂uj ∂ui

∂u 1 ∂p
+ div(uv) = − + vdivgrad(u) + SM x (2.15)
∂t ρ ∂x
∂v 1 ∂p
+ div(vu) = − + vdivgrad(v) + SM y (2.16)
∂t ρ ∂x
∂w 1 ∂p
+ div(wu) = − + vdivgrad(w) + SM z (2.17)
∂t ρ ∂x
Las ecuaciones anteriores son conocidas como la ecuación de Navier - Stokes para un fluido
incompresible y flujo newtoniano en tres dimensiones (i=1,2,3. j=1,2,3).

2.4.3. Ecuación de transporte

Se puede notar que hay similitudes significativas entre las diversas ecuaciones derivadas del
principio fundamental del flujo de fluidos, por ejemplo: conservación de masa, momentum
y la energía. Si φ es una variable general, la forma conservativa de toda la ecuación del
flujo de fluido, incluye ecuaciones para cantidades escalares tales como la temperatura y la
concentración de contaminantes, etc, se pueden escribir en forma siguiente (Versteeg, 1995):

∂ρφ
+ div(ρφu) = div(Γgradφ) + Sφ (2.18)
∂t

Dónde Γ es el coeficiente de difusión; en otras palabras, la ecuación (2.18) establece que la

suma de tasa de incremento de φ del elemento fluido más la tasa neta de flujo de φ a la salida
de un elemento fluido es igual a la suma de la tasa de incremento de φ debido a la difusión
más la tasa de incremento de φ debido al termino fuente.
 CAPÍTULO 2. MARCO TEÓRICO 22

La ecuación (2.18) también conocido como la ecuación de transporte (advección-diffusión)

para alguna propiedad φ. Destaca claramente los diversos procesos de transporte; la tasa de
cambio y el término advectivo, respectivamente, en el lado izquierdo y el término difusivo y
el término fuente,respectivamente, en el lado derecho.
En la solución de dinámica de fluidos computacional (CFD), el termino convectivo es casi
siempre de naturaleza hiperbólica, el cual es un problema sin disipación significativa. Por
otro lado, el termino difusivo es de naturaleza parabólica el cual es un problema con una
disipación significativa (Anderson & Wendt, 1995).

2.4.4. Métodos de solución numérica de las ecuaciones del flujo tur-

bulento

2.4.4.1. Simulación numérica directa (DNS)

Este método (cuyas iniciales provienen de la denominación inglesa Direct Numerical Simu-
lation) consiste, en realidad, en no utilizar ningún modelo para la turbulencia, sino realizar
discretización temporal y espacial que sean capaces de simular el flujo en un determinado
problema. La resolución directa de las ecuaciones de Navier-Stokes resulta hoy en día abor-
dable solo para un número muy limitado de problemas simples de interés académico. Los
grandes centros dedicados a la mecánica de fluidos disponen de líneas de investigación con
esta orientación, pero tanto las limitaciones en memoria de almacenamiento de las variables,
como el tiempo de cálculo hacen de momento impensable la solución generalizada de proble-
mas prácticos usando este tipo de técnica. Según Vandromme (1989), la primera solución de
este tipo se realizó en 1981 en la Universidad de Stanford.

2.4.4.2. Simulación de grandes remolino (LES)

Este tipo de técnica numérica reducen la complejidad de las ecuaciones que gobiernan con-
siderando sólo parte de los efectos de flujo turbulento. Se estudia el intercambio energético
entre las denominadas “fluctuaciones de gran escala” y se simula el efecto de las pequeñas
escalas de la turbulencia. Se trata de un tipo de modelo intermedio entre la simulación directa
y el promedia temporal de las ecuaciones de Navier-Stokes, que extiende el promedio tempo-
 CAPÍTULO 2. MARCO TEÓRICO 23

ral a la captura de ciertos efectos turbulentos básicos de forma numérica. En los modelos de
simulación de grandes remolinos, las ecuaciones no estacionarias del flujo se resuelven para el
flujo medio y para los grandes vórtices se modela el efecto de los pequeños vórtices. Aunque
sin llegar al extremo de la simulación directa, sólo es posible para problemas simplificados y
requiere capacidades de cálculo muy elevadas.

2.4.4.3. Metodologia (RANS)

Los métodologia promediada de las ecuaciones de Navier-Stokes (Reynolds Averaged Navier-

Stokes), llamado también ecuaciones de Reynolds, ha sido muy estudiada y resulta bastante
útil en la mayoría de los problemas prácticos. El procedimiento de promediar las leyes que
describen el movimiento de una partícula se introduce en las ecuaciones con el fin de obtener
el comportamiento promedio y la fluctuación de las distintas variables. Es decir se trata de
obtener una descomposición de las variables en su valor medio y su valor fluctuante. Por
ejemplo, para la velocidad, la descomposición sería:

→
− − −
→ →′
V = V +V (2.19)

dónde la componente promedio de la velocidad se obtiene haciendo la integral de la velocidad

→
−
instantánea V :
→
− ZT
1 →
−
V (t) = V (t)dt (2.20)
T
0

suponiéndose que el periodo de integración (T) es lo suficientemente grande en comparación

con la escala temporal de la turbulencia, pero lo suficientemente pequeño como para no cap-
tar cualquier fenómeno no estacionario distinto a la turbulencia.
La metodología RANS es bastante adecuado, pues la mayoría de los fenómenos no esta-
cionarios en mecánica de fluidos tiene lugar a frecuencias con rangos muy alejados de la
turbulencia.
El promedio temporal de las ecuaciones de Navier Stokes, da lugar a unos términos, denomi-
nados las tensiones de Reynolds (Reynolds stress), que involucra el producto de los términos
de fluctuación de la velocidad, cuya relación con las componentes promedios de flujo es desco-
nocida. Para obtener dicha relación es necesario introducir un modelo adicional, denominado
 CAPÍTULO 2. MARCO TEÓRICO 24

modelo de turbulencia o ecuaciones de cierre.

Las distintas posibilidades prácticas en cuanto a modelos de turbulencia son analizadas a
continuación. Habitualmente lo que interesa son los efectos de la turbulencia sobre valores
promedios de las variables: la velocidad promedio y la presión promedio en el caso del flujo
en un ducto.
Las ecuaciones de Reynolds considera la influencia de la fluctuación turbulenta en el flujo
→
−
medio con V (Con U,V y W son las componentes a lo largo de la dirección x, y e z).
En esta notación el vector fluctuación de velocidad es V ′ , con componente u’,v’,w’en las di-
recciones de los ejes x,y e z, respectivamente.
En el sistema de coordenadas cartesiana,la ecuación de continuidad para flujo medio es
(divV = 0) y las ecuaciones de Navier-Stokes para un flujo incompresible con viscosidad
constante son:

∂U 1 ∂P ∂u′2 ∂u′ v ′ ∂u′ w′

+ div(U U ) = − + νdivgradU + [− − − ] (2.21)
∂t ρ ∂x ∂x ∂y ∂z

∂V 1 ∂P ∂u′ v ′ ∂v ′2 ∂v ′ w′
+ div(V U ) = − + νdivgradV + [ − − ] (2.22)
∂t ρ ∂y ∂x ∂x ∂z

∂W 1 ∂P ∂u′ w′ ∂v ′ w′ ∂w′2
+ div(W U ) = − + νdivgradV + [ − − ] (2.23)
∂t ρ ∂z ∂x ∂y ∂z
Las ecuaciones (2.21),(2.22) y (2.23) muestra que se tienen más incognitas (10) que el número
de ecuaciones (4), situación conocida como el problema de cerradura. Para superar esta
situación la metodología RANS requiere el modelamiento de turbulencia o introducción de
modelos de turbulencia, para representar los términos de esfuerzo de Reynolds representados
por:

τxx = −ρu′2 τyy = −ρv ′2 τzz = −ρw′2 (2.24)

τxy = τyx = −ρu′ v ′ τxz = τzx = −ρu′ w′ τyz = τzy = −ρv ′ w′ (2.25)

Existen diferentes modelos de turbulencia, los modelos k-w y k-epsilon son los más conocidos.
 CAPÍTULO 2. MARCO TEÓRICO 25

Modelo de turbulencia k-ǫ estándar

Las dos ecuaciones obtenidas (uno para k y otro para (epsilon) no pueden ser resuel-
tas económicamente debido a su naturaleza de no cierre. Muchos investigadores como
(Launder & Spalding, 1983) han simplificado aún más las ecuaciones que gobierna k − ε
para producir técnicas de solución simple.Hrenya et al. (1998),(Iyer & Yavuzkurt, 1999),
informan estudios comparativos de este modelos de turbulencia κ − ε.
Entre estos modelos, el modelo estándar propuesto por (Launder & Spalding, 1983)
es popular por su robustez (Versteeg, 1995). El sistema open TELEMAC MASCA-
RET adopta el modelo original presentado por (Iyer & Yavuzkurt, 1999). De acuerdo
a TELEMAC-2D 2001, el modelo κ − ε puede ser expresado como:

Zz
1 1 ′ ′
k= u u · dz donde (i = 1 a 3) (2.26)
h 2 i j
zf

El cierre se hace por la hipótesis de Bousssinesq (Versteeg, 1995) que expresa el tensor
de Reynolds como una función de gradiente de velocidad media.

∂Ui ∂Uj 2 2
Rij = −u′i vj′ = νt ( + ) − κδij = 2νt En − κδij (2.27)
∂xj ∂xi 3 3

Donde i=1 a 2 y j=1 a 2. El último término de la ecuación (2.27) contiene k que es insig-
nificante (TELEMAC-2D,2001). Este término es una parte de la integral de la presión
en la ecuación de Navier Stokes que no se considera en TELEMAC por simplicidad. De
acuerdo a Kolmogorov (Pope, 2001), la viscosidad turbulenta es calculada como:

κ2
ν t = Cµ (2.28)
ε

En el software open TELEMAC MASCARET, la adaptación del modelo estándar κ − ε

es una extensión del modelo clásico presentado por (Rastogi & Rodi, 1978). Cabe
señalar que la integración vertical da como resultado un término de disipación debido
a la falta de uniformidad del perfil de velocidad, que no existe en el modelo clásico. El
sistema open TELEMAC MASCARET resuelve las siguientes ecuaciones:

∂κ 1 νt
+ ugrad(κ) = div(h ) + P − ε + Pkv (2.29)
∂t h σk gradκ
 CAPÍTULO 2. MARCO TEÓRICO 26

El primer término en el lado derecho es el término de difusión. Mientras que los términos
de producción se componen, por el gradiente de velocidad media horizontal y por otro
lado el término del esfuerzo vertical (TELEMAC-2D. 2001).

∂Ui ∂Uj ∂Ui

P = νt ( + ) donde i=1 a 2 y j=1 a 2 (2.30)
∂xj ∂xi ∂xj

u3∗ 1
Pkv = ck donde ck = √ (2.31)
h cf
p
u4∗ C2ε Cµ
Pεv = cε donde cε = 3,6 (2.32)
h2 cf3,4
Dónde: cf es el coeficiente de fricción, u∗ es la velocidad de corte en el lecho, el cual es
definido usando el esfuerzo cortante y el coeficiente de fricción cf

τ0
r
u∗ = (2.33)
ρ
2g
El coeficiente de fricción cf es calculado usando el coeficiente de Chezy Ch cf = Ch2
. Ese
valor también puede ser calculado usando el coeficiente de Strickler, Manning o Niku-
radse. El valor de ck y cε son determinados desde el flujo normal en el centro del canal
principal (TELEMAC-2D, 2001). Todas las constantes adimensionales cµ ,σk ,σε ,c1ε ,c2ε
son basados en casos de test clásicos. Un test prueba sobre la disminución libre de la
turbulencia permite encontrar un valor de c2ε (TELEMAC-2D 2001). La constante cµ
y c1ε son obtenidos desde la data para un flujo turbulento cerca de la pared sólida. La
constante σk y σε han sido (optimizados) en base al rendimiento de los casos de prueba.

2.5. Transporte de sedimentos

2.5.1. Características físicas del sedimento

Las características físicas del sedimento tienen un papel importante en la determinación del
movimiento de las partículas, es decir, el movimiento incipiente. Las partículas de sedimento
se pueden caracterizarse de acuerdo a su tamaño, forma, uniformidad, velocidad de asenta-
miento, cohesión y ángulo de reposo.
Las dimensiones de las partículas de sedimento cubren un gran rango, que a su vez es usado
 CAPÍTULO 2. MARCO TEÓRICO 27

para clasificar a los sedimentos (Yang & Song, 1979). El tamaño del sedimento es comúnmen-
te determinado por análisis granulométrico. Lo más común para la clasificación del sedimento
es de la Unión Geofísica Americana, que se basa en el diámetro aproximado de las partículas
de sedimento (Armitage & McGahey, 2003). Una amplia variedad de "formas"de particulas
de sedimentos se observan en los ríos. Van Rijn et al. (1993) presentó el llamado "factor
de forma"para tener en cuenta la desviación esférica. Naturalmente el cuarzo, partícula que
tiene típicamente el factor forma de 0.7.
La distribución del tamaño del sedimento es determinada por el análisis granulométrico para
determinar el tamaño del sedimento.
La uniformidad del sedimento, se determina a partir de la frecuencia de histogramas que re-
presenta el porcentaje de material que pasa por un diámetro de tamiz, pero retenido en otro
tamiz de menor tamaño. Parámetros estadísticos tales como la media, desviación estándar,
asimetría y curtosis pueden determinarse a partir de esto (Van Rijn et al., 1993).
El diámetro medio d50 es el diámetro que no excede en un 50 por ciento del material por
masa y se expresa así:
d15,9 + d84,4
d50 = (2.34)
2
El diámetro medio (d) es la media aritmética del tamaño de la partícula y se puede expresar
en términos de la mediana del diámetro medio di de cualquier fracción de tamaño de partícula
y el porcentaje de masa de esa fracción pi de la siguiente manera:
P
pi d i
d= (2.35)
100

La desviación estándar σ es determinada usando el diámetro medio ecuación (2.36). La

desviación estándar geométrica σg es calculado usando la ecuación (2.37), el cual asume la
distribución normal de las partículas, lo que es aproximadamente cierto para una clasificación
fina o arena uniforme. Si σg <1.3, el sedimento es considerado uniforme, si σg >1.3, el
sedimento es considerado no uniforme,(Melville & Raudkivi, 1996).

pi (di − d)2
P
2
σ = (2.36)
100
s
d84,4 ∼ 1 d50 d84,4
σg = = ( + ) (2.37)
d15,9 2 d15,9 d16
 CAPÍTULO 2. MARCO TEÓRICO 28

La velocidad de asentamiento de la partícula de sedimento es a veces llamada velocidad de

caída que generalmente es definido como la velocidad terminal de la partícula en un cuerpo de
agua, en general depende: del tamaño, forma, rugosidad superficial y densidad de la partícula
de sedimento; también de la densidad y la viscosidad del fluido. Además del régimen del fluido
tiene influencia en la velocidad de asentamiento.
Si la partícula es menos densa que el fluido, puede tener velocidad de asentamiento negativo
algunas veces llamado velocidad de subida (Raudkivi & Sutherland, 1981); se calcula la
velocidad de asentamiento vss usando el coeficiente de arrastre CD .
De acuerdo a Van Rijn et al. (1993) la velocidad de asentamiento de las partículas no esféricas
con: diámetro característico d, gravedad especifica s, viscosidad cinemática ν y la aceleración
de la gravedad g, pueden ser estimados de la siguiente forma:

(s − 1)gd2
vss = 1 < d ≤ 100 µm (2.38)
18ν
r
10ν 0,01(s − 1)gd3
vss = [ 1+ − 1] 100 < d ≤ 1000 µm (2.39)
d 18ν 2

vss = 1,1[(s − 1)gd]0,5 d > 1000µm (2.40)

Cheng (1997), aproximadamente la velocidad de asentamiento vss de la partícula de arena te-

niendo en cuenta la densidad relativa sumergida ∆ y el diámetro adimensional de la partícula
d∗ es: p
(25 + 1,2d2∗ − 5)(1,5)
vss = ν( ) (2.41)
d
∆g
d∗ = ( 2 )(1/3) d50 (2.42)
ν
ρs − ρ
∆= (2.43)
ρ
Generalmente la velocidad de caída es significativamente pequeña en altas concentraciones.
El efecto de la concentración en la velocidad de asentamiento es conocido como: asentamiento
obstaculizado, (Raudkivi & Sutherland, 1981). El sedimento puede ser dividido principalmen-
te en dos grupos; cohesivo (limo, arcilla) y no-cohesivo (arena).
El peso de las partículas es la fuerza determinante en la transportabilidad de sedimentos
 CAPÍTULO 2. MARCO TEÓRICO 29

no-cohesivos. Por otro lado, para sedimentos cohesivos, la fuerza fisicoquímica juega un rol
importante. Las propiedades del suelo permite que la partícula de sedimento se comporte
como una gran masa. Vanoni & Nomicos (1959) describen que una vez suspendido la masa
cohesiva se rompe y las partículas se comportan como materiales no cohesivos.
El ángulo de reposo φr , algunas veces llamado el ángulo interno de fricción, es la máxima
pendiente lateral que puede ser sostenido la partícula de sedimento que se encuentran en
el fondo o en los lados del canal antes de que comience el deslizamiento. (Van Rijn et al.,
1993); sin embargo, φr para cualquier material que no sea arena cuarcita debe determinarse
experimentalmente.

2.5.2. Fuerzas que actúan en la partícula de sedimento

En realidad, una partícula de sedimento en el lecho del canal está rodeada de partículas
vecinas que descansa sobre el lecho antes del movimiento incipiente. Tres fuerzas principales
actúan sobre la partícula expuesta: (fuerza de arrastre, suspensión y gravedad). La Figu-
ra (2.4) muestra las diferentes fuerzas que actúan sobre una partícula de sedimento (lecho
inclinado).

Figura 2.4: Fuerzas en la partícula de sedimento con lecho inclinado.

La fuerza de arrastre FD , la fuerza de sustentación FL y la fuerza de presión y al esfuerzo

cortante que actúan sobre la partícula.
El esfuerzo cortante en cualquier punto de la superficie de la partícula actúa tangencialmente
 CAPÍTULO 2. MARCO TEÓRICO 30

en la dirección del flujo y depende de la gradiente de velocidad local. La fuerza de arrastre

debido al esfuerzo cortante a menudo se denomina resistencia superficial.

∂u ∂v
τxy = (µ + µxy(t) )( + ) (2.44)
∂y ∂x

En la ecuación (2.44) µ es el coeficiente de viscosidad debido a que el esfuerzo cortante actua

en el flujo y µxy(t) se debe al flujo turbulento. La fuerza de presión en cualquier punto de
la superficie de la partícula que actúa normal a su superficie depende de la profundidad y
la velocidad local de acuerdo con la ecuación de Bernoulli. Si z es la altura sobre alguna
referencia, la ecuación de Bernoulli puede ser escrito como en la ecuación (2.44). La presión
en la profundidad del agua p = ρgy se puede determinar usando las condiciones hidrostáticas.

p u2
+ u + z = constante (2.45)
ρg 2g

La ecuación (2.45) indica claramente que para una referencia constante, el aumento en la
velocidad conducen a una disminución correspondiente en la presión local. Esta situación se
complica por la separación de flujo, que resulta en la formación de remolinos, Según Armitage
& McGahey (2003), los remolinos redistribuyen la presión de manera más uniforme sobre las
superficies separadas. Van Rijn et al. (1993) muestran que la integración de las fuerzas de
presión sobre toda la superficie de la partícula indica la fuerza de arrastre resultante, que se
conoce como: arrastre de forma.
La fuerza gravitacional FG está relacionada con el peso sumergido de la partícula. A pesar
de la ausencia del movimiento del fluido, hay una fuerza neta en la partícula y su dirección
depende del efecto de flotabilidad. La fuerza gravitacional actúa verticalmente hacia abajo si
la partícula es más densa que el agua y verticalmente hacia arriba si no lo es. Van Rijn et al.
(1993) describen la fuerza gravitatoria FG de una partícula crítica en función de un diámetro
de partícula d, densidad de partícula ρs y densidad de agua ρ:

1
FG = π(ρs − ρ)gd3 (2.46)
6

Chien & Wan (1999) describen la fuerza de arrastre FD y la fuerza de sustentación FL

usando el área proyectada ap (que se define diferente para la fuerza de arrastre y la fuerza de
sustentación), coeficiente de arrastre CD , coeficiente de sustentación CL , densidad del agua
 CAPÍTULO 2. MARCO TEÓRICO 31

y velocidad efectiva cerca de la partícula del lecho u0 como:

ρu20
FD = CD ap (2.47)
2
ρu20
FL = CL ap (2.48)
2
La presencia de una frontera en el canal hace que la distribución de la velocidad no sea lineal.
También la fuerza de sustentación es relativamente mayor en la frontera donde las fuerzas de
arrastre no actúan en el centro de la gravedad de la partícula, según Saffman (1965), para un
flujo viscoso, la fuerza de sustentación debido al movimiento de giro es menor que el efecto
de corte y es generalmente despreciado Colebrook & White (1937).

2.5.3. Papel de la rugosidad efectiva

La rugosidad efectiva k es muy importante para la simulación. En el caso de un lecho móvil,

es la suma de la rugosidad del grano (ksg ) generada por la fricción superficial y la rugosidad
de forma (ksf ) generada por la fuerza de presión que actúan sobre el lecho.

ks = ksg + ksf (2.49)

Según Van Rijn et al. (1993), se pueden usar dos enfoques para estimar la rugosidad del
lecho; parámetros relacionados con la forma (por ejemplo, longitud de la forma del lecho,
altura, inclinación y material del lecho) o parámetros integrales (por ejemplo, profundidad
media, velocidad media).
Según Young et al. (2010), la partícula en la frontera podría formar parte de una frontera
laminar (Res <5), frontera de transición (5 <Res <70) o frontera turbulenta (Res >70),
donde Res , es el número de Reynolds de la rugosidad.

u∗ ks
Res = (2.50)
ν

Dónde: ν es la viscosidad cinemática y u∗ es la velocidad de corte que se define utilizando el

esfuerzo cortante del lecho τ0 y la densidad del fluido ρ:

τ0
r
u∗ = (2.51)
ρ
 CAPÍTULO 2. MARCO TEÓRICO 32

Colebrook & White (1937) muestran cuando Res , aumenta más allá de 3.5, comienza la
separación detrás de la parte superior de la partícula y se forma una estela. De acuerdo con
Chien & Wan (1999), la fuerza de arrastre se convierte entonces en la fuerza dominante.
Otras definiciones de ks , son Yang & Song (1979), ks = d65 ; (Einstein et al., 1950), ks = d90
Meyer-Peter & Müller (1948) ks = d85 ; Simons & Richardson (1966). El valor de k disminuye
significativamente al aumentar la velocidad, probablemente porque las formas del lecho se
vuelven más redondos o se lavan (Van Rijn et al., 1993).

2.5.4. Inicio del movimiento del sedimento

El inicio del movimiento del sedimento se puede determinar a partir de las fuerzas que actúan
sobre las partículas debido a las fuerzas de arrastre, sustentación y la gravedad.
Van Rijn et al. (1993) describen el inicio del movimiento de sedimento utilizando las fuerzas
que actúan sobre la partícula de sedimento en el lecho como se muestra en la Figura (2.5)

Figura 2.5: Fuerzas en la partícula de sedimento

Fuente: Van Rijn et al. (1993)

Van Rijn et al. (1993) descuidan la relación de la fuerza de sustentación y el peso de la

partícula sumergida FL /FG = 0. Introduce el coeficiente empírico α1 ,que tiene en cuenta el
efecto de la sustentación. La condición de movimiento del sedimento, esta dado por:
 CAPÍTULO 2. MARCO TEÓRICO 33

b2
FD ≥ α1 FG tanϕ α1 = (2.52)
b1 + b2
El valor del coeficiente empírico α1 depende de la turbulencia local y, por lo tanto, del número
de Reynolds Van Rijn et al. (1993).
Yang & Song (1979) agrega más simplificaciones y muestra el movimiento incipiente del
sedimento que satisface una o más de las tres condiciones indicadas en la ecuación (2.53) y
la Figura (2.6).

FL = FG (2.53)

FD = FR (2.54)

MO = MR (2.55)

Dónde: FL es la fuerza de sustentación, FG es la fuerza de gravedad, FD es la fuerza de

arrastre, FR es el fuerza de resistencia, M0 es el momento de vuelco debido a FD ,FR y MR
es el momento de resistencia debido al FL y FG .

Figura 2.6: Fuerzas que actuan en la partícula de sedimento por Yang(1996)

Fuente: Yang & Song (1979)

En realidad, el cálculo de las fuerzas que actúan sobre una partícula es difícil. Muchos investi-
gadores, definen los criterios de movimiento incipiente del sedimento usando esfuerzo cortante
 CAPÍTULO 2. MARCO TEÓRICO 34

Shield et al., (1936), velocidad (Yang & Song, 1979). Fuerza de la corriente (Bagnold, 1966)
o probalilidad (Gessler, 1970).
Shields (1936) fue uno de los primeros investigadores en darse cuenta de la dificultad del
cálculo y análisis de las fuerzas que actúan sobre una partícula de sedimento. Él utilizó el
análisis dimensional para establecer su conocido diagrama de Shields para un movimiento
incipiente de una partícula usando el parámetro adimensional, por ejemplo:número Reynolds
de la partícula Re∗ , y el parámetro de Shield θ.

u∗d
Re∗ = (2.56)
ν
τ0 ρu2∗
θ= = (2.57)
(ρs − ρ)gd (ρs − ρ)gd
Dónde: u∗ velocidad de corte, d diámetro de la partícula, ρs es la densidad del sedimento, ρ
es la densidad del agua y g es la aceleración de la gravedad.
Shields (1936) midió la tasa de transporte en varios valores de θ al menos dos veces más grande
que el valor crítico θc , y luego extrapoló linealmente al punto cero. Este método indirecto se
utilizó para evitar la dificultad de determinar la condición de precisión en la que una partícula
de sedimento debería moverse. Si τc es el esfuerzo cortante crítico de cizallamiento del lecho
para el movimiento incipiente, luego el parámetro crítico Shields θc , está dado por:

τc
θc = (2.58)
(ρs − ρ)gd50

La Figura (2.7) muestra el parámetro Shields como una función del número de Reynolds de
la partícula. La línea sólida indica la extensión de datos por Shields, (Raudkivi & Sutherland,
1981).
 CAPÍTULO 2. MARCO TEÓRICO 35

Figura 2.7: Parámetro de Shields con respecto al número de Reynolds de la partícula.

Fuente:Raudkivi & Sutherland (1981)

Yang & Song (1979) propusieron una velocidad crítica adimensional U c/νss para describir
el movimiento incipiente del sedimento usando el número de Reynolds de la partícula Re∗ .
Mostró una variación en línea recta o parabólico de la velocidad critica adimensional U c/νss
en regimen hidráulicamente liso (Re∗ < 5) y en regimen completamente rugoso (Re∗ > 70) y
en régimen de transición (5 < Re∗ < 70), en ambos la contribución de la fricción laminar y
la fricción turbulenta deben ser considerados.
En base a los datos de campo, Yang & Song (1979) sin embargo, extendió la curva parabólica
para este régimen teniendo en cuenta los efectos de coeficientes empíricos con respecto a Re∗ .
La curva del movimiento incipiente Yang & Song (1979) se muestra en Figura (2.8).
Debe notarse que Yang & Song (1979) ignoró el efecto de la rugosidad del flujo (y/d) para el
régimen de transición debido a la falta de disponibilidad de datos. Talapatra & Ghosh (1983)
independientemente verificaron los criterios de Yang & Song (1979).
 CAPÍTULO 2. MARCO TEÓRICO 36

Figura 2.8: Variación de las dimensiones de la velocidad critica.

Fuente:Yang & Song (1979)

2.5.5. Formas de transporte de sedimentos

Una vez que se ha iniciado el movimiento de las partículas de sedimento. La partícula co-
mienza a moverse, en un canal aluvial. Se observan tres tipos diferentes de movimiento de
partículas Van Rijn et al. (1993); deslizando, rodando, saltando.
En la mayoría de los casos, la partícula de sedimento rodará o se deslizará antes de que sean
llevados a suspensión.
El transporte de partículas por laminación, deslizamiento y saltación en las proximidades del
lecho se denomina colectivamente carga de fondo del lecho, que se refiere a los sedimentos
que se mueven continuamente en contacto con la frontera, es decir, deslizando y rodando. La
saltación es cuando las partículas saltan siguiendo trayectorias de tipo parabólico cerca del
lecho.
Las partículas suspendidas son transportadas por las corrientes turbulentas y permanecen en
suspensión durante un período apreciable. Según Yang (1996), los sedimentos se transportan
principalmente a través de la carga suspendida en la mayoría de los ríos naturales mientras
 CAPÍTULO 2. MARCO TEÓRICO 37

que el transporte de fondo del lecho es solo alrededor del 5 % dependiendo de la clase del
sedimento.
Según Raudkivi & Sutherland (1981), la concentración de sedimentos en suspensión (volu-
men sólido) disminuye con la distancia desde el lecho. Esta tasa de disminución depende de
la velocidad de caída vss y de la velocidad de corte del lecho u∗ . El cual se usa para definir
el número de movilidad M n = u∗ /vss . Propone una regla empírica para la predicción de
diferentes tipos de transporte que utiliza el número de movilidad, por ejemplo:

Rodando y/o deslizándose 6 > (1/M n) > 2.

Saltando 2 > (1/M n) > 0,6.

Suspensión 0,6 > (1/M n) > 0.

Graf & Altinakar (1998) distinguen diferentes tipos de transporte de sedimentos en términos
de mezcla agua-sedimento, dependiendo de la concentración de sedimento Cs (volumen sólido
por unidad de volumen de agua) o concentración de sedimento △ρ (masa sólida por unidad
de volumen del agua), se divide en tres partes:

La carga del lecho y el transporte de carga en suspensión son comunes en una mezcla
newtoniana (Cs << 1 % o △ρ << 16kg/m3 ) que se encuentra con mayor frecuencia
en los cursos de agua.

En una mezcla cuasi newtoniano (C < 8 % o △ρ < 130kg/m3 ) usualmente se

observa la concentración en suspensión cerca del lecho.

Las corrientes de turbidez, por ejemplo, entran en esta categoría. Mezclas no newtonia-
nas (Cs > 8 % o △ρ > 130kg/m3 ) se asocian con las suspensiones hiper-concentradas
experimentadas con corrientes de turbidez y flujo de escombros.

Clásicamente, la suma de la carga de fondo y carga suspendida se conoce colectivamente

como carga total. Dependiendo de la fuente de material de sedimentos. Raudkivi & Sutherland
(1981) descompone aún más la carga total en carga de material del lecho y carga de lavado. En
general, la carga de lavado consiste en materiales que son más finos que los que se encuentran
en el lecho del río y casi nunca están en contacto con el lecho. La cantidad de carga de lavado
depende principalmente del suministro de la cuenca (hidrología), no de la hidráulica del río.
 CAPÍTULO 2. MARCO TEÓRICO 38

2.5.6. Fórmulas de transporte de sedimentos no cohesivo

Está asociado con el inicio del movimiento de sedimentos y continúa mediante el transporte
de sedimentos debido a las fuerzas hidrodinámicas a menos que la condición de equilibrio
se alcance. Se han desarrollado una gran cantidad de fórmulas semi-empíricas para calcular
la tasa de transporte sólido como una función de variables hidrodinámicas (tirante, fricción
de fondo, olas, etc.) y de las variables del sedimento como: d50 , d90 , densidad relativa, etc.
Raudkivi & Sutherland (1981).
La fórmula de transporte puede estimar solo la carga del lecho o carga suspendida o carga
total. Las fórmulas se pueden agrupar de acuerdo a su enfoque teórico, por ejemplo: (el
esfuerzo cortante, la energía de pendiente, probabilidad, potencia de las líneas de corriente
etc). Las fórmulas utilizadas en este estudio se analizan a continuación.

2.5.6.1. Meyer-Peter and Muller

La fórmula de Meyer-Peter y Muller (1948) se basa en el enfoque de la pendiente de energía y

es limitado al transporte de “carga del lecho". Llevaron a cabo extensos trabajos experimen-
tales en un canal de laboratorio con una sección transversal de 2x2m2 y una longitud de 50m.
El material del lecho que se uso fue una mezcla de partículas que tiene las siguientes caracte-
rísticas: (d = 0.4 a 29 mm, pendiente, S 0.0004 a 0.02, profundidad 0.1 a 1.2 m). De acuerdo
con Meyer-Peter & Müller (1948). El transporte de carga del lecho qb puede expresarse como
Van Rijn et al. (1993):

qb
Φb = p 3
= 8(µb θ − 0,047)(1,5) (2.59)
(s − 1)gd50

Ch log(12h/ks ) (1,5)
µb = [ ]=[ ] (2.60)
Ch
′
log(12h/d90 )
Donde:

Φb es la función de carga de fondo adimensional del sedimento.

θ es el parámetro Shields.

µb es el factor de forma del lecho o el factor de eficiencia.

 CAPÍTULO 2. MARCO TEÓRICO 39

d es el diámetro medio de la partícula.

Ch es el coeficiente de Chezy.

Ch′ es el coeficiente de Chezy relacionado con el grano.

ks es la aspereza efectiva del lecho y h es la profundidad de agua.

El valor constante 0.047 puede interpretarse como el parámetro de Shields crítico θc porque
en el experimento considera material grueso. Meyer-Peter y Muller (1948) encontraron el
valor de d = Σpi di puede ser cerca aproximadamente 1.1 a l.3 veces más grande que d50
para material casi uniforme. El diámetro medio de la partícula d50 también puede usarse por
ejemplo como d, sin embargo la rugosidad del grano está relacionado con d90 . De acuerdo
con Van Rijn et al. (1993). La fórmula muestra que una variación del 25 % del diámetro
de la partícula (d = 800 ± 200µm) resulta un 10 % de variación de la tasa de transporte
(qb = 4,1x10−5 ± 0,4x10−5 m2 /s) para arena uniforme, el coeficiente de Chezy relacionado con
el grano Ch′ es a menudo aproximado en términos de d50 así que usando la relación d90 = 3d50
(SISYPHE 2004).
Cabe señalar que el valor de Φh será cero si el valor de µb θ < θc = (0,047).SISYPHE (2004)
recomienda el uso de la fórmula de Meyer-Peter y Muller para los sedimentos finos a gruesos
en el rango 0,4 < d50 < 29mm.

2.5.6.2. Einstein y Brown

Einstein et al. (1950) fue uno de los primeros investigadores en estimar el transporte de sedi-
mentos utilizando el concepto de probabilidad. Las dificultades para definir el criterio crítico
de movimiento incipiente y la necesidad de tener en cuenta las fluctuaciones de flujo turbu-
lento en el transporte de sedimentos en lugar de las fuerzas que actúan sobre las partículas
de sedimento. Como resultado expresó la probabilidad del comienzo y cese del movimiento
del sedimento. En su experimento destacaron los siguientes puntos (Yang & Song, 1979):

Hay un intercambio constante e intenso de partículas entre el material del lecho y carga
del lecho.
 CAPÍTULO 2. MARCO TEÓRICO 40

La tasa de erosión depende del número y las propiedades de las partículas en unidad
de área y la probabilidad de la fuerza de sustentación hidrodinámica instantánea en la
partícula sea grande, lo suficiente para moverlo. Por otro lado, la tasa de deposición
por unidad de área del lecho unitario depende de la tasa de transporte más allá de
una sección determinada. Como, ver la probabilidad de que la carga hidrodinámica
permitirá que la partícula se deposite.

El enfoque de Einstein et al. (1950) se basa en la condición de lecho estable, es decir, la tasa
de deposición debe ser igual a la tasa de erosión. En consecuencia, el número de partículas de
tamaño d erosionado por unidad de tiempo y unidad de área del lecho se calcula utilizando
la probabilidad de que cualquier partícula dada se erosione. Al mismo tiempo, la cantidad de
partículas, que se depositarán por unidad de área y unidad de tiempo son calculados usando
la relación de la descarga del lecho. Por lo tanto, correlaciona esta erosión y deposición
teniendo en cuenta la fuerza de sustentación y el peso sumergido de las partículas para
obtener el transporte de carga del lecho. De acuerdo con Yang & Song (1979). Einstein et al.
(1950) usó tres diferentes figuras gráficas para calcular la carga del lecho.
Brown (1950) desarrolló aún más el enfoque de Einstein et al. (1950) teniendo en cuenta el
parámetro de Shields θ y la característica adimensional del diámetro de sedimento d∗ , Yang
& Song (1979).
Einstein et al. (1950) propusieron la siguiente ecuación para calcular la carga del lecho:

qb
Φb = p = f1 (d∗ )f2 (θ) (2.61)
g(s − 1)d350

2 36 36
f1 (d∗ ) = ( + )( 0,5) − ( )0,5 (2.62)
3 d∗ d∗

f2 (θ) = 2,15exp(−0,391/θ) θ ≤ 0,2 (2.63)

f2 (θ) = 40θ3 θ > 0,2 (2.64)

SISYPHE (2004) recomienda usar la fórmula de Einstein et al. (1950) para sedimentos rela-
tivamente gruesos en el rango 0,2 < d50 < 3,0mm
 CAPÍTULO 2. MARCO TEÓRICO 41

2.5.6.3. Engelund y Hansen

Bagnold (1966) fue uno de los primeros investigadores en desarrollar una función de transpor-
te de sedimentos teniendo en cuenta el efecto de las líneas de corriente Engelund & Hansen
(1967) luego aplicaron el enfoque de Bagnold y utilizaron el principio de similitud para obte-
ner una función adimensional de carga total de transporte de sedimento total adimensional
Φt .
El concepto de balance de energía Wr = Wd también se considera teniendo en cuenta la
altura de la forma del lecho y la longitud de la forma del lecho, Wr es el trabajo por unidad
de tiempo y ancho requerido para elevar la carga de sedimentos sobre una altura igual a
la altura del lecho. Wd es el trabajo por unidad de tiempo hecho por el fluido al mover las
partículas sobre una longitud igual a la forma del lecho.
Engelund & Hansen (1967) utilizaron 100 conjuntos de datos con el parámetro crítico de
Shields de 0.06 ( partículas finas) para establecer la relación entre la tasa de transporte y
el parámetro Shields θ. La fórmula calcula la carga total qt , puede expresar como Van Rijn
et al. (1993):
0,05U 5
qt = (2.65)
(s − 1)2 g 0,5 d50 Ch3
Donde: U es la velocidad promediada en la profundidad, s es la densidad específica, Ch es
el coeficiente Chezy, g es la aceleración de la gravedad, y d50 es el diámetro medio de la
partícula de sedimento. La ecuación (2.66) es la fórmula original de transporte de sedimentos
dada por Engelund & Hansen (1967).
Para lograr consistencia con las formulaciones del modelo y el cálculo de la fricción del fondo.
SISYPHE (2004) simplifica la fórmula original de Engelund & Hansen (1967) en términos de
la tasa de transporte no-dimensional Φt como:
qt θ2,5
Φt = p = 0,05 (2.66)
g(s − 1)d350 cf
Donde: el coeficiente de fricción cf puede ser calculado usando varios coeficientes dados en
SISYPHE(2004).
Para canales aluviales, (Chollet Cunge, 1980) modificaron aún más la fórmula de Engelund &
Hansen (1967) para explicar los efectos de la formación de dunas utilizó la fricción superficial
adimensional, basado en las pruebas de validez del modelo.
 CAPÍTULO 2. MARCO TEÓRICO 42

SISYPHE (2004) recomienda usar la fórmula de Engelund & Hansen (1967) para sedimentos
finos en el rango 0,2 < d50 < 1mm.

2.5.6.4. Van Rijn

Van Rijn et al. (1993) llevaron a cabo 130 experimentos con diámetros de partículas, d50
por lo que van desde 200 a 2000µm. Para profundidades de agua mayores a 0.1 y número de
Froude menor a 0.9. La influencia de la rugosidad de la pared lateral se eliminó utilizando
el método de Vanoni & Nomicos (1959). Usaron un factor de forma del lecho para eliminar
la rugosidad causada por la forma del lecho. Usando las características adimensionales del
diámetro de la partícula d∗ y el parámetro de esfuerzo de corte excesivo T∗ (Ecuación 2.68).
la carga del lecho se calcula como:

1,5 2,1 −0,3

qh = 0,053(s − 1)0,5 g 0,5 d50 T∗ d∗ T∗ < 3 (2.67)

1,5 1,5 −0,3

qh = 0,053(s − 1)0,5 g 0,5 d50 T∗ d∗ T∗ ≥ 3 (2.68)
θ
T∗ = −1 (2.69)
θc
Dónde: s es la gravedad específica, g es la aceleración de la gravedad. ν es la viscosidad
cinemática, θ es el parámetro de Shields y θc el parámetro crítico de Shields.
Capítulo ,,,

Modelo hidráulico del río Madre de Dios

Resumen:
En este capítulo se describe el modelo hidráulico del río Madre de Dios construido en el
Laboratorio Nacional de Hidráulica. Además, los pasos realizados para representar el lecho y
medir los datos de: tirante, velocidad y batimetría del lecho.

3.1. Similitud entre el prototipo y modelo hidráulico

Las investigaciones en un modelo hidráulico sirven para llegar a una solución de un proble-
ma específico, si hay una similitud entre el modelo y el prototipo. Teniendo en cuenta las
características de los fenómenos hidráulicos por representar y las influencias de las fuerzas
gravitacionales e inerciales. Para mantener la similitud entre el modelo y prototipo se ha
adoptado la Ley de Froude. Además, para la adopción de la escala se tuvo en cuenta los
siguientes requerimientos:

Caudales máximos y mínimos.

La extensión del prototipo a reproducirse en el modelo.

Las condiciones limites elegidas que deben permitir una representación fiel de los fenó-
menos (niveles y velocidades).

La Figura 3.1 y la Figura 3.2 muestran el modelo hidráulico y el prototipo del tramo La
Pastora en el río Madre de Dios, respectivamente.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 44

Figura 3.1: Modelo hidráulico del río Madre de Dios - Tramo la Pastora construido en el
Laboratorio Nacional de Hidráulica.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 3.2: Prototipo del río Madre de Dios - Tramo La Pastora.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 45

La escala de longitud fue 1/60, la cual permitió obtener los patrones de flujo y otras ca-
racterísticas del fenómeno a modelar sin distorsión. A partir de está escala se derivaron las
relaciones mostradas en la Tabla 3.1

Tabla 3.1: Relación de escalas Modelo - Prototipo

Magnitud Factor de escala Valor

Longitud Lp /Lm =L 60
Velocidad Vp /Vm =L0,5 7.746
Tiempo Tp /Tm =L0,5 7.746
Caudales Qp /Qm =L5/2 27885.48
Rugosidad np /nm =L1/6 1.979
Transporte de sedimentos Tp /Tm =Lx Ly Lγ 7.746
Descarga unitaria de sedimentos qp /qm =Lγ Lx Ly /(Tsp /Tm ) 464.76
Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
Nota: El subíndice ′ p′ representa al prototipo y el subíndice ′ m′ representa al modelo.

Donde:

Lp , Lm Longitudes en prototipo y modelo, respectivamente.

Vp , Vm Velocidades en prototipo y modelo, respectivamente.

Tp , Tm Tiempos en prototipo y modelo, respectivamente.

np /nm Rugosidad en prototipo y modelo, respectivamente.

Tsp /Tsm Tasas de transporte de sedimentos en prototipo y en modelo, respectivamente.

Lx Ly Lγ Relaciones de escala en la dirección x, en la dirección y, y relación de pesos

específicos (prototipo - modelo).

En 2017 se reconstruyó el lecho del modelo hidráulico del río Madre de Dios en base a la
relación de escalas modelo - prototipo las cuales se muestran en la Tabla 3.1. Para estudiar
el lecho del río Madre de Dios en el modelo hidráulico se realizo lo siguiente:
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 46

Se representó el lecho móvil de la zona de estudio a escala no distorsionada de 1/60.

Se calibró el modelo hidráulico teniendo en cuenta los niveles de agua registrado en el

río Madre de Dios.

3.2. Descripción del modelo hidráulico del río Madre de

Dios - Tramo la Pastora
A partir del informe de hidrológia (Ingeníeria de Detalle - Defensas Ribereñas Sector la
Pastora km 480+070 – km 480+910), y usando la relación de escalas que se muestra en la
Tabla 3.1, se obtuvo el caudal máximo y mínimo del modelo hidráulico, los que fueron usados
en la etapa experimental por medio del suministro de 4 bombas.

L
Caudal máximo (Qmáx) = 543 s
.

El lecho móvil del modelo hidráulico del tramo La Pastora del río Madre de Dios fue represen-
tado en el Laboratorio Nacional de Hidráulica a escala 1/60 con las siguientes características
hidráulicas:

Longitud (L) de 50 m.

Ancho medio (B) de 6.9 m.

Radio de curvatura (R) de 13.8 m.

Tirante máximo (y) de 0.25 m.

R
Por lo tanto, la relación B
resultó 2.
Para la etapa experimental se ultilizó cuatro bombas de eje vertical que permitieron sumi-
nistrar caudales de hasta 550 L/s. Este caudal fué alimentado al modelo a través de tuberías
de PVC de 356 mm de diametro (14 pulgadas).
Se realizó 6 experimentos para obtener la relación esperada entre el prototipo y el modelo hi-
dráulico de las distintas variables: tirante, velocidad y batimetría. Éstas relaciones obtenidas
fueron empleadas posteriormente para calibrar y validar el modelo numérico.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 47

3.2.1. Detalles de replanteo del lecho del río Madre de Dios en el

modelo hidráulico para el Tramo la Pastora

La información proporcionada por la empresa ODEBRECHT (Proyecto N 0 358-PNICP-

PIAP-2014, auspiciado por INNOVATE) el modelo hidráulico del río Madre de Dios, cuenta
con un total de 94 secciones transversales a lo largo del río desde la progresiva 0+120 hasta
la progresiva 4+320. Las distancias medidas desde el eje del río hacia la margen derecha son
positivas y desde el eje hacia la margen derecha son negativas, además las secciones están
separadas 1m y en la zona que describe el tramo curvo a cada 0.50 m. También se tuvo en
cuenta las cotas de cada punto y el punto de referencia BM1 = 156.891 msnm, entre otros,
por lo que fue necesario gran minuciosidad al ejecutar el replanteo del lecho.
Los pasos seguidos para replantear al lecho en el modelo hidráulico fueron:

a) Instalar el nivel de ingeniero a una distancia adecuada para barrer las progresivas de las
secciones a trabajar durante el día.

b) Seleccionar el BM más próximo de la zona a trabajar. Los BMs se ubican en la margen

derecha del modelo físico. La mira siempre se sostiene con la ayuda de un nivel cicular
(ojo de pollo).

c) Marcar las progresivas sobre los muros laterales del modelo, de tal manera que se pueda
ubicar la progresiva de la sección a trabajar, y luego colocar la regla graduada de aluminio
sobre los muros.

d) Usando la varilla de aluminio se procedió a marcar los puntos (cotas) de la sección, a partir
del eje donde los puntos ubicados a la derecha del eje son positivas (+) y a la izquierda
negativas (-) de acuerdo con los datos proporcionados.

3.2.2. Granulometría del material del lecho

Para determinar la granulometría del lecho a la escala del modelo se tuvo en cuenta que
esta no puede ser muy pequeña porque puede dejar de comportarse como material granular
y presentar fuerzas intergranulares de cohesión. Este fenómeno ocurre cuando D<0.1mm.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 48

Por esta razón, se tomó parte de la curva granulométrica del material del lecho del rio
que cumple con esa condición, que se muestra en la Tabla 3.2 y la Figura 3.2.

Los diámetros característicos, se obtuvieron de Figura 3.3, siendo estos:

D16 = 0.220 mm

D50 = 0.400 mm

D84 = 0.550 mm

D90 = 0.590 mm

Tabla 3.2: Porcentaje que pasa o es retenido por los tamices.

Tamiz Abertura de % parcial Acumulado
- Malla [mm] Retenido Retenido Pasa
”
1 1/2 38.1 - - 100
1” 25.4 20 20 80
3/4” 19.05 9.1 29.1 70.9
”
1/2 12.7 7.5 36.6 63.4
3/8” 9.525 3.5 40.1 59.9
”
1/4 6.35 3.7 43.8 56.2
N0 4 4.76 2.1 45.9 54.1
N 0 10 2 1.8 47.7 52.3
N 0 20 0.84 0.2 47.9 52.1
N 0 30 0.59 1.7 49.6 52.1
0
N 40 0.426 1.6 51.2 48.8
N 0 60 0.25 23.5 74.7 25.3
0
N 100 0.149 19 93.7 6.3
N 0 200 0.074 5.7 99.4 0.6
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 49

Figura 3.3: Curva granulométrica del modelo hidráulico del rio Madre de Dios.

d50

0.4

Fuente: Proyecto N0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 50

3.2.3. Caudal y tasa de transporte

El caudal con que se realizó la etapa experimental fue de 490 Ls . La tasa de transporte de
kg
sedimento que se usó fue de 0.01 s

3.3. Datos obtenidos durante los ensayos

Los ensayos en el modelo hidráulico se realizaron en febrero del 2017 con una duración de
un día, donde se registraron: las velocidades en las tres direcciones (x,y,z), el tirante y la
batimetría del lecho al finalizar el ensayo experimental.

3.3.1. Batimetría del lecho

Después de realizar la etapa experimental se realizó el levantamiento de la batimetría del

cauce, por medio de un Escáner Láser el cual permitió capturar una fotografía tridimen-
sional compuesta por millones de puntos. El trabajo con el Escáner Láser consistió en
reconocer el modelo hidráulico y buscar el lugar estratégico para posesionar los puntos de
control para el levantamiento de la información. Luego se procesó la información con el
software del equipo; con estos pasos se obtuvo una nube de puntos que serán exportados
al formato csv. A partir de la nube de puntos se generó un archivo binario de formato
serafín (.slf), que es un archivo de entrada para el modelo numérico,( ver detalle en el
Capitulo 5).
En la Figura 3.4 se muestra la variación del cauce después de la etapa experimental debido
a la interacción entre la fuerza centrifuga causada por la curvatura del río y la gradiente
vertical de la velocidad del flujo principal.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 51

Figura 3.4: Configuración del lecho después de realizar el experimento, se observa el cambio
de forma del lecho tipo duna.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

En la Figura 3.5 muestra la batimetría del lecho después del experimento, capturado por
el Escáner Láser. Cuatro secciones trasversales (s2,s4,s6,y s8) fueron analizados y usados
para validar y calibrar el modelo numérico.

Figura 3.5: Batimetría del lecho después de la etapa experimental, capturado por un Escaner
Laser, la flecha indica la dirección del flujo

Fuente:Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 52

En las Figuras 3.6,3.7,3.8 y 3.9 representan las secciones transversales S2,S4,S6 y S8. Las
sección S6 y S8 muestran la variación del cauce debido al proceso de erosión y sedimen-
tación que se origino por la presencia del flujo secundario, porque las líneas de corriente
de la superficie libre tienen la dirección del flujo y en el lecho tienen dirección opuesta.

Figura 3.6: Sección transversal S2, Ver Figura 3.5.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 3.7: Sección transversal S4, Ver Figura 3.5.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 53

Figura 3.8: Sección transversal S6, Ver Figura 3.5.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 3.9: Sección transversal S8,Ver Figura 3.5.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

3.3.2. Medición de las velocidades

El acopio de datos de velocidades en el modelo se realizó con el velocímetro acústico

biestático (perfilador-vectrino) que puede medir las fluctuaciones rápidas de velocidad.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 54

La medición de velocidad se basa en el efecto Doppler que consiste en enviar un pulso

acústico corto desde un transductor activo. Cuando dicho pulso atraviesa el punto donde
se orientan los transductores pasivos, el eco es registrado por cada uno de los haces. Las
señales son procesadas para determinar sus desviaciones Doppler para luego ajustar la
escala en base a la velocidad de propagación del sonido en el agua dando como resultado
la velocidad en las tres direcciones. Estas velocidades se almacenan en un computador
para luego ser analizados. Se tomaron velocidades en siete secciones transversales, medidas
desde la margen izquierda hacia la margen derecha en seis puntos desde la superficie libre
hasta 5cm encina del fondo del lecho.
En la Figura 3.10 muestra las sección sec1, sec2, sec3, sec4, sec5, sec6 y sec7 donde se
registro la velocidad en las tres direcciones con el velocímetro acústico.

Figura 3.10: Siete secciones transversales en planta y distribución de las velocidades en seis
puntos medidos desde la margen derecha hacia la margen izquierda.

Fuente: Elaboración propia.

A continuación se muestran los pasos realizados para registrar las velocidades en las
secciones (sec1, sec2, sec3, sec4, sec5, s3c6, y sec7).
En la Figura 3.11 muestra la sección sec1 donde se registraron las velocidades con el
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 55

velocímetro acústico medidos desde la margen izquierda a la margen derecha utilizando

un puente corredizo que fué alineado en cada sección transversal.

Figura 3.11: Las velocidades se midieron con Velocímetro Acústico Doppler en las siete sec-
ciones transversales

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

En la Figura 3.12 muestra el Velocímetro Acustico Biestatico que tiene 4 haces, que
permiten capturar los pulsos de las desviasiones en las tres direcciones.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 56

Figura 3.12: Velocímetro Acústico Biestático está colocado a una distancia de 5cm por encima
del lecho.

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

En la Figura 3.13 se muestra una laptop donde fueron almacenados los datos de las
velocidades en las tres direcciones (vx ,vy ,vz ), los que fueron procesadas con el lenguaje
de programación Matlab que consistio en ordenar las variables en las tres direcciones y
exportar en el archivo Excel de extensión csv.

Figura 3.13: Registro y almacenamiento de las velocidades en una laptop para su posterior
procesamiento

Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 57

3.3.2.1. Velocidad media

A partir de los datos de velocidad (vx ,vy ,vz ) medidos en los seis puntos en cada sección
transversal, se calculó la velocidad media por el método del trapecio, tal como muestra la
Tabla 3.3.
A partir de los datos calculados en seis puntos de cada sección transversal haciendo un
total de 42 datos puntuales de velocidad media los cuales fueron usados para calibrar y
validar los resultados del modelo numérico. Los datos de velocidad medidos en el modelo
se muestra en el Anexo 3.

Tabla 3.3: Velocidades medias en seis puntos de cada sección transversal

Secciones pt1 pt2 pt3 pt4 pt5 pt6
v[m/s] v[m/s] v[m/s] v[m/s] v[m/s] v[m/s]
sec1 0.324 0.339 0.350 0.367 0.375 0.281
sec2 0.390 0.423 0.443 0.483 0.461 0.359
sec3 0.409 0.431 0.462 0.443 0.408 0.050
sec4 0.406 0.401 0.394 0.369 0.250 0.027
sec5 0.420 0.351 0.336 0.284 0.234 0.003
sec6 0.479 0.382 0.338 0.334 0.282 0.024
sec7 0.256 0.317 0.296 0.322 0.286 0.081
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

A continuación en las Figuras 3.14, 3.15, 3.16, 3.17, 3.18, 3.19 y 3.20 representan las
secciones sec1, sec2, sec3, sec4, sec5, sec6, y sec7 que son los perfiles de velocidad en el
plano x-y.
Las Figuras 3.14, 3.15 y 3.16 corresponden a las secciones del tramo recto el perfil. Se
observa que la velocidad tiene mayores valores en las proximidades de la parte central de la
sección. A partir de la Figura 3.17 a 3.20 el perfil de velocidad varia debido a la curvatura
del cauce, debido a la variación de la gradiente de velocidad y la fuerza centrífuga.
 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 58

Figura 3.14: Perfil de velocidad de la Sec1 en el plano x-y.

Fuente: Elaboración propia.

Figura 3.15: Perfil de velocidad de la Sec2 en el plano x-y,ver Figura 3.10.

Fuente: Elaboración propia.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 59

Figura 3.16: Perfil de velocidad de la Sec3 en el plano x-y,ver Figura 3.10.

Fuente: Elaboración propia.

Figura 3.17: Perfil de velocidad de la Sec4 en el plano x-y,ver Figura 3.10.

Fuente: Elaboración propia.

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 60

Figura 3.18: Perfil de velocidad de la Sec5 en el plano x-y,ver Figura 3.10.

Fuente: Elaboración propia

Figura 3.19: Perfil de velocidad de la Sec6 en el plano x-y, Figura 3.10.

Fuente: Elaboración propia

 CAPÍTULO 3. MODELO HIDRÁULICO DEL RÍO MADRE DE DIOS 61

Figura 3.20: Perfil de velocidad de la Sec7 en el plano x-y, ver Figura 3.10.

Fuente: Elaboración propia

3.3.2.2. Tirantes

En el modelo hidráulico se midieron los tirantes en las siete secciones transversales desde
la margen izquierda hacia la margen derecha.(Ver Tabla 3.4)
Se midió la cota con el nivel de anteojo en la superficie libre y el fondo, por diferencia
de cotas se obtuvo el tirante en los seis puntos de la sección transversal medidos desde la
margen izquierda a la margen derecha.

Tabla 3.4: Tirantes en seis puntos en las secciones 1,2 hasta 7 medidos desde la margen
izquierda hacia la margen derecha
Ptos sec1 sec2 sec3 sec4 sec5 sec6 sec7
- y[m] y[m] y[m] y[m] y[m] y[m] y[m]
p1 0.241 0.243 0.224 0.216 0.260 0.258 0.289
p2 0.208 0.211 0.232 0.223 0.227 0.249 0.248
p3 0.205 0.185 0.200 0.194 0.202 0.205 0.222
p4 0.182 0.177 0.171 0.162 0.183 0.149 0.207
p5 0.160 0.159 0.152 0.143 0.171 0.127 0.161
p6 0.151 0.117 0.131 0.122 0.127 0.123 0.121
Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
Capítulo ,9

Sistema Open TELEMAC

MASCARET

Resumen:
En este capítulo se describe el sistema hidroinformático Open TELEMAC MASCARET
y sus componentes como: TELEMAC-2D y SISYPHE. Además, se describen los esque-
mas numéricos y los algoritmos de solución implementados para resolver el sistema de
ecuaciones.

4.1. Descripción del software

El sistema open TELEMAC MASCARET es una herramienta de modelización numérica

que tiene diferentes módulos para analizar el movimiento del agua y sedimento en el do-
minio fluvial, costa, estuario y lacustre. Para realizar una simulación numérica comprende
los siguientes pasos: "pre-procesamiento", para digitalizar los datos y describir el proble-
ma; "simulación", para calcular las propiedades del fluido en el dominio discretizado y
"post-procesamiento", para mostrar y analizar los resultados.
El Laboratoire National d’Hydraulique et Environnement (LNHE) de Electricidad de
Francia (EDF), originalmente desarrolló el programa teniendo en cuenta la comunicación
entre programas.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 63

El programa es ampliamente utilizado, comercializado y desarrollado por SOGREAH Con-

sultants, en estrecha colaboración con EDF-LNHE.
El sistema open TELEMAC MASCARET utiliza principalmente el método de discretiza-
ción de elementos finitos. También incluye el método de volúmenes finitos para el análisis
de flujo bidimensional y transporte de sedimentos no-cohesivos. Todos los módulos son
compilados y ejecutados usando el lenguaje de programación FORTRAN.
A continuación se muestra los módulos que tiene implementado el software Open Tele-
mac Mascaret para realizar una simulación numérica en aguas superficiales, costas, aguas
subterráneas, calidad de agua e hidrológia.

Figura 4.1: Módulos del software Open Telemac Mascaret.

Fuente:Elaboración propia

4.2. Características del sistema open TELEMAC MAS-

CARET

Tal como se ha mencionado previamente, el sistema Open TELEMAC MASCARET se

puede dividir en tres partes distintas: pre-procesamiento, simulación y post-procesamiento.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 64

Cada uno de estas partes implica el uso de diferentes módulos.

El pre-procesamiento implica generar la malla con el uso del módulo (MATISSE) o el
software Blue Kenue.
La simulación comprende en usar el módulo hidrodinámico (TELEMAC-2D y TELEMAC-
3D), transporte de sedimentos (SISYPHE y SEDI-3D), propagación de ondas en áreas
costeras (TOMAWAC), propagación de las olas hacia la costa o hacia los puertos (AR-
TEMIS), calidad del agua (SUBIEF-2D y SUBIEF-3D) y agua subterránea (ESTEL-2D
y ESTEL-3D).
El post-procesamiento incluye RUBENS-2D y POSTEL-3D, son módulos que permiten
visualizar los resultados bidimensionales y tridimensionales respectivamente.
Para este estudio se usó Blue Kenue, para generar la malla no estructurada de elementos
triangulares; TELEMAC-2D, está diseñado para el calculo de las variables hidrodinámicas
en un campo bidimensional; SISYPHE, calcula las variables de transporte de sedimentos y
la evolución del lecho a partir de las variables hidrodinámicas calculados con TELEMAC-
2D y ParaView para vizualización de los resultados.
El programa tiene implementado una biblioteca de elementos finitos BIEF en el sistema
Open TELEMAC MASCARET, para permitir el intercambio de resultados durante la si-
mulación con los diferentes módulos. Esta biblioteca ha sido especialmente diseñada para
varios algoritmos numéricos que resuelve las ecuaciones clásicas de mecánica de fluidos.
Otra librería UTILE es usado para controlar la entrada-salida (Input/Output).

4.3. Entorno computacional

Los principales módulos de simulación están codificados en FORTRAN 90. Por lo tanto,
se puede ejecutar en todos los entornos de trabajo que operan bajo UNIX, WINDOWS.
El software es protegido usando HASP (Hardware Against Software Piracy).
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 65

4.4. Modulo hidrodinámico TELEMAC-2D

El módulo TELEMAC-2D resuelve las cuatro ecuaciones hidrodinámicas, que son pro-
mediadas en la profundidad también conocidas como las ecuaciones de Saint Venant que
derivan de las ecuaciones de Navier-Stokes, sin embargo, necesita ciertos supuestos y apro-
ximaciones para los términos no lineales (TELEMAC-2D, 2001; BIEF, 2004).
Las siguientes ecuaciones son resueltas simultáneamente en una simulación de flujo bidi-
mensional asumiendo la distribución de la presión hidrostática (TELEMAC-2D, 2002).

Ecuación de continuidad:

∂h →
−
+ u. ▽(h) + hdiv(→
−
u)=0 (4.1)
∂t
Ecuación de momento:

∂u ∂Z 1 →
−
+ ugradu = −g + Sx + div(hvt ▽u) (4.2)
∂t ∂x h

∂v → →
− ∂Z 1 →
−
+−
u ▽(u) = −g + Sy + div(hvt ▽v) (4.3)
∂t ∂y h

∂T → →
− ∂Z 1 →
−
+−
u ▽(T ) = −g + div(hvT ▽T ) (4.4)
∂t ∂y h

Dónde:

h (m) profundidad de agua.

u y v (m/s) componentes del vector velocidad U.

T (g/l) trazador pasivo.

g (m/s2 ) aceleración de la gravedad.

vt (m2 /s) viscosidad turbulenta.

vT (m2 /s) coeficiente de difusión de tracer.

 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 66

Z (m) elevación de superficie libre.

t (s) tiempo.

Sh (m/s) fuente y sumidero del fluido.

Sx ,Sy (m2 ) termino fuente y sumidero en las ecuaciones dinámicas.

ST (g/l/s) sumidero y fuente del tracer. h,u,v y T son ecuaciones que serán calculadas.

Sx y Sy m/s2 son los términos fuente que representan al viento, la fuerza de Coriolis, la fuerza
de fricción. Los diferentes términos de las ecuaciones anteriores son resueltas en uno o más
de dos pasos (en el caso de advección por el método de las característica):

Advección de profundidad de agua (h),componentes del vector velocidad (u,v) y traza-

dor pasivo T.

propagación, difusión y fuente de las ecuaciones dinámicas.

Difusión y término fuente de las ecuaciones de transporte.

Este estudio requirió la solución de las ecuaciones anteriores junto con el modelo k −ε, usando
el método de paso fraccional y el método de las características, para calcular las componentes
de velocidad (u y v), la energía cinética turbulenta k y disipación ε. Se considera el esquema
SUPG (Streamlime Upwind Petrov-Galerkin) para la solución del termino advección de la
profundidad h.
Los términos de advección, propagación y difusión en la ecuación (4.2) son resueltos separa-
damente por TELEMAC-2D en cada nodo discretizado entre el tiempo tn y el tiempo tn+1
(denominado paso de tiempo t). Las siguientes ecuaciones son resueltas por TELEMAC-2D
durante el tiempo de paso del término advección:

h − hn
+ ugrad(h) = 0 (4.5)
∆t
u − un
+ ugrad(u) = 0 (4.6)
∆t
v − vn
+ ugrad(v) = 0 (4.7)
∆t
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 67

κ − κn
+ ugrad(κ) = 0 (4.8)
∆t

ε − εn
+ ugrad(ε) = 0 (4.9)
∆t
Las siguientes ecuaciones son resueltas por TELEMAC-2D durante la propagación y la difu-
sión:

hn+1 − h
+ hdiv(u) = 0 (4.10)
△t
un+1 − u ∂Z 1
= −g + div[hνt grad(u)] + FM x (4.11)
△t ∂x h
v n+1 − v ∂Z 1
= −g + div[hνt grad(u)] + FM y (4.12)
△t ∂y h
Las siguientes ecuaciones son resueltas por TELEMAC-2D durante la solución del modelo
κ − ε.
k n+1 − k νt
− div[ grad(k)] = P + Pkv − ε (4.13)
∆t σk
εn+1 − ε νt ε ε2
− div[ grad(ε)] = P C1ε + Pεv − C2ε (4.14)
∆t σε κ κ
Donde u, v, κ y ε son los valores intermedios de u,v k, ε con el mismo paso de tiempo durante
la advección. La viscosidad νt se calcula mediante el modelo de turbulencia utilizando el valor
de la energía cinética y disipación ε.

κ2
ν t = Cµ (4.15)
ε

Los términos: FM x y FM y son los términos fuente. Este estudio ignora toda fuente excepto
la fricción de fondo. Los valores del término fuente debido a la fricción del fondo puede ser
simulado de acuerdo con la Ley de Fricción de Chezy, que originalmente se deriva para un
flujo uniforme. Teniendo en cuenta un dominio plano y la naturaleza no lineal de la ecuación
de Chezy, que son calculados como:

g n+1 p n 2
FM x = u (u ) + (v n )2 (4.16)
hCh2
g n+1 p n 2
FM y = u (u ) + (v n )2 (4.17)
hCh2
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 68

En este estudio se considera el tratamiento numérico adicional para la solución de los términos
no lineales; es decir, hdiv (u). Esto se hace especificando el coeficiente del término implícito
de la velocidad y profundidad. La ecuación (4.18) se resuelve con el mismo tiempo de paso
teniendo en cuenta la profundidad de propagación hprop .

hdiv(u) = hprop div[θu un+1 + (1 − θu )un ] (4.18)

hprop = θh hn+1 + (1 − θh )hn (4.19)

Donde: θu es el valor implicito de la velocidad y θh es el valor implicito de la profundidad.

El término de difusión se trata de manera similar, utilizando los coeficientes para la difusión de
velocidades, el término de difusión es aproximado por la ecuación (4.19) y resuelto definiendo
el valor implicito de la velocidad θd

div[vt grad(u)] = div[vt gradθd un+1 + (1 − θd )un ] (4.20)

Finalmente, el cambio de profundidad h se calcula durante la solución de los términos: −g ∂Z

∂x

y −g ∂Z
∂y
. Esto es aproximado al resolver la ecuación (4.20).

− ggrad(Z) = −gθh grad(δh) − ggrad(hn + Zf ) (4.21)

En cada paso del tiempo, el cambio de la profundidad se agrega a la profundidad inicial para
obtener la superficie libre y continúa para el paso de tiempo sucesivo,(TELEMAC-2D, 2001;
BIEF, 2004)
Cabe señalar que las ecuaciones de Saint-Venant para dos dimensiones supone insignifican-
te velocidad vertical. Esto está relacionado con la condición hidrostática que descuida la
aceleración vertical.

4.5. Modulo de transporte de sedimentos SISYPHE

SISYPHE es el módulo de transporte de sedimentos no cohesivos y evolución del lecho. Utiliza
fórmulas empíricas que pueden expresarse en función de diversos parámetros hidrodinámicos
y de sedimento.
Este módulo no puede ser calculado en el campo de flujo independientemente, por lo tanto,
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 69

las variables hidrodinámicas se imponen en el modelo usando programación en FORTRAN,

que consiste en el intercambio de resultados de las variables hidrodinámicas y transporte de
sedimentos, que es calculado por el acoplamiento de TELEMAC-2D y SISYPHE. La evolución
de lecho es resuelta en cada nodo para calcular el cambio de espesor del sedimento utilizando
la ecuación de balance de masa.

∂Zf
(1 − np ) + div(Qs ) = 0 (4.22)
∂t

Donde np es la porosidad del lecho, Zf es la elevación del lecho y Qs es el transporte de

sedimento por unidad volumen. Las siguientes fórmulas de transporte de sedimentos son
programadas en SISYPHE para el cálculo de transporte Qs (SISYPHE, 2004).

Meyer-Peter & Muller (1948).

Einstein & Brown (1950).

Engelund & Hansen (1967).

Bijker (1968).

Bailard (1981)

Van Rijn (1984)

Dibajnia & Watanabe (1992)

Hunziker (1995)

Soulsby & Van Rijn (1997)

La descripción de estas fórmulas se puede encontrar en referencias relevantes mencionadas

anteriormente: Van Rijn (1993), Yang (1996) y SISYPHE (2004). Por simplicidad, SISYPHE
(2004) considera todas estas fórmulas en función del parámetro Shields θ. En este estudio
se utilizó la fórmula de Meyer-Peter & Müller (1948), Einstein et al. (1950), Van Rijn et al.
(1993).
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 70

4.5.1. Modificación de la magnitud y dirección de la carga del lecho

Tres aspectos claves que deben de considerarse al calcular la magnitud y dirección de la carga
del lecho, (SISYPHE, 2004).

Efecto de la pendiente del lecho.

Efecto del flujo secundario en la dirección del esfuerzo cortante, también conocidos
como flujo helicoidal.

El esfuerzo cortante del lecho se divide en dos componentes: la fricción superficial y la

fuerza de arrastre debido a la forma del lecho.

4.5.1.1. Efecto de la pendiente del lecho

Generalmente, el criterio del movimiento incipiente del sedimento y el transporte de sedimen-

tos es descrito para condiciones de lecho plano. En realidad, el lecho de los ríos son raramente
planos y las partículas normalmente se encuentran con una pendiente.
El efecto de la pendiente es aumentar la tasa de transporte del material en la dirección des-
cendente de la pendiente, y para reducirla cuesta arriba.
Este estudio consideró un factor de corrección tanto para la magnitud como para la dirección
de la tasa de transporte solido antes de resolver la ecuación de evolución del lecho.
El método de corrección se basa en Koch y Flokstra (1981), que es la intensidad de la tasa
de transporte sólido multiplicada por un factor SF para tener en cuenta el cambio en la
coordenada vertical y en la coordenada horizontal a lo largo de la dirección del lecho.

∂Zf
SF = 1 − β (4.23)
∂m

Donde: m es la coordenada en la dirección actual y β es una constante empírica que varía

entre 0.6 y 1.3.
El factor de corrección debido al cambio de dirección de transporte sólido se aplica conside-
rando la dirección del transporte sólido αa y la dirección del esfuerzo cortante de fondo δa
en relación con la dirección del flujo. Si r es la coordenada a lo largo del eje perpendicular al
flujo (transversal o vertical) entonces el cambio en la dirección del transporte sólido se toma
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 71

en cuenta usando la siguiente ecuación:

∂Zf
tanαa = tanδa − β (4.24)
∂r
El valor de la constante empírica β = 1,3 se usa de acuerdo con SISYPHE (2004).
El efecto de la pendiente del lecho es similar al agregar el término de difusión en la ecuación
de evolución del lecho (Ecuación de Exner). Este valor puede tender a suavizar los resultados
y reducir las inestabilidades. Esta formulación; sin embargo, no incluye el efecto del ángulo
de fricción y solo es aplicable para una pendiente suave.
A continuación, en la Figura 4.2 muestra que el material no-cohesivo es afectado por la
pendiente del lecho. El módulo SISYPHE considera efecto de la pendiente en el calculo de
tasa de transporte de sedimentos.

Figura 4.2: Pendiente del lecho en el transporte del material de fondo y suspensión.

Fuente:www. opentelemac. org

4.5.1.2. Corrección del esfuerzo cortante en la dirección del flujo secundario

Debido al efecto del flujo secundario la dirección de transporte de sedimentos no coincide

totalmente con la dirección del esfuerzo cortante.
Engelund (1967) propuso una expresión semi-empírica para incorporar esos efectos en modelos
bidimensionales promediados en la profundidad, se refiere a la corrección transversal de la
vorticidad. El esfuerzo cortante en la dirección del flujo principal puede ser escrito como:
Cf
τ∗ = ρ|u|u (4.25)
2
El esfuerzo cortante del lecho en la dirección transversal puede ser escrito como:

τsecn = τ ∗ tan∆δ (4.26)

 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 72

Con:
v Ah
δ = tan−1 − tan−1 = δ ∗ − ∆δ (4.27)
u rs
El término ∆δ tiene en cuenta el efecto del movimiento espiral sobre las particulas de sedi-
mento. Los ángulos δ ∗ y ∆δ indican respectivamente la dirección del esfuerzo cortante del
lecho (que coincide con la dirección de la velocidad promediada en la profundidad) y la direc-
ción debido al efecto de la corriente secundaria. El cos y sen del ángulo δ ∗ puede ser calculado
por:
τ∗ τy∗
cosδ = x∗
∗ ∗
senδ = ∗ (4.28)
|τ | |τ |
Donde: τx∗ y τy∗ son los componentes del esfuerzo cortante del lecho a partir de las ecuaciones
de momentum relacionados con los campos de velocidad promediados en la profundidad, con
q
módulo |τ | = (τx∗ )2 + (τy∗ )2 . Se han propuesto varias expresiones para el coeficiente espiral
∗

A.
Engelund & Hansen (1967) uso un coeficiente constante A = 7 y (Struiksma et al.) propuso:
√
2ε g
A = 2 (1 − ) (4.29)
k kC

Para los canales curvos, la pendiente de la superficie del agua se puede aproximar por:
′
∂zs α | u |2
g ≈ (4.30)
∂n R
′
Donde: zs es la elevación de la superficie del agua y α un coeficiente que contabiliza si el
′
lecho es liso o rugoso, 0,75 ≤ α ≤ 1.
El radio de curvatura puede expresarse entonces por:
′
α (| u |)2
R= (4.31)
g ∂z
∂n
s

La derivada normal de la elevación de la superficie del agua puede ser calculado como:
u∂zs v∂zs
∂zs ∂y
− ∂x
= n. ▽ zs = (4.32)
∂n |u|

Al reemplazar la ecuación (4.30) en ecuación (4.31):

α′ (| u |)3
R= (4.33)
g( u∂z
∂y
s
− v ∂z
∂x
s
)
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 73

Por ecuación (4.25), el esfuerzo cortante debido al flujo secundario está dado por:

Cf Ah 1
τ sec = ρ ′ g|u| (4.34)
2 α R

Con componentes:
τxsec = τsec v (4.35)

τysec = τsec (−u) (4.36)

Finalmente, las componentes del esfuerzo cortante del lecho que incluye la influencia de las
corrientes secundarias son:

τxsec = τx∗ + τxsec (4.37)

τysec = τy∗ + τysec (4.38)

con:
τ sec τysec
cos(δ) = xsec sen(δ) = sec (4.39)
τ τ
q
y |τ sec | = (τxsec )2 + (τysec )2 el módulo del esfuerzo cortante del lecho, incluye el efecto de la
corriente secundaria.

Figura 4.3: Muestra el flujo secundario helicoidal generado en la vertical debido a la fuerza
centrífuga.

Fuente:www. opentelemac. org

 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 74

4.5.1.3. La correción del flujo secundarío y la pendiente transversal del lecho

debido a la gravedad

La carga total del lecho se puede descomponer en dos componentes x, y.

qb = (qbx , qby ) = (qb cos(α), qb sin(α)) (4.40)

Donde: qb es la tasa de transporte por unidad de ancho, calculada en función de la carga de

sedimento en equilibrio qb0 y α el ángulo de dirección de la carga de fondo del lecho.
Para la carga de sedimentos en equilibrio se puede usualmente escribir como una función del
parametro de Shields (θ − θc ) para un número de Reynolds de una partícula Rp :

qb0 = f (θ − θc ; Rp ) (4.41)

Con esfuerzo cortante en el lecho θ y el parámetro critico de Shields θc . Por ejemplo, la forma
adimensional de Meyer-Peter y Mueller, esta dada por:

qb0 = 8(µθ − θc )3/2 (4.42)

Donde: µ es el coeficiente de rizos.

4.5.1.4. Dirección del transporte de sedimentos

El ángulo α que se forma entre la linea de la dirección del transporte de sedimentos y la

dirección del eje x, que es desviado desde el esfuerzo cortante es combinado por acción de la
pendiente transversal y la corriente secundarias:
1 ∂b
sinδ − f (θ) ∂y
tanα = 1 ∂b
(4.43)
cosδ − f (θ) ∂x

∂b ∂b
Los términos ∂x
y ∂y
representan las pendientes transversales y longitudinales respectiva-
mente. La forma del sedimento es f (θ) una función de ponderación que tiene en cuenta la
influencia de la pendiente transversal del lecho, expresado como función del esfuerzo cortante
adimensional o parámetro de Shields (θ) calculado como:

|τ |
θ= (4.44)
(ρs − ρ)gD50
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 75

Algunas expresiones han sido propuestas en la literatura para f (θ), por ejemplo: Koch
√
(1980)f (θ) = 3θ2
; Talmon et al., f (θ) = β2 θ con 0, 35 ≤ β2 ≤ 1, 5; y Struiksma et al.
√
(1985), f (θ) = 9( Dh )0,3 θ.

4.6. Pre-procesamiento con QGIS y Blue Kenue

El software Blue Kenue ha sido diseñado para generar la malla del dominio computacional y
las condiciones de frontera del área de estudio de una manera interactiva.
Blue Kenue construye una malla no estructurada de elementos triangulares para discretizar
el dominio. Estos elementos se utilizan para la solución de las ecuaciones hidrodinámicas
gobernantes a través del método de elementos finitos. La construcción apropiada de la malla es
de gran importancia porque permitirá entender las complejas interacciones entre la estructura
y la estabilidad numérica. Los siguientes pasos fueron considerados durante la generación de
malla:

Definir el contorno del dominio computacional.

Identificar si hay una isla dentro del dominio.

Definir la densificación local de los nodos.

Estos criterios se pueden definir manualmente o de acuerdo a los elementos externos tales
como velocidad, campos de concentración, etc. Blue Kenue genera dos archivos de salida
esenciales: GEOMETRÍA (extensión. slf) y CONLIM (archivo de condición de frontera, de
extensión. cli). Para construir estos archivos, se debe tener los siguientes datos:

Batimetría en formato .xyz

Dominio en formato i2s.

Modelo digital de terreno (DEM).

Malla con extensión .geo

condiciones de frontera .cli

 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 76

4.6.1. Características de la batimetría

La batimetría describe la elevación del lecho del río, a su vez ayuda a definir la distancia de
los puntos para la densificación de los nodos. Los datos batimétricos se definen con respecto
a un sistema de coordenadas de varias unidades (SI, CGS y Millas Náuticas). Estos datos
deberían estar en el formato .xyz.
El formato .xyz fue adoptado en este estudio.Se generó usando el editor de texto gedit y MS
EXCEL con tres columnas de datos para las coordenada x, coordenada y e coordenada z
respectivamente.

4.6.2. Modelo digital del terreno (DEM)

El Modelo Digital de Terreno (DEM) es el primer intento de triangulación para generar la

malla. El DEM comprende globalmente un mapa de densidad que proporciona un tamaño
óptimo de elementos triangulares equiláteros para garantizar la distribución adecuada de la
malla no estructurada.
La definición de los criterios de malla es una tarea importante en la construcción de la
misma. El criterio indica la dimensión del lado de los elementos triangulares en el sistema
internacional, es decir, metros m. Esto puede imponerse o definirse de muchas maneras. Este
criterio constante global define los triángulos relativamente uniforme para todo el dominio.
Para crear una malla es importante identificar que al incrementar el número de elementos
y reducir el tamaño implica mayor calculo computacional. El refinamiento del mallado es
un factor importante para determinar una solución exacta pero conlleva a mayor calculo
computacional.

4.6.3. Características de la malla

Al generar la malla triangular no estructurada basada en los criterios definidos en el modelo

digital de elevación. La información preliminar sobre la malla se inspecciona considerando
los siguientes aspectos:

La cantidad de nodos.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 77

La cantidad de islas o líneas geométricas.

La cantidad de puntos fijos.

Los valores batimétricos máximos y mínimos.

La dimensión de los segmentos más grandes y más pequeños.

Máxima distorsión de la celda.

Hacer una inspección preliminar de la calidad de la malla, usando el post-procesador Para-

View y Blue Kenue, que muestra las siguientes características en la estructura de la malla:

Batimetría.

Longitud media de los segmentos en los triángulos (relacionados con las distancias entre
nodos).

Área de los triángulos.

Distorsión de los triángulos.

Ángulo más pequeño de los triángulos.

En general, se desea una distribución uniforme de los nodos con triángulos casi equiláteros
de diferentes tamaños. Tal distribución es importante para la estabilidad numérica y para
la reducción del tiempo computacional. La estructura de la malla puede editarse manual-
mente algunos nodos antes de la generación del archivo GEOMETRY (o archivo de malla
de extensión. slf). Durante la modificación manual de la malla, se debe tener cuidado para
evitar la introducción de errores en la malla, especialmente segmentos intersecados y trián-
gulo over-constrained (fuera del rango de los criterios deseados). Cuando la malla cumple los
criterios de especificación, el archivo de GEOMETRÍA que se genera es el formato SERAFIN
(TELEMAC-2D, 2002). Este archivo de geometría es usado en el módulo hidrodinámico y de
transporte de sedimentos del sistema open TELEMAC MASCARET durante la simulación.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 78

4.6.3.1. Malla bidimensional

La malla de TELEMAC-2D consiste en triángulos de tres nodos que se construye utilizando

BIEF MESH, la subrutina FORTRAN del sistema open TELEMAC MASCARET. Cada
elemento tiene numeración local y global en el archivo de geometria en las coordenadas x e
y.
La función base para el análisis de elementos finitos TELEMAC-2D (2001); es obtenida para
los elementos usando números locales y luego se usa una tabla de conectividad para transferir
la numeración local a la numeración global para resolver todo el sistema lineal BIEF (2004).
La figura (4.4) representa solo un elemento triangular con numeración local (orden ascendente
en el sentido de las agujas del reloj).

Figura 4.4: Representación de un elemento triangular 2D y número de nodos.

Fuente:(Hervouet, 2007)

4.6.4. Condiciones de frontera

El último paso en Blue Kenue implica definir las condiciones de contorno. El tipo de límite y
los valores en los diversos nodos del dominio se definen de dos maneras: como una entidad y
como un grupo. Una entidad se refiere a las características de las condiciones de frontera en
un nodo individual, que comprenden el tipo y valor relacionado con las variables desconocidas
(profundidad del agua, componentes de la velocidad). Hay seis tipos diferentes en Blue Kenue,
para elegir y describir la naturaleza de la condición de frontera: adherencia, ola, deslizamiento,
tasa de flujo y valores impuestos. Un grupo de nodos en una linea de contorno, siguen los
siguientes pasos para crear el archivo de condición de frontera:
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 79

Asignar a todos los nodos de la línea de contorno.

Conectar para cada grupo de identidad.

En este estudio, las condiciones del archivo CONLIM del formato SERAFIN (TELEMAC-
2D), fue creado con las condiciones de frontera definida para el acoplamiento del módulo
hidrodinámico y transporte de sedimento.

4.6.4.1. Condición de frontera solida

La condición de frontera sólida se modela como una pared lisa sólida que no produzca ninguna
descarga a través del límite y que ejerce una fricción en la frontera. Cuando esta condición
de frontera se aplica a los nodos, se genera automáticamente un coeficiente generado (ac ) que
es para garantizar el gradiente proporcional a la velocidad local como:

∂u ∂v
= ac u = ac v (4.45)
∂e ∂e

Donde e es la magnitud del vector perpendicular a la pared que se deriva para un perfil de
velocidad logarítmica.

4.6.4.2. Condición de frontera liquida

El flujo y la profundidad del agua están controlados aguas arriba y aguas abajo, respec-
tivamente. La condición de frontera aguas arriba se simula como un límite abierto con un
índice de flujo prescrito, mientras que aguas abajo se simula con un límite abierto con una
profundidad prescrita.
En este caso TELEMAC-2D requiere un perfil de velocidad para inicializar el cálculo de los
componentes de velocidad. El perfil de velocidad se define normal al límite superior y su
norma es aproximada a un valor proporcional de la raíz cuadrada de la profundidad del agua
desde el fondo. Aguas abajo, se mantiene un nivel constante de la superficie libre durante la
simulación.
Una introducción repentina de una descarga puede causar inestabilidad computacional (TELEMAC-
2D), por lo tanto, la descarga se incrementa progresivamente. Se debe crear un archivo de
condición de frontera líquida (hidrograma).
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 80

4.7. Técnicas numéricas usadas por el sistema Open TE-

LEMAC MASCARET
Hay una gran cantidad de técnicas numéricas implementadas en el sistema open TELEMAC
MASACARET. Estas técnicas se utilizan para obtener la solución numérica de ecuaciones
diferenciales. En esta sección se describirá brevemente las técnicas numéricas más utilizadas,
como: el esquema numérico y algoritmos numéricos.

4.7.1. Esquemas numéricos

El sistema open TELEMAC MASCARET, tiene diferentes esquemas numéricos o técnicas

de solución de las ecuaciones del flujo de fluido con la precisión deseada. Los detalles de
aquellos se pueden encontrar en TELEMAC-3D (1998), TELEMAC-2D (2001), TELEMAC-
2D (2002), SISYPHE (2004) y BIEF (2004). En TELEMAC-2D, los esquemas numéricos
comunes son:

Esquema semi-implícito centrado + SUPG.

Esquema conservador + SUPG.

esquema de PSI

Se pueden elegir diferentes esquemas para diferentes variables (como: velocidad, profundidad,
tracer, k − ǫ) durante las simulaciones.
El esquema conservador + SUPG (Streamline Upwind Petrov Galerkin) es válido solo para
la advección de la profundidad. La advección del modelo k − ǫ solo se puede resolverse con el
método de las características. TELEMAC-2D (2002) sugiere usar el método de las caracte-
rísticas para resolver todas las variables. Dependiendo del esquema utilizado, la conservación
de masa se puede mejorar usando las llamadas sub-iteraciones dentro del mismo tiempo de
paso. Cabe señalar que el tiempo de discretización se llama semi-implícito TELEMAC-2D
(2001). Esto implica que el tiempo relacionado con los parámetros de ponderación, como: θh
y θu no deben ser menores que 0.5.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 81

4.7.2. Algoritmo de solución

Los conceptos fundamentales de los algoritmos de solución son similares en TELEMAC-

2D y SISYPHE. Los algoritmos de solución dependen de la consideración de los esquemas
numéricos. Si se elige el método de las características, las ecuaciones son tratadas en dos
pasos usando el método de pasos fraccionados:

Primer paso: Los términos de advección corresponden al transporte de las variables

físicas h,u,v y w. En la alternativa k − ǫ son tratadas teniendo en cuenta el carácter
hiperbólico de las ecuaciones.

Segundo paso: las ecuaciones restantes, es decir: propagación, difusión, términos fuen-
te y los términos de advección cuando no hay un recurso para el método de las caracte-
rísticas. Esta etapa se resuelve mediante la discretización del tiempo por el método de
elementos finitos que permite la eliminación de la no linealidad de las ecuaciones. La
discretización del espacio transforma las ecuaciones continuas en un sistema discreto
lineal donde los valores de h, u, v y w en los nodos son las variables desconocidas. El
sistema es resuelto por el método de gradiente conjugado o GMRES.

En el método de paso fraccionado, con espacio-tiempo discretizado, las incógnitas se encuen-

tran en múltiplos del tiempo de paso t tal que tn = t0 + nt. La derivada con respecto al
tiempo de una función está discretizada como:
∂f f n+1 − f n
= (4.46)
∂t ∆t
Donde: f representa las variables no conocidas h,u,v,w,k o ε. Empieza desde la solución inicial
f 0 en el tiempo t0 . La solución es obtenida por sucesivas iteraciones (tiempos de paso) por
encontrar el primer f 1 y f 2 . Sin el tiempo de paso, este método consiste de encontrar f n+1 ,
empieza desde f n y pasa a un paso intermedio para las variables h,u,v,w,k y ε tal como sigue:
f − fn
+ termino de adveccin = 0 (4.47)
△t
f n+1 − f
+ otros terminos = 0 (4.48)
△t
Donde f es el valor intermedio de f con el mismo tiempo paso, durante la advección. La
aproximación en el tiempo será de orden uno porque ciertos términos de advección no están
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 82

discretizados a tn + △t/2.
Cuando no se considera el método de las características, se resuelve el conjunto completo de
ecuaciones para los términos propagación,difusión y fuente. Esto es realizado por:

Discretización en el tiempo.

Discretización en el espacio.

Hablando en general, la "discretización en el tiempo", llamada semi-implícito. No puede ser

estrictamente aplicado. Esto se debe a que el análisis de Fourier del método semi-implícito
no tiene en cuenta la no linealidad o las condiciones de contorno. Por lo tanto, se propor-
ciona un margen de seguridad en el sistema TELEMAC expresado por una función f como:
f = θt f n+1 + (1 + θt )fn . Con θt > 0,5 permaneciendo siempre cerca de este valor mínimo 0.5
hasta 1.
La "discretización en el espacio", las funciones f se evaluarán solo en los puntos de discreti-
zación (en nodos de la malla) descomponiéndolos en las funciones base como (TELEMAC-2D
2001):
f = Σni=1 fi ψi (4.49)

Donde n es la cantidad de puntos discretos, fi es el valor de la función f en el punto i y ψi ,

es la función asociada con ese punto. Cada una de estas funciones básicas está asociada a
un grado de libertad. Esta función es evalúa en 1 en este punto y en 0 en todos los demás
puntos.
Las funciones de interpolación están dentro de un elemento de referencia, que tiene el elemento
del dominio correspondiente vinculado por transformación iso-paramétrica. Los elementos
de referencia son triángulos de tres nodos en TELEMAC-2D o SISYPHE. Las funciones
base de estos elementos están determinadas usando la transformación iso-paramétrica con
interpolaciones lineales. Esto es porque la transformación iso-paramétrico en triángulos con
interpolación lineal puede invertirse fácilmente.
 CAPÍTULO 4. SISTEMA OPEN TELEMAC MASCARET 83

4.8. Post-Procesamiento
Para el post-procesamiento se usó el software:

Blue Kenue.

ParaVIEW.

Tecplot.
Para el análisis de los resultados y la calibración del modelo.
Capítulo 9

Aplicación de TELEMAC-2D y
SISYPHE en el modelo hidráulico del
Río Madre de Dios - Tramo la Pastora

Resumen:
En este capítulo se describe los pasos que se realizaron para la configuración del modelo hi-
dráulico como: generar la malla, las condiciones de borde, generar archivos de entrada para el
modelo computacional y pasos que se realizó para el acoplamiento del módulo hidrodinámico
(TELEMAC-2D) y transporte de sedimentos morfodinamico (SISYPHE).

5.1. Delimitación del dominio

El dominio se discretizó matemáticamente para formar la malla computacional, es una malla
no estructurada adecuada para el código de elementos finitos que tiene implementado el
sistema Open TELEMAC MASCARET que esta formado por triángulos con nodos en cada
vértice donde las ecuaciones son resueltas durante la simulación para calcular las variables
hidrodinamicas como la velocidad (u,v,w), presion, etc.
Las características del dominio físico utilizado en esta investigación se muestran en la siguiente
Figura (5.1):
 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 85

Figura 5.1: Dimensiones del domino físico.

Fuente: Elaboración propia

5.1.1. Generación de la malla

Para discretizar la zona de estudio, se usó el software Blue Kenue que permite generar la
malla triangular. En el caso del modelo hidráulico del río Madre de Dios se construyó una
malla no estructurada de 0.3 m que tiene extensión .slf. Además es uno de los tres archivos
necesarios para realizar la simulación.

Blue Kenue es una herramienta que es ampliamente usado para generar el mallado. Además
permite gener otro archivo de extensión .cli, donde se define las condiciones de frontera.
Pasos para generar la malla en Blue Kenue:

a) Para generar la malla se cargo 1239922 puntos que tiene la extensión .xyz, abrir el pro-
grama Blue Kenue, click en Open y seleccionar el archivo que tiene formato .xyz.
 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 86

Figura 5.2: Un total de 1239922 puntos vizualizado en el programa Klue Kenue.

Fuente: Elaboración propia

b) Hacer click en CloseLine, dibujar el contorno del dominio. Luego hacer clic derecho y
seleccionar la opción resample para definir la distancia entre puntos a largo de todo el
dominio, en este caso es 0.3m.

Figura 5.3: Se realizó el resample (mismo distanciamiento entre nodos) del dominio

Fuente: Elaboración propia

c) Seleccionar New → T3 Mesh Generator. En la siguiente ventana definir las propiedades

de la malla, donde se tendrá que cambiar el valor de (Default Edge Length) el cual sirve
para definir la separación entre nodos. Aquí se definió la separación de 0.3 m. Click en ok.
 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 87

d) Luego seleccionar File → New → ok.

Figura 5.4: Seleccionar la opción de T3 Mesh Generator

Fuente: Elaboración propia

Figura 5.5: Se selecciono la distancia entre nodos que es de 0.3 m

Fuente: Elaboración propia

e) Arrastrar el archivo de puntos el cual contiene la elevación y el Resample Dominio a

newInterpolator2D.
 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 88

Figura 5.6: Opción interpolacion (2D Interpolator)

Fuente: Elaboración propia

f) Seleccionar New Mesh luego ir a la opción Tools → MapObject → newInterpolator2D →

ok. Luego se realizó la interpolación para que cada punto tenga valores en la coordenada
z.

Figura 5.7: Asignarle elevación a la malla generada por medio de (Map Object)

Fuente: Elaboración propia

CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 89

5.1.2. Generación de los archivos de entrada de TELEMAC-2D

En Blue Kenue se puede generar dos archivos de entrada para realizar la simulación
numérica.

El archivo de geometría que tiene la extensión .slf.

El archivo de condiciones de frontera que tiene una extensión .cli

5.1.2.1. Generación del archivo de geometría .slf

Hacer click en File → New → SELAFIN. Para crear un objeto nuevo, newSelafin, se
genera un Data Items.

Seleccionar la opción BOTTOM, que es el que reconoce open TELEMAC MASCA-

RET como batimetría del relieve de la superficie.

Arrastrar la malla generada que tiene el nombre de BOTTOM a newSelafin.

Luego seleccionar newSelafin para guardar con el nombre geo que tiene la extensión
.slf

Figura 5.8: Archivo SELAFIN Object.

Fuente: Elaboración propia

CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 90

5.1.2.2. Generación del archivo de condiciones de borde .cli

Seleccionar File → Boundary Conditions (Colims), luego elegir la opción BOTTOM.

Luego arrastrar la opción BOTTOM a la vista 2D, para definir la condición de

frontera aguas arriba y abajo.

Hacer doble clic en el nodo y presionar shift para seleccionar la condición de frontera
aguas arriba donde se definió (Open boundary with prescribed Q) y aguas abajo
(Open boundary with prescribed H)

Guardar el archivo como river.cli

Figura 5.9: Archivo Boundary Conditions (Conlim) aguas arriba.

Fuente: Elaboración propia

 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 91

Figura 5.10: Archivo Boundary Conditions (Conlim) aguas abajo.

Fuente: Elaboración propia

Figura 5.11: Condiciones de frontera (Boundary Conditions) inflow Q [m3 /s] y outflow H [m]

Fuente: Elaboración propia

CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 92

5.1.3. Archivos necesarios para ejecutar TELEMAC-2D

Para realizar la simulación numérica con TELEMAC-2D es necesario tener tres archivos:

Archivo de geometría de extensión .slf

Archivo de condiciones de borde. cli

Archivo de extensión .cas

El archivo de extensión .cas se puede descargar de la página web (www.opentelemac.org)

del programa, este archivo permite la ejecución que contiene todos los parámetros como:
esquema numérico, modelo de turbulencia. Este archivo se puede editar con un editor de
texto como Geany, gedit o kate) donde se coloca el nombre de los archivos de geometría,
condición de borde y el nombre del archivo de salida de extensión .sfl.

5.1.3.1. Descripción de los parámetros del archivo .cas

a) Equations, boundary condictions

VELOCITY PROFILES: La opción 4:4 describe que el vector velocidad es nor-

mal al límite y su norma es proporcional a la raíz cuadrada de la profundidad
del agua. Esta opción es válida solo para la tasa de flujo prescrito

PRESCRIBED FLOWRATES: Esto se utiliza para fijar el valor de la tasa de

flujo de una frontera abierta con caudal prescrito.

PRESCRIBED ELEVATIONS: Este opción se útiliza para definir la elevación

de una frontera abierta con la profundidad prescrita.

b) Initial conditions

INITIAL DEPTH: Profundidad inicial.

INITIAL CONDITIONS: es una condición donde se puede iniciar el modelo con

una misma profundidad de agua.

c) Input-output files
CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 93

BOUNDARY CONDITIONS FILE: Colocar el nombre del archivo generado en el

software BlueKenue donde se especificó las condiciones de frontera de extensión
.cli
RESULTS FILE: Colocar el nombre del archivo de salida de extensión .slf
GEOMETRY FILE: Colocar el nombre del archivo geometría de extensión .slf

d ) Secondary currentes
En un canal curvo el flujo experimenta una aceleración radial y la fuerza centrífuga
actúa en proporción a la velocidad media. A su vez, el nivel de agua esta inclinada
radialmente en el banco exterior para producir una súper elevación suficiente para
crear una gradiente de presión que permite equilibrar la fuerza centrífuga. A profun-
didades más superficiales la fuerza centrífuga excede la fuerza de presión por lo que
la fuerza resultante impulsa el fluido hacia afuera. TELEMAC-2D permite tener en
cuenta el efecto de estas corrientes secundarias con los siguientes parámetros.

SECONDARY CURRENTS: activar este parámetro con la palabra YES.

PRODUCTION COEFICIENT FOR SECONDARY CURRENTS: El término
de producción en la ecuación de advección-difusión depende linealmente de este
coeficiente que se puede calibrar.
DISSIPATION COEFFICIENT FOR SECONDARY CURRENTS: EL término
de disipación se puede modificar variando este coeficiente

e) Graphics and listing

VARIABLES FOR GRAPHICS PRINTOUTS: Esto se utiliza para especificar

la lista de variables para las que todos los valores se impriman en cada punto de
la malla y serán almacenados en el archivo de resultados.

f ) Numerical parameters

TIME STEP: Define el tiempo que separa dos instantes de cálculo.

GRAPHIC PRINTOUT PERIOD: Corrige el período para las salidas con el fin
de tener un archivo excesivamente grande.
LISTING PRINTOUT PERIOD: Esto soluciona el período entre ediciones de
time steps. El valor dado es dado en número de time steps.
CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 94

Figura 5.12: Descripción del archivo .cas de TELEMAC2D.

CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 95

Figura 5.13: Uso de los coeficientes de flujo secundario y esquema numérico

Fuente: www. opentelemac. org

5.1.4. Ejecución de TELEMAC-2D

. Para ejecutar TELEMAC-2D de debe tener los siguientes archivos:

Crear una carpeta, en donde contiene los archivos de extensión: geo.slf, river.cli y
pastora.cas.

Ejecutar en la terminal el siguiente comando: telemac2d.py pastora.cas

Se obtiene como output el archivo de resultados de extensión. slf que es utilizado

como un hort − start que es usado para inicializar todas las variables y hacer la
simulación por 24 horas.
CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 96

5.1.5. Configuración del módulo SISYPHE

En cada paso de tiempo, el modelo de transporte de sedimentos incorpora dos pasos.

Primero, se realiza la simulación del módulo hidrodinámico (TELEMAC-2D) para obtener
el archivo de inicialización de las varibles hidrodinamicas, llamado (hort-start).En segundo
lugar se realiza la simulación de forma acoplada de TELEMAC-2D y SISYPHE para
calcular la tasa de transporte apartir de las variables hidrodinámicas en cada iteración en
la conocida ecuación de continuidad de transporte de sedimentos (ecuación de Exner).
Para hacer la simulación numérica de transporte de sedimentos con el módulo SISYPHE
es necesario contar con los siguientes archivos:

Archivo de geometría .slf

Archivo de condiciones de frontera. cli

Archivos de iniciación de variables hort − start.

Archivo de extensión .cas donde especifica los parámetros del módulo hidrodinámico.

Archivo de extensión .cas que contiene los parámetros de transporte de sedimentos

como: diámetro del sedimento, formula de transporte de sedimentos, nombre del
archivo de salida.

Figura 5.14: Archivo de extensión .cas para el módulo morfodinamico.

Fuente: www. opentelemac. org

CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 97

Figura 5.15: Características físicas de la partícula de sedimento.

Fuente: www. opentelemac. org

 CAPÍTULO 5. APLICACIÓN DE TELEMAC-2D Y SISYPHE EN EL MODELO HIDRÁULICO DEL RÍO MADRE DE
DIOS - TRAMO LA PASTORA 98

Figura 5.16: Acoplamiento del módulo hidrodinámico (TELEMAC-2D) y transporte de sedi-

mentos (SISYPHE).

ACOPLAMIENTO DEL MODULO HIDRODINAMICO (TELEMAC-2D) Y

MORFODINAMICO (SISYPHE)

Condiciones
iniciales

Hidrodinámica

Transporte
de
sedimentos Nuevo
tiempo
Evolución
del lecho

Fuente:www. opentelemac. org

Figura 5.17: Algoritmo de solución para el módulo de transporte de sedimentos

Fuente:Elaboración propia
Capítulo 9,

Análisis de resultados de la simulación

numérica

Resumen:
En este capítulo se realiza la validación del modelo numérico con los datos obtenidos en el
ensayo experimental del modelo hidráulico. Se compara la velocidad que se obtuvo en la
simulación con los datos recolectados en el ensayo experimental. Finalmente se compara el
cambio de la morfología del lecho antes y después de la experiencia. Se hace un análisis de
los resultados del modelo numérico tales como: interpretar los efectos del flujo helicoidal,
configuración del proceso de sedimentación y erosión.

6.1. Calibración del módulo hidrodinámico

Los resultados numéricos son analizados en siete secciones transversales (Figura 3.5).
Primero se validó el módulo hidrodinámico TELEMAC-2D donde se usó la corrección de
flujo secundario que tiene implantado el módulo hidrodinámico (Figura 5.12). Se utiliza
esta corrección debido a que en los canales curvos el flujo experimenta una aceleración
radial y una fuerza centrífuga que actúa en proporción a la velocidad media.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 100

6.1.1. Comparación del resultado numérico y experimental de las

velocidades

Se tiene la velocidad media en siete secciones transversales y en cada sección se tiene datos
de velocidad en seis puntos, medidos de la margen izquierda hacia la margen derecha,
cuadro (3.10).
Los resultados del modelo numérico tuvieron que ser interpolados en las coordenadas en
donde fueron medidos los datos de velocidad en el modelo hidráulico.
La Tabla 6.1 muestra los datos de velocidad en seis puntos que fueron medidos desde la
margen derecha hacia la margen izquierda.

(Exp) es la velocidad experimental.

(Nu-sc) es la velocidad experimental usando la corrección de flujo secundario.

(Nu) es la velocidad numérica sin usar la opción de corrección de flujo secundario.

(Error Exp-Nu-sc) es el porcentaje de error entre los datos de la velocidad experi-

mental y numérica considerando la corrección de flujo secundario.

(Error Exp-Nu) es el porcentaje de error entre los datos experimentales y numéricos

sin considerar la corrección de flujo secundario.

Tabla 6.1: Datos de la velocidad en la sección 1 [sec1]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.281 0.271 0.12 3.6 57.3
2 0.376 0.355 0.225 5.6 40.2
3 0.368 0.358 0.22 2.7 40.2
4 0.35 0.355 0.217 1.4 38.0
5 0.339 0.351 0.232 3.5 31.6
6 0.324 0.323 0.206 0.3 36.4
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 101

Figura 6.1: Perfil de velocidad en el plano x-y, se muestra que al considerar la corrección de
flujo secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

Figura 6.2: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 102

Tabla 6.2: Parámetros Estadísticos.

Parámetros Estadísticos. Valores

Coeficiente de correlación múltiple 0.939

2
Coeficiente de determinación R 0.883
R2 ajustado 0.854

Tabla 6.3: Datos de la velocidad en la sección 2 [sec2]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.359 0.368 0.121 2.5 66.3
2 0.462 0.468 0.165 1.3 64.3
3 0.483 0.482 0.183 0.2 62.1
4 0.444 0.451 0.19 1.6 57.2
5 0.423 0.44 0.202 4.0 52.2
6 0.391 0.404 0.209 3.3 46.5
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 6.3: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 103

Figura 6.4: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales.

Fuente: Elaboración propia

Tabla 6.4: Parámetros Estadísticos.

Parámetros Estadísticos. Valores

Coeficiente de correlación múltiple 0.993

2
Coeficiente de determinación R 0.986
R2 ajustado 0.983

Tabla 6.5: Datos de la velocidad en la sección 3 [sec3]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.2 0.196 0.159 2.0 20.5
2 0.408 0.449 0.201 10.0 50.7
3 0.443 0.483 0.204 9.0 54.0
4 0.462 0.466 0.209 0.9 54.8
5 0.432 0.455 0.22 5.3 49.1
6 0.409 0.41 0.234 0.2 42.8
Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 104

Figura 6.5: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

Figura 6.6: Coeficiente de correlación entre los datos obtenidos de la simulacion numérica y
datos experimentales.

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 105

Tabla 6.6: Parámetros Estadísticos.

Parámetros Estadísticos. Valores

Coeficiente de correlación múltiple 0.987

2
Coeficiente de determinación R 0.974
R2 ajustado 0.968

Tabla 6.7: Datos de la velocidad en la sección 4 [sec4]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.19 0.155 0.127 18.4 33.2
2 0.25 0.308 0.157 23.2 37.2
3 0.369 0.483 0.218 30.9 40.9
4 0.394 0.462 0.27 17.3 31.5
5 0.402 0.454 0.288 12.9 28.4
6 0.406 0.395 0.226 2.7 44.3
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 6.7: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 106

Figura 6.8: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales.

Fuente: Elaboración propia

Tabla 6.8: Parámetros Estadísticos.

Parámetros Estadísticos. Valores

Coeficiente de correlación múltiple 0.949

2
Coeficiente de determinación R 0.902
R2 ajustado 0.877

Tabla 6.9: Datos de la velocidad en la sección 5 [sec5]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.21 0.107 0.134 49.0 36.2
2 0.235 0.221 0.193 6.0 17.9
3 0.285 0.387 0.276 35.8 3.2
4 0.336 0.423 0.299 25.9 11.0
5 0.352 0.432 0.341 22.7 3.1
6 0.42 0.419 0.293 0.2 30.2
Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 107

Figura 6.9: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

Figura 6.10: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 108

Tabla 6.10: Parámetros Estadísticos.

Parametros Estadísticos Valores

Coeficiente de correlación múltiple 0.934

2
Coeficiente de determinación R 0.872
R2 ajustado 0.840

Tabla 6.11: Datos de la velocidad en la sección 6 [sec6]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.201 0.123 0.152 38.8 24.4
2 0.282 0.162 0.278 42.6 1.4
3 0.334 0.335 0.353 0.3 -5.7
4 0.339 0.407 0.374 20.1 -10.3
5 0.383 0.463 0.402 20.9 -5.0
6 0.479 0.518 0.215 8.1 55.1
0
Fuente: Proyecto N 358-PNICP-PIAP-2014, auspiciado por INNOVATE.

Figura 6.11: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 109

Figura 6.12: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales.

Fuente: Elaboración propia

Tabla 6.12: Parámetros Estadísticos.

Parámetros Estadísticos. Valores

Coeficiente de correlación múltiple 0.922

2
Coeficiente de determinación R 0.851
R2 ajustado 0.814

Tabla 6.13: Datos de la velocidad en la sección 7 [sec7]

Velocidad en el plano x-y

Puntos Exp Nu-sc Nu Error Error
[m/s] [m/s] [m/s] Exp Nu-sc Exp Nu
1 0.202 0.119 0.156 41.1 22.8
2 0.263 0.202 0.189 23.2 28.1
3 0.309 0.413 0.25 33.7 19.1
4 0.362 0.437 0.268 20.7 26.0
5 0.376 0.432 0.303 14.9 19.4
6 0.407 0.388 0.233 4.7 42.8
Fuente: Proyecto N 0 358-PNICP-PIAP-2014, auspiciado por INNOVATE.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 110

Figura 6.13: Perfil de velocidad en el plano x-y, se muestra al considerar la corrección de flujo
secundario se aproxima a los datos experimentales.

Fuente: Elaboración propia

Figura 6.14: Coeficiente de correlación entre los datos obtenidos de la simulación numérica y
datos experimentales.

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 111

Tabla 6.14: Parámetros Estadísticos.

Parámetros Estadísticos Valores

Coeficiente de correlación múltiple 0.885

2
Coeficiente de determinación R 0.784
R2 ajustado 0.730

Figura 6.15: Se muestra las velocidades máximas en el tramo curvo, donde la velocidad está
en el rango de 0-0.5m/s

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 112

6.2. Calibración del módulo morfodinamico

Del acoplamiento del módulo hidrodinámico y transporte de sedimentos (TELEMAC-2D

y SISYPHE), se consideró la corrección de flujo secundario, el efecto de la pendiente del
lecho y la formula de transporte de sedimentos (Meyer-Peter y Müller) para tener una
mejor aproximación a los resultados como la tasa de transporte de sedimentos, evolución
del lecho. Por medio de la ecuación de continuidad de transporte de sedimentos (Ecuación
de Exner). Al considerar esos parámetros se muestra una marcada erosión en el banco
exterior y sedimentación en el banco interior.
Se tiene cuatro secciones transversales Figura (3.5), donde se compara con la batimetría
antes y después del experimento en el modelo hidráulico del río Madre de Dios y con los
resultados obtenidos del modelo numérico. Se observa el cambio de la morfología del lecho
debido al efecto de la pendiente del lecho y a la ocurrencia del flujo secundario en el tramo
curvo.

(expe1) es la batimetría del lecho antes del ensayo experimental.

(exp2) es la batimetría del lecho después del ensayo experimental en modelo hidráu-
lico.

(nume) resultado numérico obtenido por el acoplamiento del módulo (TELEMAC-2D

y SISYPHE).
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 113

Figura 6.16: Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S2)
SECCIÓN TRANSVERSAL SEC2

-2.6

-2.7

-2.8
Profundidad(m)

-2.9

-3

-3.1

-3.2 expe1
exp2
nume
-3.3
-6 -5 -4 -3 -2 -1 0 1 2 3 4

d(m)

Fuente: Elaboración propia

Figura 6.17: Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S4)
SECCIÓN TRANSVERSAL SEC4

-5.4

-5.6

-5.8
Profundidad (m)

-6

-6.2

-6.4 expe1
exp2
nume
-6.6
-5 -4 -3 -2 -1 0 1 2 3 4

d(m)

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 114

Figura 6.18: Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Sección S6)
SECCIÓN TRANVERSAL SEC6

-7.3

-7.4

-7.5

-7.6
Profundidad (m)

-7.7

-7.8

-7.9

expe1
-8 exp2
nume
-8.1
-7 -6 -5 -4 -3 -2 -1 0 1 2

d(m)

Fuente: Elaboración propia

Figura 6.19: Comparación de secciones transversales con la batimetría antes y después del
experimento y el resultado del modelo numérico. (Seccción S8)
SECCIÓN TRANSVERSAL SEC8

-1.1

-1.2

-1.3

-1.4
Profundidad(m)

-1.5

-1.6

-1.7

expe1
-1.8 exp2
nume
-1.9
-6 -5 -4 -3 -2 -1 0 1 2 3 4

d(m)

Fuente: Elaboración propia

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 115

6.3. Análisis de sensibilidad del parámetro de flujo se-

cundario

Se realizo la simulación bidimensional del modelo hidráulico del río Madre de Dios, de-
bido a que un modelo bidimensional de profundidad promediada no puede representar la
velocidad en la vertical en tramos curvos porque se genera un movimiento helicoidal que
es altamente tridimensional, para ello se uso la opción de corrección de flujo secundario
en el módulo de transporte de sedimentos, SECONDARY CURRENTS ALPHA COEF-
FICIENT α = 0,7.
Además un análisis de sensibilidad para el parámetro alpha para valores de 0.7, 0.9, 1.2.
El valor que representa mejor el proceso de erosión y sedimentación es el valor de SE-
CONDARY CURRENTS ALPHA COEFFICIENT α = 0,7.

Figura 6.20: Coeficiente de flujo secundario α=0.7

Fuente:Elaboración propia.
CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 116

Figura 6.21: Coeficiente de flujo secundario α=0.9

Fuente:Elaboración propia.

Figura 6.22: Coeficiente de flujo secundario α=1.2

Fuente:Elaboración propia.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA VI 117

Figura 6.23: Para el coeficiente de flujo secundario α = 0.7, se muestra una marcada erosión en la margen derecha y sedimentación en la margen izquierda.

Coeficiente de Flujo Secundario α=0.7

Sedimentación

Erosión

Fuente: Elaboración propia.

 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 118

6.4. Discusiones y comentarios

6.4.1. Corrección por la corriente secundaria

El análisis experimental se realizó en el Laboratorio Nacional de Hidráulica – LNH, al rea-

lizar un modelo numérico fue necesario hacer un análisis de sensibilidad de los parámetros
que son los siguientes.

6.4.1.1. El coeficiente de producción para la corriente secundaria As

Determina la tasa de producción de vorticidad. La evidencia experimental sugiere que

este coeficiente debe ser mayor que 5. Este coeficiente está asociado con la advección de
vorticidad. En el presente caso de estudio se usó el valor de 7 que llego a representar mejor
el problema físico como se muestra en el capítulo 6. Además se calculó el error porcentual
donde muestra que al usar ambos coeficientes el error porcentual es mínimo con respecto
a los datos medidos en la etapa experimental.

6.4.1.2. El coeficiente de disipación para la corriente secundaria Ads

Este coeficiente se establece en general en 0,5. Este coeficiente está asociado con la difusión
lateral de la vorticidad fluvial debido a la turbulencia. En el presente estudio se usó dicho
coeficiente con el valor de 0.5 y la difusión lateral (VELOCITY DIFFUSIVITY = 1.E-6)
de la vorticidad con el modelo de turbulencia K-epsilon como se muestra en la Figura
5.12. Ambos parámetros fueron sugeridos por Bernard y Schneider (1992) que los valores
de ambos coeficientes empíricos deben establecerse mediante la comparación de prueba y
error de las velocidades predichas y medidas. Tal como se desarrolló en el presente trabajo.

6.4.2. Corrección por efecto de pendiente (Effect Slope)

A continuación se muestra la fórmula que tiene implementada el software Open Telemac

Mascaret para usar la opción de corrección por efecto de pendiente. FÓRMULA PARA
EL EFECTO DE LA PENDIENTE = 1 es la formula de Koch y Flokstra donde aumenta
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA 119

dz
la velocidad de transporte cuesta abajo y la reduce por el factor (1 − B dx ), en lugar de
cambiar el esfuerzo cortante crítico. La palabra clave Beta (B) tiene un valor predetermi-
nado de 1.3. Debido a que la batimetría del modelo hidráulico tiene una alta resolución y
que el problema a estudiar es altamente tridimensional se usó la opción 1 de corrección de
pendiente. Finalmente al usar la corrección de flujo secundario, efecto de pendiente del le-
cho y el cálculo acoplado de las variables hidrodinámicas (velocidad, tirante) y transporte
de sedimentos (erosión, sedimentación y evolución del lecho) se obtuvo mejores resultados
como se muestra en el ítem 6.2. Donde se muestra resultados de variación del lecho con
la corrección de flujo secundario (rojo) y sin corrección de flujo secundario (verde) y bati-
metría inicial (verde). En la figura 6.24 se muestra los resultados de simulación numérica
después de acoplamiento de la hidrodinámica y transporte de sedimentos obteniendo una
tasa de erosión en el banco exterior. Es importante considerar ambos parámetros en un
software bidimensional que resuelven las ecuaciones Saint Venat en aguas poco profundas.
Finalmente los resultados presentados son los que se plantearon en los objetivos.
Conclusiones

Se simulo numéricamente las características del flujo helicoidal que se forma en la

sección transversal del tramo curvo del modelo hidráulico del rio Madre de Dios que
experimenta una aceleración radial y la fuerza centrífuga que actúa en proporción a
la velocidad media. A su vez la superficie del agua esta inclinada radialmente en el
banco exterior para producir súper-elevación que genera un gradiente de presión para
equilibrar el promedio de la fuerza centrífuga. Además por los resultados de la simula-
ción permite concluir que el perfil de velocidades en el plano x-y debido a la formación
del flujo helicoidal las velocidades son mayores en la margen derecha llegando a ser el
máximo de 0.5m/s, como se muestra en la Figura 6.15.

Se observó el efecto del flujo helicoidal en las riberas del tramo curvo, debido a la
dirección de las líneas de corriente cerca del lecho y de la superficie libre que son
opuestas y generan recirculación cerca a las riberas que genera el colapso.
Cabe resaltar que al usar la corrección de flujo secundario que son los coeficientes As
Producción de Flujo Secundario (Production for secondary currents) que su valor es
7.071 y el coeficiente Ads Coeficiente de Disipación para Flujo secundario (Dissipation
coefficient for secondary flow) que su valor es 0.5 para obtener mejores resultados en
una simulación bidimensional fue presindible usar ambos coeficientes porque el flujo
en el tramo de estudio es tridimensional. Se obtuvieron mejores resultados como se
muestra en las Figuras 6.1, 6.3 el perfil de velocidades obtenidos de la simulación se
aproximan a las velocidades medidas en la etapa experimental.

Se analizó el comportamiento de flujo helicoidal que ocurre en tramo curvo, para ello
se incorpora la corrección de flujo secundario y se consiguió una adecuada representa-
 

ción del campo de velocidades, disminuyendo los errores de manera significativa. Por
lo tanto, es necesario incorporar la corrección de flujo secundario para mantener la
no-uniformidad del campo de velocidades en la vertical en modelos de profundidad
promediada.

Se interpretó e incorporo el efecto del flujo helicoidal, en el módulo de transporte de

sedimentos (SISYPHE), para analizar la dinámica del transporte de sedimento. La tasa
de transporte y la erosión incrementada debido al efecto del flujo helicoidal cerca a las
riberas y en el lecho.
Cabe resaltar que al usar la corrección de flujo secundario en un modelo bidimensional
al realizar el acoplamiento de módulo hidrodinámico y transporte de sedimentos se
obtuvieron mejores resultados de la erosión y sedimentación como se muestra en las
Figuras 6.18 y 6.19 donde los resultados de la simulación numérica de cambio de lecho
se aproxima a los resultados obtenidos en la etapa experimental donde el levantamiento
topografíco de la batimetría se realizo con un Escáner Láser de alta presición. Además
en la Figura 6.20 se muestra los resultados de la erosión originada por el efecto del flujo
secundario.

Se configuro el modelo numérico para hacer la simulación hidrodinámica y transporte

de sedimentos con el módulo TELEMAC-2D y SISYPHE para interpretar el efecto del
flujo helicoidal que tiene un efecto directo en los mecanismos de transporte de erosión
y sedimentación. Como se muestra en la Figura 6.21 y 6.24 el cambio de la morfología
del lecho debido al proceso de erosión en las margenes de las riberas.
 

Recomendaciones
Se realizó simulaciones con el módulo hidrodinámico TELEMAC-2D; donde es difícil repre-
sentar parámetros tridimensionales debido a que el flujo en canales curvos induce a una
fuerza centrífuga que genera corriente secundaria y súper-elevación de la superficie de agua.
Se recomienda realizar una simulación tridimensional. Con una parametri-zación relevante
de efectos tridimensionales que pueden ser incorporados en un modelo hidrodinámico de pro-
fundidad promediada (bidimensional) en este casos no pueden capturar algunos parámetros
tridimensionales. Estas limitaciones pueden ser representadas en un modelo tridimensional.
Una simulación numérica tridimensional con TELEMAC-3D, donde se puede obtener una
solución numérica de las ecuaciones de Navier Stokes el cual representara mejor los paráme-
tros hidromorfodinamicos. Actualmente, la velocidad de cómputo y los programas de pre y
post-procesamiento hacen posible la simulación de casos reales con una alta resolución, que
es posible estudiar fenómenos altamente complejos
ANEXOS
 124

Fotografias del modelo hidráulico del Río Madre de

Dios -Tramo la Pastora

Figura 6.24: Se empleo partículas de tecnopor para visualizar la dirección del flujo

Fuente:Elaboración propia.

Figura 6.25: Configuración del lecho despues del experimento, se muestra el cambio de la
morfologia del lecho.

Fuente:Elaboración propia.
 125

Figura 6.26: Configuración del lecho despues del experimento, se muestra el cambio de la
morfologia del lecho.

Fuente:Elaboración propia.

Fotografias del resultado del modelo numérico

Figura 6.27: Abrir los resultados en de la simulación en el software Blue Kenue.

Fuente:Elaboración propia.
 126

Figura 6.28: Hacer clic derecho en el espacio trabajo para quitar la grilla del espacio de
trabajo.

Fuente:Elaboración propia.

Figura 6.29: Se muestran las isolineas en el modelo experimental, donde el color rojo es la
máxima velocidad.

Fuente:Elaboración propia
 127

Figura 6.30: Se usó fluoresceína para ver el comportamiento del fluido en el modelo hidráulico, y se muestra que tiene el mismo patrón en los resultados del
modelo numérico.

Fuente: Elaboración propia.

 128

Datos de velocidad medidos en el modelo hidráulico del

río Madre de Dios - Tramo la Pastora

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
210 362.71 6.07 1.78 362.76
209 364.95 3.85 1.11 364.97
208 366.45 4.91 1.70 366.49
207 365.86 4.97 1.84 365.90
206 365.95 6.59 2.27 366.02
205 366.77 6.58 2.23 366.83
204 366.95 6.89 2.51 367.02
203 366.43 7.61 2.72 366.52
202 365.62 8.24 3.02 365.72
201 365.63 8.50 3.28 365.74
200 365.78 7.67 3.44 365.87
199 365.43 9.09 3.79 365.56
198 364.60 8.43 3.99 364.72
197 364.59 9.29 4.27 364.73
196 363.93 9.99 4.38 364.10
195 363.69 9.70 4.49 363.84
194 363.37 9.80 4.73 363.53
193 362.23 10.11 5.05 362.40
192 361.46 10.05 5.19 361.64
191 360.71 10.43 5.50 360.90
190 360.46 10.53 5.67 360.65

Tabla 6.15: Tabla de velocidades de la sección 1 en el punto 1

Fuente:Estudio De Línea de Base – FINCyT.
 129

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
189 360.62 11.49 5.76 360.85
188 361.56 11.65 6.05 361.80
187 359.67 11.77 5.75 359.91
186 359.34 12.05 6.28 359.60
185 360.32 12.01 6.49 360.58
184 360.58 12.07 6.98 360.85
183 359.44 12.84 7.10 359.74
182 359.69 12.69 7.24 359.98
181 358.86 13.27 7.29 359.18
180 360.17 13.93 7.09 360.51
179 359.79 14.54 7.00 360.15
178 358.04 14.59 6.73 358.40
177 357.30 15.37 7.40 357.71
176 357.45 15.85 7.01 357.87
175 347.26 10.84 -1.39 347.44
174 351.31 10.01 -1.32 351.45
173 350.81 11.52 -1.31 351.00
172 352.05 10.35 -1.18 352.20
171 352.79 10.66 -0.82 352.95
170 352.49 11.13 -0.95 352.67

Tabla 6.16: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 130

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
169 353.68 12.26 -1.05 353.89
168 353.43 12.01 -1.13 353.64
167 352.54 10.98 -1.33 352.71
166 352.64 10.97 -1.39 352.81
165 353.13 10.53 -1.49 353.29
164 353.60 10.85 -1.65 353.77
163 354.36 10.25 -1.50 354.51
162 354.49 10.18 -1.30 354.64
161 355.36 10.11 -1.37 355.50
160 354.98 10.00 -1.26 355.12
159 355.35 9.65 -1.18 355.49
158 354.75 9.35 -0.88 354.88
157 354.15 9.12 -0.79 354.27
156 354.73 8.84 -0.61 354.84
155 353.44 9.28 -0.47 353.56
154 353.24 9.42 -0.42 353.37
153 352.45 9.79 -0.63 352.59
152 354.37 10.08 -0.59 354.51
151 351.82 10.47 -0.43 351.97
150 352.03 10.96 -0.38 352.21
149 354.68 11.59 0.19 354.87

Tabla 6.17: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA VI 131

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
148 354.46 11.61 0.42 354.65
147 353.52 12.06 0.52 353.73
146 355.80 13.39 0.76 356.05
145 356.00 12.91 0.99 356.23
144 352.88 14.75 0.77 353.19
143 352.84 14.75 0.82 353.15
142 354.36 15.05 1.18 354.68
141 353.11 14.48 1.23 353.41
140 355.55 25.52 -9.55 356.60
139 356.50 6.55 -9.45 356.68
138 355.87 17.34 -9.91 356.43
137 355.38 20.60 -9.54 356.10
136 354.85 14.04 -8.54 355.23
135 355.29 18.39 -9.26 355.89
134 356.57 17.61 -9.42 357.13
133 355.14 18.30 -9.05 355.72
132 353.71 16.13 -9.41 354.20
131 353.45 22.06 -9.30 354.26
130 352.90 21.75 -9.39 353.69
129 353.61 23.97 -9.12 354.54
128 353.54 23.65 -9.45 354.46
127 351.40 24.73 -8.12 352.37

Tabla 6.18: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 CAPÍTULO 6. ANÁLISIS DE RESULTADOS DE LA SIMULACIÓN NUMÉRICA VI 132

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
126 351.16 23.02 -8.46 352.01
125 350.41 23.81 -8.18 351.31
124 349.53 23.30 -7.92 350.40
123 349.11 23.44 -7.48 349.97
122 346.36 22.49 -7.35 347.17
121 346.14 22.52 -6.62 346.94
120 345.92 20.82 -6.39 346.61
119 345.22 32.91 -7.61 346.87
118 343.52 42.27 -9.22 346.23
117 345.31 49.00 -10.05 348.92
116 344.06 17.42 -5.83 344.55
115 343.50 32.65 -7.81 345.13
114 341.36 36.40 -7.16 343.37
113 342.93 16.66 -3.97 343.36
112 343.65 30.42 -5.76 345.04
111 340.74 34.89 -5.08 342.56
110 339.40 22.31 -5.80 340.18
109 333.90 41.76 -4.10 336.53
108 336.00 14.53 -9.29 336.44
107 334.47 38.80 -4.99 336.76
106 338.56 21.67 -5.81 339.30

Tabla 6.19: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 133

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
105 339.40 29.36 -1.06 340.67
104 338.29 31.25 -0.93 339.73
103 339.30 32.45 -0.84 340.85
102 339.31 33.63 -0.29 340.97
101 338.97 33.57 -0.05 340.63
100 340.14 35.47 -0.03 341.99
99 339.93 35.89 0.19 341.82
98 338.69 35.52 0.40 340.55
97 338.36 36.53 0.51 340.32
96 337.50 36.71 1.03 339.49
95 337.69 36.47 1.21 339.66
94 338.29 36.82 1.45 340.29
93 338.38 37.10 1.58 340.41
92 337.48 37.33 1.77 339.55
91 337.39 38.90 2.19 339.63
90 336.96 39.53 2.32 339.28
89 335.72 39.95 2.60 338.09
88 334.01 41.01 2.46 336.53
87 335.14 40.51 2.60 337.59
86 331.91 39.47 2.61 334.26
85 331.29 40.61 2.36 333.78

Tabla 6.20: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 134

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
84 332.28 41.17 2.05 334.82
83 330.98 42.39 1.88 333.69
82 328.44 44.51 1.82 331.44
81 328.50 44.27 1.99 331.47
80 329.71 45.11 2.26 332.79
79 329.33 45.87 2.84 332.52
78 329.99 46.62 3.25 333.29
77 329.74 46.82 3.72 333.07
76 328.18 47.98 3.57 331.69
75 323.29 47.81 3.49 326.82
74 321.73 47.84 3.53 325.28
73 318.62 47.60 3.46 322.18
72 313.41 47.25 3.54 316.97
71 306.16 47.25 3.91 309.81
70 311.86 20.31 3.17 312.53
69 323.68 22.23 3.92 324.47
68 330.36 22.19 2.70 331.12
67 340.19 22.85 3.09 340.97
66 342.02 23.35 3.61 342.83
65 344.92 22.88 3.83 345.70
64 346.27 23.24 3.79 347.07
63 346.14 22.06 4.28 346.87

Tabla 6.21: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 135

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
62 345.58 22.48 4.60 346.34
61 345.89 23.75 4.75 346.74
60 346.07 24.36 4.82 346.96
59 345.73 25.82 5.09 346.73
58 346.06 25.05 5.55 347.01
57 345.62 26.21 5.69 346.66
56 345.02 26.23 5.93 346.07
55 344.16 26.55 6.13 345.24
54 343.20 26.87 6.39 344.31
53 343.10 28.16 6.64 344.32
52 341.45 27.94 6.80 342.66
51 340.38 27.51 7.15 341.56
50 339.13 28.05 7.31 340.36
49 338.97 27.63 7.54 340.18
48 339.65 26.62 7.48 340.78
47 338.70 27.35 7.64 339.88
46 338.48 25.66 7.51 339.53
45 336.85 25.70 7.67 337.91
44 335.93 25.22 7.91 336.97
43 335.87 24.13 8.38 336.84
42 334.38 24.80 8.43 335.41

Tabla 6.22: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
 136

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
41 335.08 23.69 8.70 336.03
40 335.18 25.21 8.71 336.24
39 333.60 25.25 8.63 334.66
38 331.96 25.76 8.69 333.07
37 330.68 25.37 8.65 331.77
36 329.05 24.38 9.25 330.08
35 234.63 41.22 3.45 238.25
34 261.43 42.16 3.57 264.83
33 289.49 49.75 1.72 293.74
32 294.32 48.33 2.10 298.27
31 295.17 48.91 2.26 299.21
30 295.60 48.36 2.74 299.54
29 294.30 48.73 3.23 298.33
28 290.07 49.52 3.81 294.30
27 287.40 47.73 4.49 291.37
26 283.39 46.37 5.26 287.21
25 280.23 47.40 6.06 284.27
24 277.42 46.36 6.74 281.35
23 273.49 46.19 7.51 277.47
22 269.31 45.38 8.20 273.23
21 266.32 44.89 9.28 270.24
20 262.45 43.80 10.35 266.28

Tabla 6.23: Tabla de velocidades de la sección 1 en el punto 1

Fuente:Estudio De Línea de Base – FINCyT.
 137

Tirante Vx Vy Vz Mag. V
[mm] [mm/s] [mm/s] [mm/s] [mm/s]
19 258.78 41.48 11.43 262.33
18 253.15 39.63 12.59 256.54
17 247.90 38.41 13.60 251.23
16 240.51 36.51 14.17 243.68
15 231.85 33.75 14.97 234.77
14 222.44 31.23 15.08 225.13
13 212.29 29.28 15.56 214.87
12 201.30 27.44 15.90 203.78
11 172.05 16.61 13.56 173.39
10 63.41 -8.94 -1.57 64.06
9 15.48 -11.59 0.27 19.34
8 4.16 -0.06 0.56 4.20
7 0.33 0.02 -0.01 0.33
6 0.26 0.14 -0.03 0.30
5 2.60 0.31 0.09 2.62
4 0.37 1.26 -0.05 1.31
3 0.05 0.43 0.01 0.44
2 -0.41 1.68 -0.24 1.75
1 -2.84 0.41 -0.07 2.87

Tabla 6.24: Tabla de velocidades de la sección 1 en el punto 1

.
Fuente:Estudio De Línea de Base – FINCyT.
Bibliografía

Anderson, J. D. & Wendt, J. (1995). Computational fluid dynamics, vol. 206. Sprin-
ger.

Armitage, N. & McGahey, C. (2003). A unit stream power model for the prediction
of local scour in rivers. Water Research Commission. WRC Report (1098/1), 03.

Bagnold, R. A. (1966). An approach to the sediment transport problem from general

physics. US government printing office.

Brown, C. B. (1950). Sediment transportation.

Chanson, H. (1999). The hydraulics of open channel flow: An introduction. physical

modelling of hydraulics .

Cheng, N.-S. (1997). Simplified settling velocity formula for sediment particle. Journal
of hydraulic engineering 123(2), 149–152.

Chien, N. & Wan, Z. (1999). Mechanics of sediment transport. American Society of

Civil Engineers.

Colebrook, C. & White, C. (1937). Experiments with fluid friction in roughened

pipes. Proc. R. Soc. Lond. A 161(906), 367–381.

Einstein, H. A. et al. (1950). The bed-load function for sediment transportation in open
channel flows, vol. 1026. Citeseer.

Engelund, F. & Hansen, E. (1967). A monograph on sediment transport in alluvial

streams. Technical University of Denmark 0stervoldgade 10, Copenhagen K. .
 BIBLIOGRAFÍA 139

Gessler, J. (1970). Self-stabilizing tendencies of alluvial channels. Journal of the Wa-

terways, Harbors and Coastal Engineering Division 96(2), 235–249.

Graf, W. H. & Altinakar, M. S. (1998). Fluvial hydraulics: Flow and transport

processes in channels of simple geometry. 551.483 G7.

Hervouet, J.-M. (2007). Hydrodynamics of free surface flows: modelling with the finite
element method. John Wiley & Sons.

Hrenya, C., Miller, S., Mallo, T. & Sinclair, J. (1998). Comparison of low
reynolds number k-ε turbulence models in predicting heat transfer rates for pipe flow.
International journal of heat and mass transfer 41(11), 1543–1547.

Iyer, G. R. & Yavuzkurt, S. (1999). Comparison of low reynolds number k–ε models
in simulation of momentum and heat transport under high free stream turbulence.
International journal of heat and mass transfer 42(4), 723–737.

Julien, D. J., Pierre y Anthony (2002). Bed load motion and grain sorting in a
meandering stream. Journal of Hydraulic Research 40(2), 125–133.

Koch, F. C., FG (1980). Bed Level Computations for Curved Alluvial Channels: Prepa-
red for the 19th IAHR Congress, New Delhi, India, February 1981. Waterloopkundig
Laboratorium.

Launder, B. E. & Spalding, D. B. (1983). The numerical computation of turbulent

flows. In: Numerical Prediction of Flow, Heat Transfer, Turbulence and Combustion.
Elsevier, pp. 96–116.

Lohnes, R., RA y Handy (1968). Slope angles in friable loess. The Journal of Geology
76(3), 247–258.

Melville, B. W. & Raudkivi, A. J. (1996). Effects of foundation geometry on bridge

pier scour. Journal of Hydraulic Engineering 122(4), 203–209.

Meyer-Peter, E. & Müller, R. (1948). Formulas for bed-load transport. In: IAHSR
2nd meeting, Stockholm, appendix 2. IAHR.
 BIBLIOGRAFÍA 140

Morell, M., Tassi, P. & Vionnet, C. (2014). Flow pattern at a river diffluence at the
alluvial system of the paraná river. Journal of Applied Water Engineering and Research
2(2), 140–156.

Pope, S. B. (2001). Turbulent flows.

Puhakka, M., Kalliola, R., Rajasilta, M. & Salo, J. (1992). River types, si-
te evolution and successional vegetation patterns in peruvian amazonia. Journal of
Biogeography , 651–665.

Rastogi, A. K. & Rodi, W. (1978). Predictions of heat and mass transfer in open
channels. Journal of the Hydraulics division 104(3), 397–420.

Raudkivi, A. & Sutherland, A. (1981). Scour at bridge crossings: Wellington, new

zealand, national roads board, road research unit. Bulletin 54, 100.

Riesterer, J., Wenka, T., Brudy-Zippelius, T. & Nestmann, F. (2016). Bed

load transport modeling of a secondary flow influenced curved channel with 2d and 3d
numerical models. Journal of Applied Water Engineering and Research 4(1), 54–66.

Schlichting, H., Gersten, K., Krause, E., Oertel, H. & Mayes, K. (1955).
Boundary-layer theory, vol. 7. Springer.

Shields, A. (1936). Application of similarity principles and turbulence research to bed-

load movement .

Simon, A. (1989). Shear-strength determination and stream-bank instability in loess-

derived alluvium, west tennessee, usa. Applied Quaternary Research, EJ DeMulder
and B. l Hageman (Editors). AA Balkema Publishers, Rotterdam, The Netherlands ,
129–146.

Simons, D. B. & Richardson, E. V. (1966). Resistance to flow in alluvial channels.

US Government Printing Office.

Struiksma, N., Olesen, K., Flokstra, C. & De Vriend, H. (1985). Bed deforma-
tion in curved alluvial channels. Journal of Hydraulic Research 23(1), 57–79.
 BIBLIOGRAFÍA 141

Talapatra, S. & Ghosh, S. (1983). Incipient motion criteria for flow over a mobile
bed sill. In: Proceedings of the 2nd International Symposium on River Sedimentation.
Water Resources and Electric Power Press Nanjing.

Tassi, K. E. k. A., D. Wang (2014). Flow and bed morphology in large-amplitude

meanders, comparison of two and three-dimensional numerical models. Water Resources
Research 45(2).

Terzaghi, K., Peck, R. B. & Mesri, G. (1996). Soil mechanics in engineering practice.
John Wiley & Sons.

Van Rijn, L. C. et al. (1993). Principles of sediment transport in rivers, estuaries and
coastal seas, vol. 1006. Aqua publications Amsterdam.

Vanoni, V. A. & Nomicos, G. N. (1959). Resistance properties of sediment-laden

streams. Journal of the Hydraulics Division 85(HY5), 77–107.

Versteeg, W., HK Malalasekera (1995). ,“an introduction to computational fluid

dynamics: The finite volume method” addison wesley longman, ltd., harlow, england. .

Villaret, C., Hervouet, J.-M., Kopmann, R., Merkel, U. & Davies, A. G.

(2013). Morphodynamic modeling using the telemac finite-element system. Computers
& Geosciences 53, 105–113.

White, F. M. (1991). Viscous fluid flow mcgraw-hill book company. New York 19(1),
400.

Yang, C. T. & Song, C. (1979). Theory of minimum rate of energy dissipation. Journal
of the Hydraulics Division 105(7), 769–784.

Young, D. F., Munson, B. R., Okiishi, T. H. & Huebsch, W. W. (2010). A brief

introduction to fluid mechanics. John Wiley & Sons.