Machine Translated by Google

Cuerpo de Ingenieros
del Ejército de EE. UU.
Centro de Ingeniería Hidrológica

Aplicación de modelador
Orientación para constante vs.
Inestable y 1D vs 2D vs
Modelado Hidráulico 3D

agosto 2020

Aprobado para publicación pública. Distribución ilimitada. TD­41

Machine Translated by Google

Formulario Aprobado OMB No. 0704­0188

PÁGINA DE DOCUMENTACIÓN DEL INFORME
Se estima que la carga de presentación de informes públicos para esta recopilación de información es de un promedio de 1 hora por respuesta, incluido el tiempo para revisar las instrucciones,
buscar fuentes de datos existentes, recopilar y mantener los datos necesarios y completar y revisar la recopilación de información. Envíe comentarios sobre esta estimación de carga o cualquier otro
aspecto de esta recopilación de información, incluidas sugerencias para reducir esta carga, a la Dirección de Comunicaciones y Servicios Ejecutivos del Departamento de Defensa (0704­0188).
Los encuestados deben tener en cuenta que, independientemente de cualquier otra disposición legal, ninguna persona estará sujeta a ninguna sanción por no cumplir con una recopilación de
información si no muestra un número de control de la OMB actualmente válido.

POR FAVOR NO DEVUELVA SU FORMULARIO A LA ORGANIZACIÓN ANTERIOR.

1. FECHA DEL INFORME (DD­MM­AAAA) 2. TIPO DE INFORME 3. FECHAS CUBIERTAS (De ­ Hasta)
agosto 2020 Documento de formación
4. TÍTULO Y SUBTÍTULO 5a. NÚMERO DE CONTRATO
Guía de aplicación de Modeler para Steady versus
5b. CONCEDER NÚMERO
Inestable y 1D versus 2D versus 3D hidráulico
Modelado 5c. NÚMERO DE ELEMENTO DEL PROGRAMA

6. AUTOR(ES) 5d. NÚMERO DE PROYECTO

Gary Brunner, PE, [Link]., M. ASCE Gaurav
5e. NÚMERO DE TAREA
Savant, Ph.D., PE
Ronald E. Heath 5F. NÚMERO DE UNIDAD DE TRABAJO

7. NOMBRE(S) Y DIRECCIÓN(ES) DE LA ORGANIZACIÓN REALIZADORA 8. INFORME DE ORGANIZACIÓN REALIZADORA NÚMERO TD­41

Centro de Ingeniería Hidrológica del
Instituto de Recursos Hídricos del
Cuerpo de Ingenieros del Ejército de EE. UU. (CEIWR­HEC)
609 Segunda Calle
Davis, CA 95616­4687

Cuerpo de Ingenieros del Ejército de EE. UU.

Centro de Investigación y Desarrollo de Ingenieros
Laboratorio Costero e Hidráulico
3909 Halls Ferry Road
Vicksburg, MS 39180

9. NOMBRE(S) Y DIRECCIÓN(ES) DE LA AGENCIA PATROCINADORA/SEGUIMIENTO 10. ACRÓNIMO(S) DEL PATROCINADOR/MONITOR

11. NÚMERO(S) DE INFORME DEL PATROCINADOR/MONITOR

12. DISTRIBUCIÓN / DECLARACIÓN DE DISPONIBILIDAD

Aprobado para publicación pública. Distribución ilimitada.

13. NOTAS COMPLEMENTARIAS

14. RESUMEN
Todos los modelos, numéricos o físicos a escala, son representaciones simplificadas del mundo real (prototipo). Afortunadamente, existen
numerosos problemas prácticos de ingeniería para los cuales los modelos numéricos simplificados del prototipo son suficientes para
proporcionar descripciones utilizables del comportamiento del sistema. El desafío para el modelador es seleccionar un modelo apropiado para resolver
su problema de ingeniería particular reconociendo al mismo tiempo que el modelo no es una representación perfecta del prototipo. La selección de un
modelo comienza con el desarrollo de una comprensión de qué aspectos del complejo sistema del mundo real son más importantes para el
problema de ingeniería que se aborda. El propósito de este documento es brindar orientación a los ingenieros hidráulicos de nivel medio sobre
cuándo usar el modelado de flujo inestable en lugar del modelado de flujo estacionario; y cómo seleccionar entre modelado unidimensional (1D),
bidimensional (2D) o tridimensional (3D) para un problema determinado.

15. TÉRMINOS DEL TEMA

perfiles de superficies de agua, hidráulica de ríos, flujo estacionario y no estacionario, hidrodinámica unidimensional, bidimensional y
tridimensional, programa de computadora, modelo numérico, HEC­RAS, AdH.

16. CLASIFICACIÓN DE SEGURIDAD DE: 17. LIMITACIÓN DEL 18. NÚMERO DE 19a. NOMBRE DE LA PERSONA RESPONSABLE
a. INFORME b. ABSTRACTO C. ESTA
PÁGINA RESUMEN PÁGINAS
Ud. Ud.
Ud. UU. 114
19b. NÚMERO DE TELÉFONO

Formulario Estándar 298 (Rev. 8/98)

Prescrito por ANSI Std. Z39­18
Machine Translated by Google

Aplicación de modelador
Orientación para constante vs.
Inestable y 1D vs 2D vs
Modelado Hidráulico 3D

agosto 2020

Cuerpo de Ingenieros del Ejército de EE. UU. Cuerpo de Ingenieros del Ejército de EE. UU.
Instituto de Recursos Hídricos Centro de Investigación y Desarrollo de Ingenieros
Centro de Ingeniería Hidrológica Laboratorio Costero e Hidráulico
609 Segunda Calle 3909 Halls Ferry Road
Davis, CA 95616 Vicksburg, MS 39180

(530) 756­1104 (601) 634­2502

(530) 756­8250 FAX [Link]
[Link] [Link]
TD­41
Machine Translated by Google
Machine Translated by Google

TD­41 Tabla de contenido

Tabla de contenido

Lista de tablas ............................................... ................................................. ................................iii

Lista de Figuras ............................................... ................................................. ................................v
Abreviaturas ................................................. ................................................. .......................... ix

Capítulo 1 ................................................. ................................................. ................................ 1­1

Introducción ............................................... ................................................. ...............1­1

Capítulo 2 ................................................. ................................................. ................................ 2­1

Conocimiento del sistema hidráulico y propósito del modelo hidráulico .................2­1

Capítulo 3 ................................................ ................................................. ................................ 3­1

Requerimientos de datos ................................................ ................................................. ..3­1

Capítulo 4 ................................................. ................................................. ................................ 4­1

Salida/Resultados del modelo ................................. ................................................. ...4­1

Capítulo 5 ................................................ ................................................. ................................ 5­1

Modelado de flujo estacionario versus flujo no estacionario ................................. ........................5­1

Capítulo 6 ................................................ ................................................. ................................ 6­1

Modelado unidimensional versus bidimensional ................................. ............6­1

i
Machine Translated by Google

Tabla de contenido TD­41

Capítulo 7 ................................................ ................................................. ................................ 7­1

Modelado bidimensional versus tridimensional ................................. .......... 7­1

Capítulo 8 ................................................ ................................................. ................................ 8­1

Modelos físicos hidráulicos ................................................. .......................................... 8­1
Capítulo 9 ................................................ ................................................. ................................ 9­1
Resumen ................................................. ................................................. ................ 9­1
Referencias ................................................. ................................................. ................................ 1

ii
Machine Translated by Google

TD­41 Lista de tablas

Lista de tablas
Tabla 4­1. Salidas de modelos hidráulicos y nivel de detalle 1D, 2D y 3D. ................................ 4­1
Tabla 7­1. Esquemas sugeridos para el cierre de turbulencias ................................. ................ 7­20
Tabla 9­1. Modelado recomendado para varios sistemas comúnmente modelados. ........................ 9­2

III
Machine Translated by Google

Lista de tablas TD­41

IV
Machine Translated by Google

TD­41 Lista de Figuras

Lista de Figuras
Figura 2­1. Múltiples rutas de flujo para el agua que se mueve dentro de un sistema de diques después de una brecha. .. 2­4
Figura 2­2. Ejemplo de LIDAR detallado y datos de canal (izquierda) versus datos DEM de 10 m (derecha).
................................................. ................................................. ........................................ 2­5

Figura 3­1. Modelo de terreno sin datos de canales submarinos. ................................................ 3 ­2

Figura 3­2. Modelo de terreno con batimetría de canal grabada en el modelo de terreno. ........................ 3­2
Figura 3­3. Ejemplo de uso del suelo y polígonos definidos por el usuario para definir la rugosidad de un modelo 2D.
................................................. ................................................. ........................................ 3­3

Figura 3­4. Terreno detallado y malla de modelado 2D del canal de emisario de la calle 17 en Nueva
Orleans, Luisiana. ................................................. ................................................. ................... 3­4

Figura 4­1. Ejemplo de gráfico de elevación de la superficie del agua 1D versus 2D. ........................................ 4­3
Figura 4­2. Gráfico de velocidad unidimensional (1D) en una sección transversal de ejemplo. ........................ 4­4
Figura 4­3. Gráfico de velocidad bidimensional (2D) en una sección transversal de ejemplo. ........................ 4­4
Figura 4­4. Ejemplo de gráficos de velocidad en 2D y 3D a través de aberturas de puertas. ................................ 4­5
Figura 4­5. Resultados del modelo unidimensional (1D) para un área interior con una ruptura de dique. El color verde a rojo indica el
terreno (elevación baja a alta) y los azules indican la profundidad del agua (el azul oscuro indica mayor
profundidad). ................................................. ................................... 4­5
Figura 4­6. Resultados del modelo bidimensional (2D) para un área interior con una ruptura de dique.
El color verde a rojo indica el terreno (elevación baja a alta) y los azules indican la profundidad del agua (el azul oscuro indica
mayor profundidad). ................................................. ........................ 4­6

Figura 5­1. Corriente empinada (pendiente = 5 pies/milla) con perfiles calculados utilizando flujos máximos y flujos
instantáneos. ................................................. ................................................. ........ 5­3

Figura 5­2. Corriente plana (pendiente = 0,5 pies/milla) con perfiles calculados utilizando flujos máximos y flujos
instantáneos. ................................................. ................................................. ........ 5­3

Figura 5­3. Fuerzas que actúan sobre una masa de agua desde la sección transversal 2 hasta la sección transversal 1 ................. 5­5
Figura 5­4. Familia de ejemplo de curvas de calificación precalculadas para un puente. ................................ 5­8
Figura 5­5. Modelo de flujo estacionario calibrado y modelo de flujo inestable con y sin
Se agregaron pérdidas por contracción y expansión. ................................................. ........................ 5­9
Figura 5­6. Ejemplo de diseño de áreas de flujo ineficaces (polígonos rellenos con líneas diagonales negras) para modelado
1D. ................................................. ................................................. ................ 5­10
Figura 5­7. Hidrograma que entra y sale del tramo de un río con áreas de flujo ineficaz que actúan como 5­11
almacenamiento. ................................................. ................................................. ..........................

Figura 6­1. Definición de símbolos utilizados en las ecuaciones de movimiento 1D y 2D. ........................ 6­3
Figura 6­2. Ejemplo de rotura de un sistema de diques con agua yendo en muchas direcciones. ......... 6­13
Figura 6­3. Ejemplo Rompe presa que sale a una zona extremadamente plana y se extiende.
Las profundidades del agua se muestran en tonos de azul (el azul oscuro indica una mayor profundidad del agua). ..........6­13
Figura 6­4. Bahía inferior del río Columbia con profundidades de agua mostradas en tonos de azul. .......... 6­14
Figura 6­5. Ejemplo de superelevación de la superficie del agua alrededor de una curva pronunciada... 6­15
Figura 6­6. Modelo detallado 2D promediado lateralmente de gradientes de velocidad vertical. El
las direcciones mostradas son x y z. ................................................. ........................................ 6­15

v
Machine Translated by Google

Lista de Figuras TD­41

Lista de Figuras
Figura 6­7. Modelo 2D detallado del flujo que circula alrededor de los muelles desde el andén de una estación de ferrocarril. ... 6­
dieciséis

Figura 6­8. Ejemplo de un sistema fluvial de flujo altamente unidimensional (Allegheny ­

Ríos Monongahela, confluencia en Pittsburgh, PA). El color verde a rojo indica el terreno (elevación baja a alta) y los azules
indican la profundidad del agua (el azul oscuro indica mayor
profundidad). ................................................. ................................................. ................................ 6­17
Figura 6­9. Ejemplo de modelo 1D/2D del río Truckee cerca de Reno, NV. El color verde a rojo indica el terreno (elevación
baja a alta) y los azules indican la profundidad del agua (el azul oscuro indica mayor
profundidad). ................................................. ................................................. 6­18
Figura 6­10. Ejemplo de malla 2D con malla detallada del canal principal, con rejillas alineadas al flujo. El color verde a rojo indica
el terreno (de baja a alta elevación). ........................ 6­21
Figura 6­11. Ejemplo de malla 2D con una malla detallada de mayor resolución alrededor de las estructuras.
El color verde a rojo indica el terreno (de baja a alta elevación).................... ... 6­22
Figura 6­12. Ejemplo de salida de velocidad para un modelo 2D detallado de un puente (la velocidad se superpone al terreno
donde el color verde a rojo indica una elevación baja a alta). ................................ 6­24
Figura 6­13. Ejemplo de un modelo 1D calibrado para el Bajo Río Columbia. ................ 6­27

Figura 7­1. Definición de símbolos utilizados en ecuaciones 3D. ................................................. ....... 7­3
Figura 7­2. Canal profundo rodeado de áreas poco profundas, el rojo indica más profundo y el verde
menos profundo. ................................................. ................................................. ........................ 7­4
Figura 7­3. Estratificación de la salinidad y diferencias de velocidad en un estuario. ................................ 7­4
Figura 7­4. Estratificación de temperatura en un lago. ................................................. ........................ 7­5
Figura 7­5. Campo de velocidad promediado en profundidad calculado en el río Mississippi con vorticidad
transporte. ................................................. ................................................. ........................ 7­6
Figura 7­6. Cambio en la velocidad calculada producida por el método de transporte de vorticidad. ................ 7­7
Figura 7­7. Flujo helicoidal en el recodo de un río. ................................................. ................................ 7­7
Figura 7­8. Diferencia de velocidad entre la superficie y el fondo en el recodo de un río. ........................ 7­8
Figura 7­9. Comparación de velocidad promedio entre 2D (profundidad promediada) y 3D­NHMP
Modelos (sin presión) para un tramo recto de un río con nueve compuertas del mismo tamaño. Colores escalados para
ilustrar patrones, no valores exactos. ................................................. ................ 7­10
Figura 7­10. Comparación de velocidad promedio en profundidad entre modelos 2D y 3D­NHMP (presurizado) para un tramo
recto de un río con nueve compuertas donde la compuerta central penetra 0,5 metros en la columna de agua. Colores
escalados para ilustrar patrones, no valores exactos. ...... 7­11
Figura 7­11. Batimetría y velocidad 2D para un aliviadero inclinado que debe transmitir una probabilidad
Inundación máxima (PMF) de 6.179 m3 /seg. ................................................. ........................ 7­12
Figura 7­12. Los resultados son un aliviadero inclinado (que debe transmitir un PMF de 6179 m3 /seg), velocidad 3D (los colores
representan la velocidad) y superficie del agua (el espesor representa la superficie del agua). ......... 7­13
Figura 7­13. Batimetría, velocidad 2D y perfil de la superficie del agua para el aliviadero inclinado.
(que debe transmitir un PMF de 6.179 m3 /seg). ................................................. ................ 7­13
Figura 7­14. Batimetría para un aliviadero escalonado que debe pasar un PMF de aproximadamente
3.145 m3 /s. ................................................. ................................................. ........................ 7­14
Figura 7­15. Resultados de la simulación hidráulica 3D­NHMP para un aliviadero escalonado (que
debe pasar un PMF de aproximadamente 3.145 m3 /s)................................. ........................ 7­14

vi
Machine Translated by Google

TD­41 Lista de Figuras

Lista de Figuras
Figura 7­16. Resultados de la simulación hidráulica 2D para un aliviadero escalonado (que debe pasar un PMF de
aproximadamente 3.145 m3 /s) para diferentes rugosidades de Manning (n). El espesor indica
profundidad. ................................................. ................................................. ............. 7­15
Figura 7­17. Estrategias de mallado para 3D. ................................................. ................................ 7­17
Figura 7­18. Mallado horizontal para modelos 3D, el color verde a rojo indica la batimetría
características (de más profundas a menos profundas). ................................................. ........................................ 7­18
Figura 7­19. Ejemplo de mallado vertical Arbitrario­Lagrangiano Euleriano (ALE) para 3D
modelos. ................................................. ................................................. ........................ 7­19

viii
Machine Translated by Google

Lista de Figuras TD­41

viii
Machine Translated by Google

TD­41 Abreviaturas

Abreviaturas
1D unidimensional

2D bidimensional

3D tridimensional

3D­NH 3D­No hidrostático

3D­NHMP 3D­Multifásico no hidrostático

adh Hidráulica adaptativa

CERVEZA INGLESA
Euleriano arbitrario­lagrangiano

CEIWR Instituto USACE de Recursos Hídricos

CELRH USACE Distrito de Huntington

cfs pies cúbicos por segundo (ft3 /seg)

cms metros cúbicos por segundo (m3 /seg)

UPC unidad Central de procesamiento

DEM modelo de elevación digital

MDT modelo de terreno digital

[Link], Abreviatura latina estilo APA para "por ejemplo"

CERD Centro de Investigación y Desarrollo de Ingenieros del USACE

ERDC­CHL Laboratorio de Hidráulica y Costas del USACE­ERDC

pie pies

GPU unidades de procesador de gráficos

HEC Centro de Ingeniería Hidrológica CEIWR

HEC­HMS Software del sistema de modelado hidrológico HEC

HEC­RAS Software del sistema de análisis de ríos HEC

HPC computadoras de alto rendimiento

es decir, Abreviatura latina estilo APA para "eso es"

LIDAR Detección de luz y alcance

kg kilogramo

metro metro

NOAA Administración Nacional Oceánica y Atmosférica

Servicio Nacional de Meteorología Servicio Meteorológico Nacional

PMF inundación máxima probable

ix
Machine Translated by Google

Abreviaturas TD­41

Abreviaturas
ppt partes por mil

RAM memoria de acceso aleatorio

segundo segundo

USACE Cuerpo de Ingenieros del Ejército de EE. UU.

WS superficie del agua

WSE elevación de la superficie del agua

X
Machine Translated by Google

TD­41 Capítulo 1

Capítulo 1

Introducción

Todos los modelos, numéricos o físicos a escala, son representaciones simplificadas del mundo real
(prototipo). Cuando se aplican correctamente, la mayoría de los modelos numéricos modernos proporcionarán
soluciones razonablemente precisas de ecuaciones de conservación básicas (conservación de
energía, masa y momento). Sin embargo, la caracterización de sistemas hidráulicos complejos
generalmente requiere parametrizaciones, a menudo puramente empíricas, para describir procesos físicos que
no pueden determinarse directamente a partir de la solución de las ecuaciones de conservación básicas. El
ejemplo clásico de parametrización empírica en el modelado hidráulico es el uso del valor n de Manning
para definir la rugosidad hidráulica.

Afortunadamente, existen numerosos problemas prácticos de ingeniería para los cuales modelos numéricos
simplificados del prototipo son suficientes para proporcionar descripciones utilizables del comportamiento del
sistema. El desafío para el modelador es seleccionar un modelo apropiado para resolver su problema de
ingeniería particular reconociendo al mismo tiempo que el modelo no es una representación perfecta del prototipo.
La selección de un modelo comienza con el desarrollo de una comprensión de qué aspectos del complejo
sistema del mundo real son más importantes para el problema de ingeniería que se aborda.

El propósito de este documento es brindar orientación a los ingenieros hidráulicos de nivel medio sobre
cuándo usar el modelado de flujo inestable en lugar del modelado de flujo estacionario; y cómo seleccionar
entre modelado unidimensional (1D), bidimensional (2D) o tridimensional (3D) para un problema
determinado. Como este documento pretende ser una guía práctica para las aplicaciones de modelos
hidráulicos, no se presentarán derivaciones/discusiones teóricas detalladas de las ecuaciones 1D,
2D y 3D.

Este documento cubrirá lo siguiente:

Conocimiento del Río, u otro sistema hidráulico y Propósito del Hidráulico

Modelado.
Requisitos de datos para modelos estables versus inestables y 1D versus 2D versus
3D. Salida/resultados proporcionados por modelos 1D, 2D y
3D. Modelado de flujo estacionario versus flujo inestable (1D y
2D). Modelado 1D vs 2D
Modelado 2D vs 3D.

El software utilizado en el desarrollo de este documento incluye:

• HEC­RAS: utilizado para ejemplos de modelado 1D y 2D. • AdH: se

utiliza para ejemplos de modelado 2D, 3D y 3D no hidrostático. • OpenFOAM: utilizado
para ejemplos de modelado multifase 3D no hidrostático. • PROTEUS: se utiliza para ejemplos
de modelado multifase 3D no hidrostático.

1­1
Machine Translated by Google

Capítulo 1 TD­41

Agradecimientos
Los capítulos 1 a 6 de este documento fueron escritos por el Sr. Gary W. Brunner (CEIWR­HEC).
Los capítulos 7 y 9 de este documento fueron escritos por el Dr. Gaurav Savant (USACE­ERDC).
El capítulo 8 fue escrito por el Sr. Ronald E. Heath (USACE­ERDC­CHL). Los tres autores
colaboraron en la revisión inicial de los borradores de los capítulos. Los resultados de OpenFOAM
fueron proporcionados por Brian Hall y Nicholas Koutsunis (ambos del CELRH). Los resultados de
PROTEUS fueron proporcionados por el Dr. Christopher E. Kees (USACE­ERDC) y el Sr. Jason H.
Collins (USACE­MVN).

1­2
Machine Translated by Google

TD­41 Capitulo 2

Capitulo 2

Conocimiento del Sistema Hidráulico y

Propósito del modelo hidráulico

Para responder a las preguntas de los enfoques de modelado de flujo estable versus flujo inestable y
1D versus 2D versus 3D, el modelador debe tener conocimiento del sistema hidráulico que se va a modelar, así
como una comprensión clara del propósito del modelo que se va a desarrollar. Cada sistema es único y tendrá
información específica del sitio que se debe considerar para tomar una decisión de modelado adecuada. La
siguiente es una lista de algunas de las cosas que normalmente se deben considerar antes de intentar tomar una
decisión sobre el enfoque de modelado.

Propósito del modelo

El propósito del modelo y el nivel esperado de precisión requerido pueden dictar significativamente el enfoque de
modelado y el nivel requerido de precisión de los datos de origen. Los modelos hidráulicos se desarrollan para todo
tipo de propósitos, tales como: desarrollar un modelo para producir respuestas aproximadas rápidamente; un
estudio de planificación detallado utilizado para evaluar alternativas de estudio; un estudio de diseño en el que el
modelo se utilizará directamente para diseñar una estructura; modelado y mapeo en tiempo real, mapeo
de consecuencias para fallas de presas o diques, etc.

Los modelos que deben desarrollarse rápidamente para proporcionar respuestas aproximadas tenderán a ser
modelos 1D o 2D que no son muy detallados. Estos tipos de modelos se pueden ejecutar en modo de flujo
constante o flujo inestable, para calcular superficies y velocidades del agua que sean aproximadas.
A veces se desarrollan modelos aproximados para operaciones de emergencia cuando no existe ningún
modelo para un área que experimenta un evento importante. Para un caso como este, el desarrollo de un
modelo 2D simple es a veces más rápido que el desarrollo de un modelo 1D, especialmente si ya existe un
modelo digital del terreno y es lo suficientemente detallado como para capturar características hidráulicas
importantes. En otras palabras, diseñar un área de flujo 2D (crear una malla simple y adjuntar algunas
condiciones de contorno (flujo y/o precipitación) a la malla) es muy fácil. Sin embargo, este tipo de modelo 2D no
es un modelo detallado y no debe verse como tal simplemente porque resuelve las ecuaciones de flujo 2D sobre
una malla computacional. Un modelo 2D detallado generalmente requiere tanto esfuerzo como un modelo 1D
detallado, debido al hecho de que el usuario necesita dedicar una cantidad significativa de tiempo a crear una
malla computacional adecuada, definir la rugosidad de la superficie del terreno para toda el área espacial y
calibrar el modelo.

Los estudios de planificación detallados generalmente se realizan con modelos 1D, 2D o combinados 1D/2D que
se calculan en modo de flujo estacionario o no estacionario. Para un nivel de estudio de planificación, las
características específicas del sistema fluvial y la llanura aluvial (como se describe anteriormente en la sección
“Conocimiento del sistema fluvial”), así como los resultados del modelo requeridos para el estudio, dictarán la elección del modelo.
Generalmente, los modelos 3D no se utilizan en estudios de planificación a menos que todo el estudio se dedique
al diseño/análisis de una única estructura hidráulica.

2­1
Machine Translated by Google

Capitulo 2 TD­41

Para estudios detallados a nivel de diseño, es común utilizar un modelo 2D o incluso 1D como herramienta de selección
preliminar (es decir, para reducir el número de alternativas de diseño para modelar en detalle). Sin embargo, los
modelos 3D y los modelos físicos se utilizan con mayor frecuencia para diseñar estructuras hidráulicas, tales como:
aliviaderos, vertederos, aberturas de compuertas, estanques amortiguadores/disipadores de energía; estructuras
de toma de flujo; sistemas de tuberías presurizadas; cruces de arroyos complejos; pasajes para peces; pilares y
estribos para puentes; espigones de ríos/estructuras de entrenamiento; curvas complejas de ríos; y similares. Además,
los modelos 3D también se pueden utilizar para realizar análisis detallados de estructuras existentes.

La previsión de ríos en tiempo real es otra área común que requiere modelos hidráulicos. Debido al hecho de que el
pronóstico en tiempo real requiere modelos que se ejecuten rápidamente, la mayoría de los sistemas de pronóstico
utilizan modelos 1D o combinados 1D/2D. Sin embargo, si el área de pronóstico no es demasiado grande, también
se pueden aplicar modelos 2D más detallados, ya que los requisitos computacionales pueden no ser tan grandes.
La necesidad de respuestas rápidas que sean razonablemente precisas es la principal fuerza impulsora.
Además, los sistemas de pronóstico en tiempo real generalmente solo necesitan modelos hidráulicos para producir
superficies de agua, caudales y mapas de inundaciones, y no distribuciones de velocidad detalladas en 2D
o 3D.

Los modelos de transporte, como los modelos de calidad del agua y sedimentación, pueden ser sensibles a la
precisión de la hidráulica calculada. Los modelos hidráulicos utilizados junto con los modelos de transporte pueden
requerir validación adicional para resolver con precisión los fenómenos de transporte. Por ejemplo, una combinación
específica de coeficientes de rugosidad del canal y de la llanura aluvial en un modelo 1D puede reproducir las
etapas observadas del río. Sin embargo, si la variación modelada de la velocidad del canal con el nivel difiere del
prototipo, los cálculos del transporte de sedimentos pueden producir socavación o deposición excesiva en el canal. En
este caso, es posible que se requieran ajustes adicionales a la distribución del flujo entre el canal y la llanura
aluvial para obtener resultados razonables del modelo de sedimentación.

Conocimiento del sistema fluvial

El modelador debe conocer la descripción física de los canales, áreas de llanuras aluviales, puentes/alcantarillas,
presas, diques, caminos y otras estructuras hidráulicas a las que se aplicará el modelo, antes de elegir un enfoque de
modelado. Algunas preguntas típicas a responder son:

¿Cuál es el tamaño/longitud de los sistemas a modelar? ¿La extensión del sistema que se va a modelar es 1 milla,
10, 50, 100, 500 o 1000 millas? La longitud del sistema a modelar puede dictar qué nivel de modelado se puede utilizar.
Por ejemplo, una aplicación es un modelo de pronóstico para un sistema muy grande. El modelador no debe desarrollar
un modelo 2D detallado de 100 millas de sistema fluvial y esperar que ese modelo se ejecute en un período de tiempo
razonable en una estación de trabajo. Por lo tanto, es posible que sea necesario modelar sistemas
extremadamente grandes en 1D, o en una combinación de 1D y 2D, para tener un modelo que se ejecute en una
cantidad de tiempo computacional razonable. En este caso, puede ser deseable un modelo 2D para informar al
modelo 1D que necesita ejecutarse más rápidamente para soportar un modelo de pronóstico.

¿Cuál es la complejidad del sistema a modelar? Hay muchos factores que forman parte de la cuestión de la complejidad.
¿El sistema es hidráulicamente empinado o tiene partes empinadas del sistema? Los sistemas hidráulicamente empinados
pueden ser más difíciles de modelar que los sistemas fluviales planos. En

2­2
Machine Translated by Google

TD­41 Capitulo 2

En un sistema hidráulicamente empinado hay velocidades más altas, cambios más rápidos en profundidad, área
y velocidad, que ocurren en distancias muy cortas. También puede haber flujo supercrítico y transiciones de flujo
de sub a supercrítico y de supercrítico a subcrítico. Estas condiciones pueden hacer que el cálculo de
soluciones estables a las ecuaciones de conservación sea un desafío y es posible que se requieran técnicas
especiales para resolver las ecuaciones.

Otro factor que forma parte de la cuestión de la complejidad es el número y tipo de estructuras hidráulicas
en el sistema. Las preguntas típicas son:

• ¿Hay muchos puentes que impactan las elevaciones de la superficie del agua durante los flujos más
altos? • ¿Hay puentes muy torcidos o puentes en llanuras aluviales muy amplias con elevaciones
variables de la superficie del agua?
• ¿Existen puentes de flujo sin presión con pilares o estribos complejos? • ¿Hay presas
y vertederos que modelar?
• ¿Existen roturas de presas o diques que deban simularse? • ¿Hay
estructuras de portones que requieren formas únicas en las que se deben operar las puertas?
durante los eventos?

Si hay numerosas estructuras hidráulicas en el sistema que necesitan ser modeladas, entonces se debe
seleccionar un enfoque de modelado que pueda capturar con precisión los aspectos más importantes de cómo
esas estructuras afectan la hidráulica del sistema. El enfoque de modelado también dependerá del nivel de
detalle necesario para modelar una estructura particular. Si se requiere modelar un sistema más grande con
muchas estructuras, puede que no sea posible tener un modelo único que sea muy detallado en el enfoque
para modelar cada estructura individual (es decir, no será posible modelar cada estructura como un modelo 3D
detallado). o incluso estructura 2D). Sin embargo, eso puede no ser tan importante. Si el objetivo del modelo es
predecir las elevaciones de la superficie del agua y los caudales dentro del sistema, entonces el conocimiento
detallado de la distribución de velocidades 2D/3D a través de estructuras individuales no es
necesariamente importante.

¿Se conoce generalmente la trayectoria del flujo del agua para toda la gama de eventos? Comprender la
trayectoria del flujo del agua para toda la gama de eventos es muy importante en el proceso de toma de
decisiones sobre el modelo. Si la trayectoria del flujo del agua está bien definida para toda la gama de eventos,
entonces se puede utilizar un enfoque de modelado 1D, siempre que sea válido para los demás aspectos del modelo.
Sin embargo, si se desconoce la trayectoria del flujo del agua para algunos de los eventos que se van a modelar,
o si el agua puede dividirse e ir en varias direcciones (es decir, el agua que pasa por encima o a través de un
dique puede extenderse en varias direcciones una vez que ingresa al dique), área interior, Figura 2­1),
entonces usar un enfoque de modelado 2D para esas áreas es más apropiado que el modelado 1D. Además, si
la trayectoria del flujo del agua puede cambiar significativamente durante el evento, los enfoques de modelado
2D pueden manejar esto, mientras que los enfoques de modelado 1D no.

2­3
Machine Translated by Google

Capitulo 2 TD­41

Figura 2­1. Múltiples rutas de flujo para el agua que se mueve dentro de un sistema de diques después de una brecha.

¿Existen aspectos únicos del sistema que afectarán significativamente los resultados calculados?
Al estudiar el sistema a analizar, el modelador debe considerar aspectos únicos del sistema que son importantes
para representar con precisión el movimiento del agua y las elevaciones de la superficie del agua/límites de
inundación resultantes. Algunos ejemplos de características únicas del sistema que afectarán
significativamente los resultados del modelo son: el sistema está influenciado por las mareas, de modo
que las mareas oceánicas tienen un impacto significativo en las elevaciones de la superficie del agua; la
velocidad y dirección del viento históricamente han afectado las elevaciones de la superficie del agua; el río
se ve afectado por hielo flotante o atascos de hielo; tiende a haber problemas de escombros durante las
inundaciones y los escombros tienden a acumularse en las estructuras hidráulicas (puentes, alcantarillas, presas,
etc.); hay sistemas de diques que pueden rebasarse o romperse, donde es necesario abordar el recorrido del
flujo interior; y existen estructuras hidráulicas únicas que requieren modelado especializado u operaciones de compuerta.

Fuentes y exactitud de los datos

La fuente y el nivel de precisión de los datos que se utilizan para desarrollar el modelo son muy importantes para la
decisión del enfoque de modelado y la precisión de los resultados del modelo. Específicamente, el nivel de detalle
y precisión de los datos del terreno, los datos batimétricos, los datos de la sección transversal, la información
de los diques y los datos de la estructura hidráulica son importantes para decidir qué tan detallado se puede
desarrollar un modelo. Por ejemplo, si no existen datos detallados del terreno y datos batimétricos (es decir, solo
hay disponibles DEM de 10 metros, pero hay secciones transversales topográficas, consulte el ejemplo en la
Figura 2­2), entonces el aumento percibido en la precisión al usar un modelo 2D sobre un modelo 1D puede que en realidad no

2­4
Machine Translated by Google

TD­41 Capitulo 2

ser realizado. El tipo de datos y su nivel de precisión influirán en la calidad de la elección del modelo.

Además, en el proceso de selección del modelo también se debe considerar el nivel de precisión de las
condiciones hidrológicas/de límites utilizadas para impulsar el modelo. Si sólo se dispone de estimaciones de los
caudales máximos y no se conocen los hidrogramas completos en las ubicaciones externa e interna del modelo, es
posible que no sea posible modelar el flujo inestable.

Figura 2­2. Ejemplo de LIDAR detallado y datos de canal (izquierda) versus datos DEM de 10 m (derecha).

Duración de los Eventos a Modelar

Los tipos de eventos varían desde: un caudal constante de una magnitud específica; un tipo de evento de
escorrentía de lluvia normal; Inundaciones repentinas; ruptura de presas y/o diques; liberaciones de flujo
desde estructuras hidráulicas; etc. La duración de un evento depende del tamaño de la cuenca/sistema fluvial
que se modela, así como del propósito del estudio. Es posible que algunos propósitos de estudio solo requieran la
modelización de flujos máximos para una variedad de eventos. De hecho, todos los modelos 1D, 2D y 3D se pueden
utilizar en un modo de flujo constante. Sin embargo, los modelos 1D se utilizan generalmente para modelar grandes
extensiones de sistemas fluviales basándose en caudales máximos derivados de modelos hidrológicos o datos
observados. En general, los modelos 2D y 3D se utilizan en un “modo de flujo estacionario” (ejecutando un flujo
inestable usando un flujo constante y/o condiciones de contorno del hidrograma de etapa) para tramos cortos de río,
o para el diseño y análisis de estructuras hidráulicas.

Para el modelado de flujo inestable, la duración de los eventos puede tener un impacto en el tipo de enfoque
de modelado hidráulico. Si los eventos que se modelan tienen una duración más corta (es decir, de 1 a 3 días, o
menos de una semana), entonces los modelos 1D, 2D o incluso 3D pueden seguir siendo una opción viable, siempre
que el área que se modela sea pequeña. A medida que el dominio del modelo se hace más grande, es posible que los
modelos 3D ya no sean una opción viable para ubicaciones con eventos de más de una semana. Como la
duración del evento varía de unas pocas semanas a meses, si el sistema fluvial tiene una longitud significativa,
entonces incluso un enfoque de modelado 2D puede no ser viable debido a la duración del tiempo computacional
requerido para ejecutar un evento de este tipo en un sistema fluvial. de tamaño significativo.

2­5
Machine Translated by Google

Capitulo 2 TD­41

En este tipo de situación es cuando el modelado 1D o el modelado combinado 1D/2D pueden ser una mejor
opción. Luego, para los análisis de períodos de registro, en los que se requieren uno o más años de simulación,
generalmente se utilizan modelos 1D, pero posiblemente modelos combinados 1D/2D, si las áreas de flujo 2D son
pequeñas o solo entran en juego durante eventos de flujo grandes. dentro del período de registro (es decir, las áreas
de flujo 2D se utilizan para modelar las áreas detrás de los sistemas de diques).

Por lo general, los modelos 1D, con grandes extensiones espaciales y temporales, generalmente se ejecutan en
procesadores únicos en máquinas de escritorio. Por otro lado, los modelos 2D y 3D, con extensiones
espaciales y temporales superiores a decenas de millas y unos pocos días, requieren la utilización de
máquinas multiprocesador, como computadoras personales a nivel de estación de trabajo (con muchos núcleos) y
computación de alto rendimiento ( sistemas HPC). Los recientes avances en modelos 2D y 3D han permitido
la simulación de cientos de kilómetros y años de simulaciones, pero esto requiere la utilización de importantes
recursos computacionales (es decir, supercomputadoras).

Resultados del modelo requeridos

Casi todos los estudios que requieren un modelo hidráulico necesitan elevaciones, profundidades y caudales de la
superficie del agua calculados. Todos los enfoques de modelado pueden producir este tipo de resultados, pero con
distintos niveles de precisión. Específicamente, los modelos 1D calculan elevaciones promedio de la superficie del
agua en cada sección transversal y área de almacenamiento dentro del modelo. Por el contrario, los modelos 2D y
3D tienen superficies de agua que varían espacialmente según el tamaño y la cantidad de celdas/elementos/
nodos utilizados en la malla computacional. Las profundidades se pueden calcular a partir de cualquiera de los mapas
de inundación y elevaciones de la superficie del agua resultantes del modelo, pero dependen de la precisión del
modelo de terreno subyacente. Los caudales de los modelos 1D generalmente se informan como flujo total en
cada sección transversal o como caudal en el canal principal, sobrebanco izquierdo y sobrebanco derecho (aunque el
caudal se puede dividir aún más según el transporte a través de la sección transversal). y la suposición de
que el flujo es perpendicular a la sección transversal). Los caudales de los modelos 2D se pueden adquirir a lo largo
de cualquier línea definida por el usuario dentro de la malla computacional.

Muchos estudios también requieren información sobre la velocidad por diversas razones. Específicamente, los
modelos 1D solo producen velocidades promediadas horizontal y verticalmente. Estas velocidades a
menudo se informan por separado para el canal principal y las áreas de sobrebanco izquierda y derecha. Al igual
que con el flujo, las velocidades se pueden discretizar aún más según el transporte de la sección
transversal y la suposición de que el flujo es perpendicular a la sección transversal. Sin embargo, en zonas de
contracciones y expansión detalladas, por ejemplo el flujo a través de la abertura de un puente, las velocidades
producidas por los modelos 1D no son tan precisas como las de los modelos 2D y 3D. Las distribuciones de velocidad
detalladas para canales/llanuras de inundación normales, así como las velocidades detalladas a través de
contracciones/expansiones, alrededor de curvas cerradas y alrededor de estructuras hidráulicas, requieren
enfoques de modelado 2D y posiblemente incluso 3D. Sin embargo, el modelador debe desarrollar una
malla computacional con suficientes celdas/caras/elementos para producir una distribución de velocidad detallada
para las ubicaciones y tipos de estructura deseados.

También es posible que se requiera otro tipo de información sobre los modelos hidráulicos. Generalmente, la
mayoría de los modelos hidráulicos 1D generan una amplia gama de variables hidráulicas en cada sección transversal
y estructura hidráulica (HEC­RAS genera cerca de 300 variables hidráulicas diferentes en cada sección transversal
para cada paso de caudal/tiempo). Por otro lado, los modelos 2D y 3D generalmente no producen esto.

2­6
Machine Translated by Google

TD­41 Capitulo 2

tipo de salida directamente, pero puede tener formas de llegar a la salida a partir del posprocesamiento de los resultados del
modelo básico de profundidades, elevaciones de la superficie del agua y velocidades.

Otros tipos de resultados, como los tiempos de llegada; duraciones de flujo/profundidad; porcentaje de tiempo inundado;
tiempos de residencia; etc., todos requieren un modelado de flujo inestable, que puede ser en forma de modelo de flujo
inestable 1D, 2D o 3D.

Experiencia del modelador

La cantidad y el tipo de experiencia que tenga el modelador también afectarán la elección del modelo que se utilizará. En
general, los modelos de flujo estacionario 1D son los más fáciles de usar y comprender los resultados. Pasar al modelado de
flujo inestable requiere más conocimiento de hidráulica y también más conocimiento de algoritmos de soluciones numéricas,
como técnicas de soluciones de diferencias finitas, volúmenes finitos y elementos finitos. El modelado de flujo
inestable requiere más conocimiento de la propagación de las olas y de cómo cambiará de forma un hidrograma a
medida que se mueve de un punto del sistema a otro. Las técnicas de solución numérica para el modelado de flujo
inestable 1D tienden a ser metodologías de diferencias finitas o de volúmenes finitos. Además, los modelos 2D y 3D son en
su mayoría enfoques de elementos finitos o volúmenes finitos.

Comprender estas técnicas de solución numérica es importante cuando se utilizan dichos modelos.
Elegir un paso de tiempo computacional, un espaciamiento de sección transversal (modelado 1D) o un tamaño de celda
(modelado 2D y 3D) apropiado es importante para lograr soluciones numéricas que no atenúen artificialmente (a
menudo llamado difusión numérica) el hidrograma a medida que se mueve a través del sistema. . Posteriormente, el modelado
de flujo inestable 2D y 3D también requiere que el usuario comprenda el modelado de turbulencias y sus efectos en el
campo de flujo. Además, las fuerzas externas sobre el sistema, como la rotación de la Tierra (efecto Coriolis) y las tensiones
del viento, pueden ser importantes y sólo pueden tenerse en cuenta en enfoques de modelado 2D y 3D.

2­7
Machine Translated by Google

Capitulo 2 TD­41

2­8
Machine Translated by Google

TD­41 Capítulo 3

Capítulo 3

Requerimientos de datos

Los requisitos de datos pueden variar significativamente para los enfoques de modelado unidimensional (1D),
bidimensional (2D) y tridimensional (3D), así como para el modelado de flujo estacionario versus no estacionario.
La cantidad y calidad de los datos disponibles pueden dictar el tipo de modelado que se puede
lograr. Las principales áreas en las que los requisitos de datos pueden ser diferentes son: información
topográfica (datos del terreno); vegetación de canales y llanuras aluviales o uso de la tierra (definiendo
valores de rugosidad); infraestructura de drenaje subterráneo; información de la estructura de la
superficie; hidrología requerida; condiciones de borde; y datos de calibración.

Datos del terreno

Los requisitos del terreno pueden variar desde la definición de un modelo solo con secciones transversales
(modelado 1D) hasta modelos de terreno detallados de todo el canal y la llanura aluvial, así como
características tales como: caminos, diques, muros de contención, estructuras de formación de canales,
etc. Los modelos se pueden impulsar teniendo solo secciones transversales en las ubicaciones necesarias para
calcular una superficie de agua precisa y la ruta del hidrograma. Sin embargo, el modelado 2D o 3D requiere
un modelo de terreno (Digital Elevation Model, DEM o Digital Terrain Model, DTM) de todo el sistema.

Además, la precisión de ese modelo de terreno afectará directamente la precisión del enfoque de modelado
2D o 3D. Por ejemplo, si se desarrolla un modelo de terreno detallado a partir de datos LIDAR, pero los
datos del canal submarino (batimetría) no están definidos con precisión, es mucho más fácil modificar las
secciones transversales en un marco de modelado 1D que modificar todo el modelo de terreno para incorporar
el batimetría del canal.

En algunos casos, una tarea necesaria puede ser fusionar datos batimétricos del sonar multihaz con datos
LIDAR terrestres y aéreos para producir un modelo de terreno perfecto. Específicamente, los modelos 2D y 3D
solo serán precisos si el terreno incluye una representación precisa del terreno submarino del canal
principal y cualquier estructura que influya en el campo de flujo. Consulte la Figura 3­1 y la Figura 3­2 para ver
un ejemplo de datos de terreno con y sin datos de canal grabados en el modelo de terreno.

3­1
Machine Translated by Google

Capítulo 3 TD­41

Figura 3­1. Modelo de terreno sin datos de canales submarinos.

Figura 3­2. Modelo de terreno con batimetría de canal grabada en el modelo de terreno.

Coeficientes de rugosidad
Los requisitos de datos para definir la rugosidad también pueden variar. En general, se requiere
conocimiento de la vegetación y el uso de la tierra para todo el dominio de modelado para todos los enfoques de modelado.

3­2
Machine Translated by Google

TD­41 Capítulo 3

Sin embargo, la rugosidad para los enfoques de modelado 1D solo debe definirse en cada sección transversal.
Para un enfoque de modelado 1D, los coeficientes de rugosidad se pueden definir sección transversal por sección
transversal, o el modelador puede diseñar información espacial sobre vegetación/uso de la tierra. Las cuadrículas y
polígonos de uso del suelo se relacionan con los valores de rugosidad y luego los valores de rugosidad se extraen en
la intersección de las secciones transversales y las capas de rugosidad espacial.

En general, los enfoques de modelado 2D y 3D requieren que el modelador diseñe información

espacial de vegetación/uso de la tierra y la relacione con los valores de rugosidad. Luego, la rugosidad se define
espacialmente para cada celda de cálculo/cara de elemento de la malla computacional 2D/3D a lo largo del límite
de la superficie del terreno. Además, la rugosidad del canal principal debe definirse con polígonos separados
definidos por el usuario. Esto generalmente es necesario porque la mayoría de los conjuntos de datos de uso/
cobertura del suelo no definen el canal en detalle y solo lo definen con un único tipo de uso del suelo. En la Figura
3­3 se muestra un ejemplo de un modelo 2D con rugosidad definida con el uso del suelo en las áreas sobre el banco
y polígonos definidos por el usuario para el canal principal.

Figura 3­3. Ejemplo de uso del suelo y polígonos definidos por el usuario para definir la rugosidad de un modelo 2D.

Sistemas de drenaje subterráneo

El modelado detallado de los sistemas de drenaje subterráneos puede o no ser necesario para un estudio en
particular. Para modelos detallados de áreas urbanas, a menudo se requiere modelar los sistemas de drenaje
subterráneo. Debido a las complejas interconexiones de los sistemas de tuberías subterráneas, el modelado
de flujo constante generalmente no es una opción. Muy a menudo, el modelado de flujo inestable 1D se utiliza
para modelar sistemas de drenaje subterráneo. Esta generalización se debe a que los flujos y velocidades en este
tipo de sistemas son de naturaleza muy unidimensional.
Por lo tanto, rara vez se realiza el modelado 2D o 3D de los sistemas de drenaje subterráneo existentes.
Independientemente del enfoque de modelado, los datos necesarios para modelar sistemas de drenaje subterráneos
son los mismos.

3­3
Machine Translated by Google

Capítulo 3 TD­41

Estructuras Hidráulicas
Los requisitos de datos para modelar estructuras hidráulicas también variarán entre los enfoques de modelado 1D, 2D
y 3D. Para los modelos 1D, las estructuras pueden definirse con ecuaciones semiempíricas o curvas de
clasificación e insertarse como condiciones de contorno internas entre secciones transversales.
Los datos necesarios para modelar estructuras hidráulicas se basan en el modelo computacional hidráulico que se
utiliza para modelar la estructura (por ejemplo, una ecuación de vertedero solo necesita un perfil de línea central de la
parte superior del vertedero, la forma del vertedero y un coeficiente del vertedero. Mientras que un la puerta sólo necesita
un ancho, una altura, una elevación invertida y un coeficiente de puerta).

Además, los enfoques de modelado 2D también pueden utilizar un modelado de estructuras hidráulicas similar dentro del
dominio 2D; sin embargo, este es un enfoque 1D para modelar la estructura dentro del área 2D. El verdadero modelado
de flujo 2D/3D de una estructura hidráulica requerirá un modelado detallado del terreno/superficie de la
estructura hidráulica desde la entrada aguas arriba, a través de la estructura, hasta la salida aguas abajo. Además,
es necesario definir con mayor precisión la rugosidad de toda la superficie de la estructura. En la Figura 3­4 se muestra
un ejemplo de un modelo 2D muy detallado de la estructura del emisario de la Calle 17 en Nueva Orleans, Luisiana. Este
modelo se utilizó como herramienta de selección para reducir el número de diseños posibles a modelar con más
detalle con un modelo 3D.

Figura 3­4. Terreno detallado y malla de modelado 2D del canal de desagüe de la calle 17 en Nueva Orleans, Luisiana.

Los flujos de campo cercano en estructuras hidráulicas pueden ser no hidrostáticos. En tales casos, la aplicación de
modelos hidrostáticos tenderá a sobreestimar las pérdidas de energía en las proximidades de la estructura y puede no
reproducir con precisión los patrones de flujo prototipo inmediatamente aguas abajo de la estructura.
estructura.

3­4
Machine Translated by Google

TD­41 Capítulo 3

Hidrología

La hidrología utilizada para impulsar los modelos de flujo estable frente a los de flujo inestable es diferente. El flujo constante
requiere que el usuario defina los caudales para todo el sistema utilizando un modelo hidrológico o datos medidos. El modelado
hidrológico utiliza técnicas de enrutamiento más simples que el enrutamiento hidráulico (como: Muskingum, Modified Puls y
Muskingum­Cunge). Los modelos de flujo inestable utilizan hidrogramas en todos los límites externos aguas arriba, así como
para cualquier ubicación lateral e interna requerida. Y los modelos de flujo inestable utilizan métodos de enrutamiento
hidráulico más detallados para resolver cómo se mueve el agua a través del sistema. La hidrología requerida para los modelos
1D, 2D y 3D es prácticamente la misma, ya que se basa mucho más en el tamaño/extensión del dominio de modelado y no
en el enfoque de modelado. Sin embargo, la forma de ingresar los datos de flujo en el sistema puede ser diferente,
dependiendo del enfoque de modelado (es decir, condiciones de contorno internas 1D versus 2D).

Condiciones de borde

Las condiciones de contorno internas y aguas abajo también pueden requerir diferentes cantidades de datos para los enfoques
de modelado de flujo estacionario versus inestable. Los modelos de flujo estacionario generalmente solo necesitan elevaciones
de la superficie del agua aguas abajo para calcular cada perfil. Los modelos de flujo inestable pueden requerir hidrogramas de
etapa completa o curvas de clasificación de caudal de etapa. La cantidad de datos depende de la ubicación y del tipo de
condición de contorno que se aplica. Si un modelo está influenciado por las mareas, entonces se requieren hidrogramas de
etapa completa para enfoques de modelado de flujo inestable.
Sin embargo, en ubicaciones de ríos que están más controladas por el caudal del río y las fuerzas de gravedad/fricción, las curvas
de clasificación o los enfoques de condiciones de contorno de profundidad normal (ecuación de Manning) se pueden aplicar para
flujo estable y flujo inestable de la misma manera.

Datos de calibración

Es necesario calibrar el modelo independientemente del tipo de modelo. La cantidad y el tipo de datos observados
necesarios para la calibración del modelo también pueden variar entre los enfoques de modelado.
Los modelos de flujo estacionario utilizan elevaciones máximas de la superficie del agua y, opcionalmente, magnitudes de
velocidad para compararlas con los valores calculados. Los modelos de flujo inestable necesitan hidrogramas de etapa
y de flujo completo en los medidores, así como marcas de agua máxima cuando estén disponibles.

Generalmente, los enfoques de modelado 1D solo utilizan superficies de agua observadas e hidrogramas de flujo en los
medidores, y luego marcas de agua altas entre las ubicaciones medidas. Sin embargo, los modelos 2D y 3D también pueden
necesitar información observada sobre la velocidad o la distribución del flujo. Además, para el modelado 2D, la velocidad
observada debe medirse espacialmente a través del río y la llanura aluvial para calibrar la distribución de velocidad calculada. A
menudo también se requieren extensiones de inundación en flujos particulares para calibrar modelos 2D más detallados. Los
usuarios deberán consultar la documentación específica del modelo que están utilizando para obtener más información sobre la
calibración de ese modelo.

3­5
Machine Translated by Google

Capítulo 3 TD­41

3­6
Machine Translated by Google

TD­41 Capítulo 4

Capítulo 4

Salida/resultados del modelo

Los requisitos para los resultados del modelo hidráulico, así como el nivel de detalle, influirán en el tipo de
modelo utilizado para un estudio. Por ejemplo, si se necesitan velocidades detalladas en la base de un
dique, o alrededor de un pilar o estribo de un puente, entonces el modelado 1D no puede proporcionar ese tipo
de detalle y será necesario utilizar el modelado 2D o 3D. Entonces, las preguntas que los modeladores deben
hacerse al comienzo de un estudio son: cuáles son todos los resultados hidráulicos requeridos para este
estudio, qué nivel de detalle se necesita para cada uno de los resultados hidráulicos y qué nivel de
precisión se espera o se desea. para cada resultado hidráulico?

Diferentes niveles de detalle en las salidas hidráulicas

La siguiente es una tabla de salidas/resultados hidráulicos comunes que a menudo se solicitan en

estudios de modelado hidráulico (Tabla 4­1). Además, la Tabla 4­1 describe el nivel de detalle de los modelos
1D, 2D y 3D para esa salida hidráulica.

Tabla 4­1. Salidas de modelos hidráulicos y nivel de detalle 1D, 2D y 3D.

Hidráulico Flujo inestable 1D Modelado de flujo inestable Modelado de flujo inestable

Salida/Resultados Modelado 2D WSE con 3D WSE con

variación horizontal. Un variación horizontal. Un

Superficie máxima del agua WSE promedio único por sección
WSE para cada celda. WSE para cada celda/nodo.
Elevación (WSE) transversal y área de almacenamiento.

WSE promedio versus tiempo para Variación horizontal de Variación horizontal de

Hidrogramas de etapa secciones transversales y almacenamiento WSE frente al tiempo para WSE frente al tiempo para
áreas. cada celda/nodo computacional. cada celda/nodo computacional.

Flujos máximos en Flujos máximos en

Flujo máximo en cada sección
líneas de salida/perfil líneas de salida/perfil
Tasas de flujo máximo transversal y estructuras
definidas por el usuario definidas por el usuario
hidráulicas.
y estructuras hidráulicas. y estructuras hidráulicas.

Flujo versus tiempo en Flujo versus tiempo en

Flujo vs tiempo en cada sección líneas de salida/perfil líneas de salida/perfil
Hidrogramas de flujo transversal, condiciones de contorno definidas por el definidas por el
y estructuras hidráulicas. usuario, condiciones de contorno usuario, condiciones de contorno
y estructuras hidráulicas. y estructuras hidráulicas.

Velocidades promedio para el canal

Velocidades que varían Velocidades variables
principal, sobrebanco izquierdo y
horizontalmente pero horizontal y
sobrebanco derecho. Una
Velocidades promedian verticalmente. Una verticalmente. Una velocidad
mayor discretización se basa en
velocidad promedio por cara de malla computacional.
subdivisiones basadas
para cada cara de celda/elemento.
en transporte.

4­1
Machine Translated by Google

Capítulo 4 TD­41

Hidráulico Flujo inestable 1D Modelado de flujo inestable 2D Flujo inestable 3D

Salida/Resultados Modelado La dirección del Modelado
flujo horizontal se calcula en función

de los detalles del terreno y la malla

El modelador debe definir la Las direcciones
computacional.
direcciones de flujo y dirección del flujo al trazar los tramos de tridimensionales y los

patrones los ríos y las áreas de almacenamiento. patrones de flujo se calculan

Se pueden determinar patrones
directamente.
de circulación horizontal (remolinos).

Los tiempos de llegada de las

inundaciones se
basan en velocidades
Los tiempos de llegada de las inundaciones
Los tiempos de llegada de las tridimensionales y patrones
se basan en los cálculos de
inundaciones se de flujo, así como en las
velocidades promedio 1D y la
basan en velocidades elevaciones de la superficie del
interpolación de superficies de agua entre
Horas de llegada de las inundaciones bidimensionales y patrones agua dentro de cada
secciones transversales. El enrutamiento
de flujo, así como en celda vertical. Actualmente, el
del conjunto de niveles no se puede
elevaciones de la superficie del modelado 3D no se utiliza
utilizar para estimar los tiempos de llegada
agua dentro de cada celda/nodo. con tanta frecuencia para los
de inundaciones en áreas de almacenamiento.
tiempos de llegada en
situaciones fluviales, debido
a los elevados requisitos/tiempos computacionales.

La profundidad se calcula a partir La profundidad se calcula a partir

La profundidad se calcula a partir
de las elevaciones de la superficie del de la superficie del agua de la
de la superficie del agua de la
agua interpoladas espacialmente menos celda menos las elevaciones del
celda/nodo menos las elevaciones
Mapeo de peligros la elevación del terreno en esa terreno en cada ubicación.
del terreno en cada ubicación.
Profundidad x Velocidad ubicación. La velocidad se La velocidad se interpola a
La velocidad se promedia
interpola a partir de la interpolación de partir de velocidades
verticalmente en cada
velocidades promediadas 1D calculadas espacialmente
ubicación.
descritas anteriormente. en 2D en cada cara de celda/nodo.

El límite de la superficie del agua se

El límite de profundidad cero se El límite de profundidad cero se
calcula en cada sección
calcula para cada celda/nodo calcula para cada celda/nodo
transversal, luego se crea una superficie
Inundación que esté parcialmente húmedo. que esté parcialmente húmedo.
de interpolación y se cruza
Límites Estos límites se fusionan para Estos límites se fusionan para
con el terreno para encontrar el límite del
formar polígonos continuos. formar polígonos continuos.
agua (elevación de profundidad cero).

Para las secciones transversales 1D,

la sección transversal se divide en Las propiedades hidráulicas se
Para celdas/nodos 2D, es el
sectores definidos por el usuario y promedian verticalmente, luego se
Esfuerzo cortante esfuerzo cortante promedio en
luego se calculan los valores promedio calcula el esfuerzo cortante
calculado como: cada cara, luego se interpola
para cada sector. Los valores se promedio en cada cara y luego se
(γ R Sf). entre las caras.
interpolan entre secciones transversales interpola entre las caras.
utilizando la superficie de
interpolación de secciones transversales.

Para las secciones transversales 1D,

Las propiedades hidráulicas se
la sección transversal se divide en Para celdas/nodos 2D, es la
La potencia de la promedian verticalmente, luego la
sectores definidos por el usuario y velocidad promedio
corriente se calcula potencia promedio de la
luego se calculan los valores promedio multiplicada por el esfuerzo
como: velocidad promedio corriente se calcula en
para cada sector. Los valores se cortante promedio en cada cara,
multiplicada por el esfuerzo cada cara y luego se interpola
interpolan entre secciones transversales luego se interpola entre
cortante promedio entre las caras.
utilizando la superficie de las caras.

interpolación de secciones transversales.

4­2
Machine Translated by Google

TD­41 Capítulo 4

En la Figura 4­1 se muestra un gráfico de la elevación de la superficie del agua (WSE), en la misma ubicación,
a partir de un modelo 1D y 2D (mismo caudal). Como se muestra en la Figura 4­1, el modelo 1D tiene una
línea horizontal (azul) que atraviesa toda la sección transversal de la superficie del agua. Sin embargo, la
superficie del agua varía respecto al modelo 2D (línea verde). Esta ubicación de ejemplo está en el extremo aguas
arriba de una curva hacia la izquierda, razón por la cual los resultados del modelo 2D muestran una superficie de
agua más alta en el lado derecho del perfil del terreno.

Figura 4­1. Ejemplo de gráfico de elevación de la superficie del agua 1D frente a 2D.

Los gráficos de velocidad de los resultados del modelo 1D (gráfico superior) y 2D (gráfico inferior) se muestran
en la Figura 4­2 y en la Figura 4­3 para el mismo ejemplo de ubicación y caudal. Como se indicó anteriormente,
esta ubicación de ejemplo está directamente aguas arriba de una curva (a la izquierda) en el río. La Figura 4­2 y
la Figura 4­3 ilustran que la superficie del agua y las velocidades para el resultado del modelo 2D contienen más
detalles en esta ubicación, donde el resultado 1D muestra una distribución de velocidad mucho más
uniforme, debido al enfoque aplicado para distribuir un resultado de velocidad 1D en el espacio.

4­3
Machine Translated by Google

Capítulo 4 TD­41

Figura 4­2. Gráfico de velocidad unidimensional (1D) en una sección transversal de ejemplo.

Figura 4­3. Gráfico de velocidad bidimensional (2D) en una sección transversal de ejemplo.

En la Figura 4­4 y la Figura 4­5 se muestran gráficos de velocidad de un modelo 2D y 3D para el mismo conjunto
de compuertas, con el mismo caudal. Como puede ver en los gráficos a continuación, los gráficos de velocidad 3D
son una descripción más precisa del movimiento real del fluido a través de las compuertas. Sin embargo, ese nivel
de detalle puede ser necesario o no para cualquier estudio en particular.

4­4
Machine Translated by Google

TD­41 Capítulo 4

Figura 4­4. Ejemplo de gráficos de velocidad 2D (arriba) y 3D (abajo) a través de las aberturas de las puertas.

Los resultados proporcionados en la Figura 4­5 y la Figura 4­6 son ejemplos de mapas de áreas inundadas
para un área interior protegida por un dique. Específicamente, la Figura 4­5 muestra los resultados de un modelo 1D
en el que el área interior se modeló con áreas de almacenamiento interconectadas. Por otro lado, la Figura 4­6
proporciona los resultados de un modelo 2D de la misma área y el mismo flujo que llega a través de una
brecha aguas arriba del dique.

Figura 4­5. Resultados del modelo unidimensional (1D) para un área interior con una ruptura de dique. El color verde a rojo indica el
terreno (elevación baja a alta) y los azules indican la profundidad del agua (el azul oscuro indica mayor profundidad).

Como se muestra en la Figura 4­5, los resultados 1D muestran agua desconectada. Esto se debe a que las áreas de
almacenamiento se llenan automáticamente desde la elevación más baja hasta la más alta, con una superficie
de agua horizontal. Para la misma ubicación, los resultados del modelo 2D (Figura 4­6) muestran trayectorias de flujo
terrestre que conectan todas las áreas interiores. Entonces, para este área de ejemplo, el modelo 2D proporciona

4­5
Machine Translated by Google

Capítulo 4 TD­41

resultados que son más realistas sobre cómo ocurriría la inundación, así como para calcular las
profundidades, velocidades y tiempos de llegada de las aguas de la inundación.

Figura 4­6. Resultados del modelo bidimensional (2D) para un área interior con una ruptura de dique. El color verde
a rojo indica el terreno (elevación baja a alta) y los azules indican la profundidad del agua (el azul
oscuro indica mayor profundidad).

Nivel de precisión esperado/deseado

Los diferentes tipos de estudios tienen diferentes niveles de precisión esperada/deseada. Por ejemplo, en
situaciones de emergencia, puede ser necesario desarrollar muy rápidamente un modelo aproximado con datos
limitados. El nivel esperado de precisión de dicho modelo no es alto y, por lo tanto, el enfoque de modelado
puede ser menos detallado (es decir, modelado 1D o 2D con menos detalle). Proporcionar algún nivel de resultados
hidráulicos (mapas de inundaciones, tiempos de llegada y velocidades) es mejor que no tener ninguna
información. Por el contrario, si el modelador está realizando un estudio de diseño muy detallado, entonces el nivel
de precisión esperado/deseado es muy alto. Para este tipo de estudio, se debería dedicar más tiempo a adquirir
datos detallados; realizar modelado detallado; completar análisis de calibración/verificación; y completar análisis de
riesgo e incertidumbre.

El nivel de precisión esperado/deseado también afectará el enfoque de modelado. Por ejemplo, si se espera que el
estudio tenga velocidades bidimensionales y tridimensionales precisas en una ubicación, entonces será necesario
utilizar enfoques de modelado 2D y 3D, ya que el modelado 1D no puede proporcionar altos niveles de precisión
para ese tipo de salida hidráulica. Por lo tanto, el conocimiento de todas las salidas hidráulicas requeridas, así
como el nivel de precisión esperado/deseado, es muy importante para tomar una decisión sobre el enfoque de
modelado.

4­6
Machine Translated by Google

TD­41 Capítulo 5

Capítulo 5

Modelado de flujo estacionario versus flujo inestable

Este capítulo analiza las diferencias entre el modelado de flujo estacionario y no estacionario.
Específicamente, este capítulo describe: la definición de flujo estacionario e inestable; supuestos utilizados en el
modelado de flujo estacionario; enrutamiento hidrológico versus hidráulico; diferencias en los cálculos hidráulicos; diferencias
en las estrategias de calibración; y limitaciones del modelado de flujo constante.

Definiciones

El modelado de flujo constante se basa en el uso espacial de un conjunto específico de caudales y luego en el cálculo de
elevaciones de la superficie del agua, velocidades, etc., en función de esos caudales. Por una definición estricta:

Flujo constante: el flujo (es decir, profundidad, velocidad, descarga) no cambia con el tiempo.

Debido a que los flujos de los ríos suelen ser turbulentos y la velocidad en cualquier punto fluctúa constantemente,
la definición de flujo estacionario puede ampliarse para incluir flujos en los que la velocidad media y la profundidad en
cualquier punto pueden considerarse constantes durante el período de tiempo que se modela.
Asimismo, una suposición de flujo constante puede ser razonable para los cálculos hidráulicos si el flujo cambia
gradualmente con el tiempo. Algunos ejemplos de flujo constante son:

Flujo en un canal con descarga constante desde aguas arriba.

Caudal natural de un río en el que el caudal cambia gradualmente con respecto al tiempo. Modelar un tramo de

río muy corto, de modo que el caudal sea efectivamente constante

en todo el tramo en cualquier momento.

Modelar una red fluvial en la que el flujo en cada segmento de la red sea efectivamente
constante.

El modelado de flujo inestable se basa en proporcionar hidrogramas completos en todos los puntos aguas arriba del
sistema fluvial, así como para cualquier punto de afluencia lateral, luego las ecuaciones de flujo inestable se utilizan para
encaminar los hidrogramas y al mismo tiempo se calculan las elevaciones de la superficie del agua. La definición estricta
de flujo inestable es:

Flujo inestable: el flujo (es decir, profundidad, velocidad, descarga) cambia con el tiempo.

Algunos ejemplos de flujo inestable son:

• Represas y diques rompen las olas de

inundación • Bahías, estuarios y arroyos afectados por las

mareas • Liberaciones de energía hidroeléctrica en embalses

5­1
Machine Translated by Google

Capítulo 5 TD­41

• Inundaciones

repentinas • Inversiones del flujo tributario debido a

remansos • Inundaciones naturales por eventos de escorrentía de lluvia

Supuestos de flujo constante

El flujo constante supone que un caudal determinado persistió durante un tiempo suficientemente largo, por lo que una
suposición de flujo constante es válida. En general, esta suposición es cierta para tramos más cortos de un sistema fluvial y
para eventos en los que la superficie del agua sube y baja lentamente. Sin embargo, a medida que el dominio del modelado
se hace más grande y/o los caudales aumentan y disminuyen rápidamente, los caudales máximos no ocurren al mismo
tiempo espacialmente. Para estas condiciones, la suposición de un caudal constante comienza a fallar y puede no ser
apropiada ni siquiera para calcular las elevaciones de la superficie del agua.

La pendiente de la corriente también es un factor muy importante al asumir el supuesto de flujo constante.
Para arroyos de pendiente media a pronunciada, la superficie del agua calculada se basa en el terreno, la rugosidad
y los caudales en las inmediaciones de donde se calcula la superficie del agua (excepto cuando hay un remanso
significativo de una estructura o constricción aguas abajo). Por lo tanto, calcular la elevación de la superficie del
agua basándose en los caudales máximos en todos los lugares es una suposición válida incluso para sistemas muy
grandes, en los que el flujo máximo no se produjo simultáneamente.

En la Figura 5­1 se muestran los perfiles de la superficie del agua para una corriente moderadamente empinada
(5 pies/milla), calculados con un modelo de flujo constante. Ambos perfiles tienen un caudal de 9.000 pies cúbicos por
segundo (cfs) aguas arriba. Uno de los perfiles trazados se calculó con los flujos máximos ingresados en todas las
ubicaciones (Figura 5­1). El segundo perfil se basa en los caudales instantáneos en el momento en que el flujo máximo
(9,000 cfs) estaba en el extremo aguas arriba del sistema. Observe que debido a que la pendiente es relativamente
pronunciada, la superficie de agua resultante en el extremo aguas arriba (donde el flujo es de 9,000 cfs para ambos
tramos) es la misma (Figura 5­1).

Sin embargo, a medida que la pendiente del arroyo se aplana, las elevaciones/tasas de flujo de la superficie del agua
aguas abajo afectarán el cálculo de las elevaciones de la superficie del agua aguas arriba. Cuanto más plana sea la
pendiente, mayor será la distancia aguas arriba que se verá afectada por las superficies de agua aguas abajo. Para este
tipo de situaciones, la suposición de un flujo constante produciría elevaciones de la superficie del agua demasiado altas.
La sobreestimación de la elevación se debe al hecho de que una superficie de agua calculada aguas arriba, basada
en un caudal máximo, estará sesgada por las superficies de agua aguas abajo, que también se calcula en función de
caudales máximos. Sin embargo, si esos caudales máximos no ocurrieron al mismo tiempo, la suposición de flujo constante
no es válida y conducirá a una sobreestimación de las elevaciones de la superficie del agua.

En la Figura 5­2 se muestra un gráfico del perfil de la superficie del agua para un arroyo plano (0,5 pies/milla). Ambos
perfiles tienen un flujo de 9,000 cfs aguas arriba (Figura 5­2). Uno de los perfiles trazados se calculó con los flujos máximos
ingresados en todas las ubicaciones (Figura 5­2). El segundo perfil se basa en los caudales instantáneos en el
momento en que el flujo máximo (9,000 cfs) estaba en el extremo aguas arriba del sistema (Figura 5­2). Observe que
debido a que la pendiente es muy plana, la superficie de agua resultante en el extremo aguas arriba no es la misma, aunque
se utilizan 9,000 cfs para ambos perfiles en el extremo aguas arriba. Hay más de 0,5 pies de diferencia en la
superficie del agua resultante en la zona aguas arriba.

5­2
Machine Translated by Google

TD­41 Capítulo 5

Al final, el modelo de flujo constante de flujos máximos simultáneos en todas partes da respuestas más
altas que el modelo de flujo instantáneo (Figura 5­2).

Figura 5­1. Corriente empinada (pendiente = 5 pies/milla) con perfiles calculados utilizando flujos máximos
y flujos instantáneos.

Figura 5­2. Corriente plana (pendiente = 0,5 pies/milla) con perfiles calculados utilizando flujos máximos y
flujos instantáneos.

5­3
Machine Translated by Google

Capítulo 5 TD­41

Enrutamiento hidrológico versus hidráulico

Una aplicación exitosa de cualquier modelo de flujo estacionario requiere que los caudales ya hayan sido
calculados con precisión mediante un modelo hidrológico o medidos mediante un conjunto completo y preciso
de medidores de caudal (o algún otro método apropiado). El ruteo hidrológico consiste en resolver la ecuación
de continuidad y una relación entre el almacenamiento en el río y el vertido en la desembocadura del tramo del
ruteo. Algunos ejemplos de métodos de enrutamiento hidrológico son Modified Puls, Muskingum y
Muskingum­Cunge.

Si se utiliza un modelo hidrológico no sólo para calcular la precipitación­escorrentía sobre la cuenca, sino
también para realizar todo el recorrido dentro del sistema, entonces los caudales utilizados en el modelo de
flujo constante son tan precisos como el modelo hidrológico. Por lo tanto, el uso de un modelo hidráulico de
flujo estacionario se basa en el hecho de que se consideró que un modelo hidrológico era apropiado no
sólo para desarrollar el caudal a partir de los cálculos de precipitación­escorrentía, sino también para enrutar
todos los flujos a través del sistema durante el período. evento. Por lo tanto, una gran parte de la decisión entre
el modelado hidráulico de flujo estacionario versus el de flujo inestable se reduce a la pregunta: ¿es el
enrutamiento del flujo hidrológico lo suficientemente preciso como para producir caudales que puedan usarse en
los modelos hidráulicos de flujo estacionario correspondientes?

El enrutamiento hidráulico (enrutamiento de flujo inestable) resuelve las ecuaciones de continuidad y

momento juntas para enrutar los hidrogramas y calcular las elevaciones de la superficie del agua.

Diferencias computacionales

Para comprender mejor las diferencias entre el modelado de flujo estacionario y el modelado de flujo no
estacionario, el modelador debe conocer todas las diferencias computacionales entre los dos enfoques. La
siguiente es una descripción de las principales diferencias computacionales entre el enrutamiento de flujo 1D
estable y 1D inestable.

Las ecuaciones de flujo inestable (enrutamiento hidráulico) tienen una base más física porque se derivan de
la ecuación de continuidad y la segunda ley del movimiento de Newton:

∑=

dónde:
F = Suma de todas las fuerzas que actúan sobre un cuerpo de
agua m = Masa del cuerpo de agua a
= Aceleración (o desaceleración) del fluido

Ecuaciones de flujo estacionario e inestable

Como se mencionó anteriormente, las ecuaciones de flujo inestable se derivan de la segunda ley del movimiento
de Newton. En la Figura 5­3 se muestra un diagrama de las fuerzas que actúan sobre un cuerpo de agua en
una dimensión (es decir, flujo unidimensional).

5­4
Machine Translated by Google

TD­41 Capítulo 5

Figura 5­3. Fuerzas que actúan sobre un cuerpo de agua desde la sección transversal 2 hasta la sección transversal 1.

Aplicando la segunda ley del movimiento de Newton a un cuerpo de agua encerrado por dos secciones
transversales en las ubicaciones 1 y 2, se puede escribir la expresión para el cambio de impulso en una unidad de tiempo.
como:

−+−=Δ

dónde:
P = Fuerza debida a la presión hidrostática
Wx = Fuerza debida al peso del agua en dirección X
Ff = Fuerza debida a la fricción en el límite externo de 2 a 1
Q = Descarga ρ
= Densidad del agua
ΔVx = Cambio en la velocidad de 2 a 1 en la dirección X

La ecuación de momento unidimensional y la ecuación de continuidad se pueden escribir en forma de ecuación

diferencial parcial, con respecto a la descarga (Q), el área (A) y la profundidad (h), y comúnmente se
muestran en los libros de texto de hidráulica de la siguiente manera:

Ecuación del momento:

( )
+ + −+=0

Ecuación de continuidad:

+=

5­5
Machine Translated by Google

Capítulo 5 TD­41

dónde:
Q = Descarga β =
Coeficiente de distribución de velocidades
A = Área de la sección
transversal t
= Tiempo x = Distancia en la dirección del
flujo h = Profundidad del agua
S0 = Pendiente del lecho
Sf = Pendiente de fricción, de la ecuación de Manning
ql = Entradas laterales

Para un flujo estable, los términos basados en el tiempo en las ecuaciones de impulso y continuidad llegan a cero.
Por lo tanto, la forma de flujo estacionario de las ecuaciones unidimensionales de momento y continuidad se puede
escribir de la siguiente manera:

Forma de flujo estacionario de la ecuación del momento:

( /)
+ −+=0

Forma de flujo estacionario de la ecuación de continuidad:

Si bien la forma anterior de las ecuaciones de continuidad y momento se puede utilizar para resolver el flujo
estacionario unidimensional, en general, la mayoría de los programas de flujo estacionario 1D resuelven la
ecuación de energía unidimensional. La ecuación de energía unidimensional de flujo estacionario (a
menudo llamada ecuación de Bernoulli) se escribe como:

++ =++ +
2 2

dónde:
Z = Elevación del canal principal se invierte en las secciones transversales 1 y
2 Y = Profundidad del agua en las secciones
transversales 1 y 2 V = Velocidad promedio
del agua (Q/A) α = Coeficientes de
ponderación de velocidad g =
Aceleración gravitacional he = Pérdidas de energía de sección transversal 2 a 1.
(pérdidas por fricción (hf) y pérdidas
por contracción/expansión (hce)) hf = Pérdidas por fricción hf = LSf
L = distancia promedio ponderada entre secciones transversales
−
hce = Pérdidas por contracción y expansión =C
= Coeficiente de pérdida por contracción o expansión.

5­6
Machine Translated by Google

TD­41 Capítulo 5

Para flujo estacionario unidimensional, la ecuación de energía se resuelve iterativamente de una sección
transversal a la siguiente. Para flujo inestable, las ecuaciones de continuidad y momento se resuelven
simultáneamente, generalmente en un esquema de solución matricial que resuelve todo el espacio
(todos los puntos computacionales) en cada paso de tiempo (esquema de solución implícito). Sin embargo,
existen otros esquemas de solución que resuelven una sección transversal a la vez (esquemas de solución explícitos).

Además, debido a que existen derivadas basadas en el tiempo en las ecuaciones de flujo inestable, el
modelador debe seleccionar un paso de tiempo computacional apropiado para resolver las ecuaciones. El
paso de tiempo computacional se selecciona en función de la resolución de los hidrogramas a trazar, así como de
la precisión numérica y la estabilidad de la resolución de las ecuaciones matemáticas no lineales. Un enfoque
común para seleccionar el intervalo de cálculo es utilizar un criterio numérico de precisión/estabilidad llamado
condición de Courant:

∆
= ≤ 1,0
∆
Por lo tanto:
∆
∆≤

dónde:
C = Número de Courant
Vw = Velocidad de la onda de inundación (celeridad de la
onda) (pies/s) ΔT = Paso (s) de tiempo
computacional ΔX = Distancia promedio entre secciones transversales (1D) o tamaño de celda/elemento de cálculo
(2D y 3D)

Propiedades hidráulicas

Específicamente para el modelado 1D, al resolver la ecuación de energía 1D, todas las propiedades
hidráulicas (área; perímetro mojado, transporte, almacenamiento, etc.) se resuelven exactamente para cada
sección transversal según sea necesario. Debido a que las simulaciones de flujo inestable requieren
mucho más cálculo (las ecuaciones a menudo se resuelven de forma iterativa durante miles de pasos de tiempo),
las propiedades hidráulicas a menudo se calculan previamente para todas las elevaciones posibles de la superficie
del agua en cada sección transversal o puente/alcantarilla. Luego, las propiedades hidráulicas se interpolan a
partir de las curvas precalculadas durante los cálculos del flujo inestable. Generalmente se utilizan métodos de
interpolación lineal, por lo que existe cierto error dependiendo del número de curvas y del número de puntos en cada curva.
En la Figura 5­4 se muestra un ejemplo de una familia de curvas de elevación de flujo versus altura de agua
que se calcula previamente para un cruce de puente típico.

5­7
Machine Translated by Google

Capítulo 5 TD­41

Figura 5­4. Familia de ejemplo de curvas de calificación precalculadas para un puente.

Pérdidas por fricción

Para la ecuación de energía 1D, el término hf mide la energía interna disipada en toda la masa de agua entre
las dos secciones transversales, mientras que el elemento hf en la ecuación de momento mide las pérdidas
debidas a fuerzas externas ejercidas sobre el agua por las paredes del canal. La distinción inherente entre los dos
principios radica en el hecho de que la energía es una cantidad escalar mientras que el momento es una
cantidad vectorial. Haciendo caso omiso de las pequeñas diferencias en los coeficientes de ponderación
de velocidad α y β, para un flujo gradualmente variado, las pérdidas de energía interna, calculadas a partir de la
ecuación de energía, son prácticamente idénticas a las fuerzas externas en la ecuación de momento (Chow,
VT, Open Channel Hydraulics, McGraw ­Hill, 1959, pág.51.).

Tanto el flujo estacionario como el flujo inestable utilizan ecuaciones de pérdida por fricción, como la ecuación de
Manning, para describir las pérdidas de energía internas y las fuerzas externas debidas a la fricción. En
ambos métodos, la ecuación de Manning se utiliza para calcular el término de pendiente de fricción Sf en un punto
(es decir, sección transversal o cara de celda 2D). Además, tanto el flujo estacionario como el flujo inestable
requieren que se calcule una pendiente de fricción promedio entre las dos secciones transversales, para calcular
con precisión la pérdida por fricción a lo largo de la longitud entre las secciones transversales. Como
existen diferentes formas de calcular una pendiente de fricción promedio, esto puede ser una fuente de
diferencias entre los dos enfoques computacionales. Por ejemplo, HEC­RAS utiliza un método llamado
ecuación de "Transporte promedio" para calcular la pendiente de fricción promedio para un flujo estacionario
1D. Sin embargo, la pendiente de fricción promedio para flujo inestable se calcula con la ecuación de
“Pendiente de fricción promedio” en HEC­RAS (Manual de referencia hidráulica de HEC­RAS, Capítulo 2, septiembre de 2016).

5­8
Machine Translated by Google

TD­41 Capítulo 5

Pérdidas por contracción y expansión

El enfoque del momento integra fuerzas que actúan sobre las superficies y extremos de un volumen de
control; por lo tanto, los impactos de las contracciones/expansiones del flujo se capturan en las fuerzas en
los extremos aguas arriba y aguas abajo de ese volumen de control. Por supuesto, para que este enfoque
funcione correctamente, se necesita una selección adecuada de las áreas de flujo mediante una
contracción y expansión (ubicación de la sección transversal). Por otro lado, el enfoque energético integra
trabajo/energía para el volumen de control; Los coeficientes empíricos multiplicados por el cambio en la
carga de velocidad se utilizan para describir las pérdidas asociadas con el gasto de energía turbulenta
asociado con la contracción/expansión del flujo. Sin embargo, la investigación ha descubierto que el uso
de un modelo calibrado para flujo constante dentro del solucionador de flujo inestable HEC­RAS puede
dar como resultado superficies de agua calculadas más bajas debido a que se omiten las pérdidas completas
por la turbulencia de contracción/expansión. Debido a que esta es una diferencia computacional conocida en
HEC­RAS entre flujo estacionario y flujo inestable, las pérdidas empíricas por contracción y expansión
generalmente se agregan a los algoritmos computacionales de flujo inestable como una opción. Para ilustrar
esta opción, la Figura 5­5 muestra los resultados del modelo HEC­RAS de flujo constante calibrado;
resultados del modelo de flujo inestable con exactamente la misma geometría y datos de flujo; y luego un
modelo final de flujo inestable calibrado. El modelo de flujo inestable se calibró activando las fuerzas
empíricas de contracción/expansión y ajustando los coeficientes para que coincidan con el modelo de flujo estable ya calibrado pa

Figura 5­5. Modelo de flujo estacionario calibrado y modelo de flujo inestable con y sin pérdidas por contracción
y expansión agregadas.

5­9
Machine Translated by Google

Capítulo 5 TD­41

Áreas de almacenamiento/flujo ineficaz

Para el modelado de flujo unidimensional (1D), es necesario definir partes de las secciones transversales como áreas
de flujo ineficaces para obtener la cantidad correcta de área de flujo activa (efectiva). Si no se cumple este requisito,
el área de flujo será incorrecta en muchos lugares, lo que producirá elevaciones y velocidades de la superficie del
agua demasiado bajas o demasiado altas, lo que a su vez afectará el cálculo de las pérdidas por fricción y las
pérdidas por contracción/expansión. Para el modelado de flujo constante, las áreas de flujo ineficaz en realidad solo se
usan para describir qué partes de la sección transversal tienen agua en movimiento y cuáles no. La contabilidad del
volumen no se realiza en flujo estacionario, por lo que los efectos que el almacenamiento de la llanura aluvial
tiene en el hidrograma se realizan fuera del modelo hidráulico (enrutamiento hidrológico). En la Figura 5­6 se muestra
un ejemplo de diseño de áreas de flujo ineficaces para el modelado de flujo estable 1D y flujo inestable 1D.

Figura 5­6. Ejemplo de diseño de áreas de flujo ineficaces (polígonos rellenos con líneas diagonales negras) para 1D
modelado.

Para el modelado de flujo inestable 1D, las áreas de flujo ineficaz no solo se utilizan para definir el área de
flujo ineficaz, sino que también se utilizan para calcular los volúmenes de almacenamiento entre las secciones
transversales. A medida que un hidrograma se encamina a través de un tramo que contiene áreas de flujo
ineficaces (y por lo tanto áreas de volumen de almacenamiento), el agua saldrá del canal para llenar estos
volúmenes de almacenamiento. A medida que pase la ola de inundación, el agua volverá de las áreas de almacenamiento al

5­10
Machine Translated by Google

TD­41 Capítulo 5

canal como flujo en la rama descendente del hidrograma. En la Figura 5­7 se muestra un ejemplo del efecto del
almacenamiento en sección transversal debido a áreas de flujo ineficaces.

En la Figura 5­7 hay dos hidrogramas, uno en el extremo aguas arriba del tramo y otro en el extremo aguas abajo.
Este ejemplo de tramo de río contiene una cantidad significativa de áreas de flujo ineficaz. A medida que el hidrograma
pasa por el tramo, el agua ingresa al área de almacenamiento. El agua en el lado ascendente del hidrograma y el pico
fluye hacia el área de almacenamiento. Luego, a medida que pasa el pico, el agua sale del almacenamiento y
agrega flujo a la rama descendente del hidrograma.

FLUJO SUPERIOR DEL BARCO DE VAPORCR

3500
Leyenda

17943.58

AFLUENCIA
3000

2500

2000
sLfCcF(
O)JS)UF

1500

1000

500
2400 0600 1200 1800 2400 0600 1200 1800 2400
01 de enero de 1997 02 de enero de 1997
Tiempo

Figura 5­7. Hidrograma que entra y sale del tramo de un río con áreas de flujo ineficaces que actúan como
almacenamiento.

Estrategias de calibración

Las estrategias de calibración de modelos varían entre el modelado de flujo estacionario y el modelado de flujo inestable.
En general, se requieren diferentes enfoques debido al hecho de que hay menos información calculada para
comparar con los datos observados cuando se realiza un enfoque de modelado de flujo estacionario.

En general, la estrategia de calibración para el modelado de flujo constante 1D es comparar las elevaciones calculadas
de la superficie del agua con los datos medidos y cualquier marca de agua alta disponible. Generalmente se
realizan ajustes a los valores n de Manning, a los coeficientes de contracción y expansión, y a cualquier
coeficiente de estructura hidráulica (coeficientes de vertedero, coeficientes de puentes y alcantarillas, etc.).
Es posible que se necesiten ajustes adicionales para la ubicación de estaciones/elevaciones de áreas
de flujo ineficaces y elevaciones de diques de sección transversal. Para el modelado 2D y 3D, se realizaría una
calibración adicional de los perfiles de velocidades si las mediciones de velocidad observadas están disponibles.

5­11
Machine Translated by Google

Capítulo 5 TD­41

Para el modelado de flujo inestable 1D, la estrategia de calibración comienza de la misma manera que para el modelado de flujo
estacionario; sin embargo, el modelador también compara hidrogramas completos (elevación de la superficie del agua y flujo versus
tiempo) con los resultados calculados. Entonces, además de observar las etapas máximas, se evalúa toda la forma, el tiempo, el
volumen y la magnitud de los hidrogramas y se comparan con los valores observados. El proceso de calibración es más complejo
con el modelado de flujo inestable.
Las diferencias en el momento y la magnitud podrían deberse a coeficientes de rugosidad incorrectos, definiciones deficientes del
almacenamiento en llanuras aluviales o incluso coeficientes de estructura hidráulica incorrectos. Los ajustes realizados al modelo
aguas arriba afectarán el tiempo y la magnitud del flujo/etapa aguas abajo.
Los ajustes realizados aguas abajo pueden afectar las elevaciones de la superficie del agua aguas arriba. Además, es muy
importante realizar un seguimiento del volumen dentro del sistema. A veces, las condiciones de contorno de flujo proporcionadas al
modelo pueden no describir con precisión el volumen que ingresa al sistema o el momento de los picos de inundación.

Limitaciones del modelado de flujo constante

Los modelos de flujo estacionario, o incluso ejecutar un modelo de flujo inestable (1D, 2D o 3D) en un modo de flujo estable (flujo
constante), generalmente no deben usarse cuando existen las siguientes situaciones en el sistema fluvial que se analiza (esto no es
un análisis exhaustivo). lista):

El río está influenciado por las mareas, y la marea tiene un efecto significativo en la superficie del agua.
elevaciones para el área de interés.

Los eventos que se modelan son muy dinámicos con respecto al tiempo (es decir, olas de inundación que rompen represas;
inundaciones repentinas; sistemas fluviales en los que el caudal máximo aumenta muy rápidamente, se mantiene alto por
un tiempo muy corto y luego retrocede rápidamente).

Durante el evento se producen redes de flujo complejas y/o inversiones de flujo.

Durante el evento se producen eventos dinámicos, como desbordamientos y rupturas de diques.

Sistemas fluviales extremadamente planos, donde la gravedad, la presión hidrostática y la fricción no son necesariamente
las únicas fuerzas significativas que actúan sobre el flujo (es decir, fuerzas de aceleración locales y convectivas).

Sistemas con Estaciones de Bombeo que mueven una cantidad importante de agua.

Sistemas con estructuras que tienen operaciones de compuerta complejas basadas en etapas y flujos en
el sistema.

Sistemas con una enorme cantidad de almacenamiento en los que el hidrograma se atenuará
significativamente.

Incluso considerando todo lo dicho anteriormente. Todavía hay muchas áreas en las que se puede utilizar un buen modelo
hidrológico (uno que sea representativo de la cuenca y que haya sido bien calibrado) junto con un modelo hidráulico de flujo constante
para realizar estudios de cuenca.

5­12
Machine Translated by Google

TD­41 Capítulo 5

que requieren resultados del modelo hidráulico. Sin embargo, corresponde al modelador decidir cuándo no es
apropiado utilizar un enfoque de modelado de flujo estacionario.

Limitaciones del modelado de flujo inestable

Condiciones iniciales.

Los modelos de flujo inestable generalmente requieren la especificación de estimaciones iniciales de la elevación y
velocidad de la superficie del agua (es decir, condiciones iniciales) en todo el dominio del modelo. Estas estimaciones
pueden obtenerse de diversas formas. Algunos modelos calcularán las condiciones iniciales a partir de un análisis de flujo
estacionario basado en los flujos iniciales y las etapas especificadas como condiciones límite en la simulación. El
modelador puede especificar entradas cero y una superficie de agua constante y luego drenar el modelo reduciendo las
etapas aguas abajo mientras se aumentan las entradas, un cálculo de flujo inestable, hasta que se calculen las
condiciones iniciales deseadas (es decir, un “arranque en frío”). Otra alternativa es obtener las condiciones iniciales
de una simulación de modelo anterior (es decir, un “arranque en caliente”) o de otro modelo (es decir, un “arranque
en caliente”).

En la mayoría de los casos, las condiciones iniciales serán una representación imperfecta del estado del sistema al
inicio de una simulación de flujo inestable. Estas imperfecciones pueden variar desde inconsistencias localizadas a corto
plazo en la velocidad y profundidad calculadas hasta errores significativos en el volumen de agua presente en el modelo al
comienzo de la simulación. Por lo tanto, los modelos de flujo inestable generalmente requieren un período de
"calentamiento" o "aceleración" al comienzo de la simulación para "eliminar" estas imperfecciones del sistema. En el
caso de un arranque en caliente que simplemente extiende la simulación de un conjunto dado de condiciones límite (es
decir, los hidrogramas de flujo y de etapa especificados en los límites del modelo), el período de aceleración puede ser
insignificante. Para dominios de modelos grandes y complejos o ríos relativamente planos, eliminar errores de
volumen significativos del sistema puede requerir un período de calentamiento del orden de tiempo necesario para
encaminar un hidrograma de inundación a través del sistema.

Cabe señalar que los modelos de transporte de calidad del agua generalmente requieren conjuntos adicionales de datos
sobre la condición inicial. Los modelos de transporte pueden ser más sensibles a las imperfecciones en las
condiciones iniciales que los modelos hidrodinámicos y requieren períodos de puesta en marcha más largos
para calcular condiciones iniciales aceptables para una simulación.

Selección de eventos.

Como se describió anteriormente, los modelos de flujo estacionario tienden a generar estimaciones conservadoras de
los perfiles de la superficie del agua en comparación con los modelos de flujo inestable. Asimismo, se puede generar
una estimación conservadora de las velocidades máximas del canal con un modelo de flujo estable reduciendo los
coeficientes de rugosidad hidráulica en el canal dentro de límites razonables. En el modelado de flujo inestable, la
selección de los hidrogramas de simulación (condiciones de contorno) utilizados para desarrollar estas estimaciones
puede tener una influencia significativa en los perfiles de velocidad o superficie máxima del agua calculados.

Considere el siguiente ejemplo. En una simulación de flujo inestable del paso de una inundación a través de un río
típico, las etapas calculadas para una descarga determinada serán mayores en el tramo descendente del hidrograma
que en el tramo ascendente del hidrograma (hay excepciones a este caso típico).

5­13
Machine Translated by Google

Capítulo 5 TD­41

Un modelo conceptual simple de este fenómeno de “curva de calificación en bucle” es que los niveles del
río aguas abajo de un lugar de observación serán más bajos en el tramo ascendente que en el tramo descendente.

Ahora, consideremos el caso en el que el hidrograma de inundación se compone de múltiples eventos muy
próximos entre sí con el mismo caudal máximo. Si la etapa aguas abajo calculada nunca recupera su condición
original, podemos esperar razonablemente que cada pico de descarga sucesivo en nuestro punto de
observación produzca una etapa pico más alta (y una velocidad pico más baja). Por lo tanto, una estimación
del nivel máximo basada en una simulación de flujo inestable de un hidrograma de inundación simple con un solo
pico podría subestimar el riesgo de inundación real en el prototipo. En los ríos aluviales, la “curva de clasificación
del circuito” observada puede ser mayor (bucle más ancho) que el “bucle hidrodinámico” predicho por un modelo
de flujo inestable. Esta diferencia normalmente se atribuye a cambios mal cuantificados en el transporte del
canal impulsados por procesos de sedimentación.

5­14
Machine Translated by Google

TD­41 Capítulo 6

Capítulo 6

Modelado unidimensional versus bidimensional

La cuestión del modelado hidráulico 1D versus 2D es mucho más difícil que la del flujo estacionario versus el
inestable. Definitivamente hay algunas áreas en las que el modelado 2D puede producir mejores resultados que
el modelado 1D, y también hay situaciones en las que el modelado 1D puede producir resultados que son
tan buenos como los modelos 2D con menos esfuerzo y requisitos computacionales.
Desafortunadamente, existe una gran variedad de situaciones que caen en un área gris, y se podrían enumerar
los aspectos positivos y negativos de ambas metodologías para aplicaciones específicas.

Este capítulo analiza las diferencias entre el modelado hidráulico unidimensional (1D) y bidimensional (2D), de
flujo estacionario y no estacionario. Específicamente, este capítulo describe la definición de modelado hidráulico
1D y 2D y los supuestos utilizados en cada uno; las ecuaciones utilizadas para el modelado de flujo inestable 1D
y 2D; diferencias computacionales; ejemplos de aplicaciones; modelo de desarrollo; calibración de modelos;
modelador, conocimientos, destrezas y habilidades; y un resumen de las ventajas y desventajas del
modelado 1D y 2D.

Definiciones

En general, casi todo el movimiento de fluidos es tridimensional. Sin embargo, las ecuaciones del movimiento
del agua a menudo se derivan en forma unidimensional y bidimensional, para una amplia gama de aplicaciones
prácticas.

Cuando las ecuaciones de movimiento se derivan en forma unidimensional, se supone que las fuerzas
que actúan sobre un cuerpo de agua son predominantes en una dirección, x, a lo largo de la línea central del canal
del río. Esto supone que las fuerzas verticales y laterales que actúan sobre el cuerpo de agua son pequeñas en
comparación con la dirección x y, por lo tanto, se puede suponer que son insignificantes.
Entonces, la forma unidimensional de las ecuaciones de movimiento solo tiene en cuenta las fuerzas en la
dirección x.

Cuando las ecuaciones se derivan en forma bidimensional, se asume que las fuerzas que actúan sobre una masa
de agua son predominantes en x (a lo largo de la línea central del canal del río) y en y (lateralmente a
través del canal o llanura aluvial). , que es la forma más común para el modelado 2D, o en x y z (verticalmente
a lo largo de la profundidad). Para la forma de las ecuaciones escritas en términos de x e y, el promedio vertical
(z, o fuerzas verticales) supone que las fuerzas verticales que actúan sobre un cuerpo son pequeñas en
comparación con las fuerzas horizontales (x, y) y, por lo tanto, se consideran ser insignificante. Para la forma de
las ecuaciones escritas en términos de x y z, el promedio lateral (y, o fuerzas laterales) supone que las fuerzas
laterales son insignificantes.

6­1
Machine Translated by Google

Capítulo 6 TD­41

Ecuaciones unidimensionales y bidimensionales

Las leyes físicas que gobiernan el flujo del agua en una corriente son: (1) el principio de
conservación de la masa (continuidad) y (2) el principio de conservación del momento. Estas
leyes se expresan matemáticamente en forma de ecuaciones diferenciales parciales, a las
que en adelante nos referiremos como ecuaciones de continuidad y de momento. En la
derivación de las ecuaciones de flujo inestable unidimensional (1D) y bidimensional (2D) (a
menudo llamadas ecuaciones de aguas poco profundas), se hacen varias suposiciones sobre el flujo:

1. El agua es un fluido incompresible.

2. La distribución de presiones se considera hidrostática.
3. La aceleración vertical del agua se considera despreciable.
4. La pendiente del lecho se considera suave (menos de 1:10 de pendiente).
5. Los efectos de la fricción límite se pueden tener en cuenta con leyes de resistencia al flujo derivadas
para flujos estacionarios (es decir, la ecuación de Manning).
6. La aproximación de Boussinesq es válida (ignorando las fuerzas causadas por diferencias en el agua
densidad).

Ecuaciones unidimensionales (1D)

Las ecuaciones unidimensionales de continuidad y momento se pueden escribir en forma de ecuación diferencial
parcial, con respecto a la profundidad (h) y la velocidad (u), y comúnmente se muestran en los libros de
texto de hidráulica de la siguiente manera:

Ecuación de continuidad:
( )
+ −=0

Ecuación del momento:

++ −++=0

dónde:
tu = velocidad en la dirección x =
h profundidad del
gt agua =

X gravedad = tiempo = distancia en la dirección del flujo (plano x)

q = Término de entrada lateral (fuente/sumidero)
S0 = Pendiente del
sf lecho = Pendiente de fricción, de la ecuación de Manning.
sh = término de fuerza agregada (pérdidas menores adicionales)

La Figura 6­1 proporciona un diagrama que define los símbolos utilizados en las ecuaciones de movimiento
1D y 2D.

6­2
Machine Translated by Google

TD­41 Capítulo 6

Figura 6­1. Definición de símbolos utilizados en las ecuaciones de movimiento 1D y 2D.

Ecuaciones bidimensionales (2D)

La ecuación bidimensional de continuidad y momento se puede escribir en forma de ecuación diferencial

parcial, con respecto a la profundidad (h) y la velocidad (u, v, U, V), y comúnmente se muestran en los libros de
texto de hidráulica de la siguiente manera:

Ecuación de continuidad promediada verticalmente:

( ) ( )
+ + −=0

Ecuaciones de momento promediadas verticalmente:

+ +=−+ + −+

++ =­ + + −−

dónde:
v = velocidad en la dirección
y y = distancia en la dirección lateral (plano y)
H = elevación de la superficie del agua (z + profundidad)
Vermont
= coeficiente de viscosidad del remolino
cf horizontal = coeficiente de fricción del fondo
F = parámetro de Coriolis

6­3
Machine Translated by Google

Capítulo 6 TD­41

Ecuación de continuidad promediada lateralmente:

+ −=0

Ecuaciones de momento promediadas lateralmente:

UB UB B
Ud. BH PA   g 1 B 1 B
UW   x 
xx
 xz

t z X X X z

1PH


z z
gramo

dónde:
v = velocidad en la dirección
y y = distancia en la dirección lateral (plano y)
H = elevación de la superficie del agua (z +
Vermont
Profundidad) = coeficiente de viscosidad del
cf remolino horizontal = coeficiente
de fricción del fondo f = parámetro de Coriolis
U = velocidad lateralmente promediada en la dirección x
W = velocidad lateralmente promediada en dirección
B z = ancho
P = presión lateral promediada
xx xz = tensiones turbulentas en las direcciones xx y xz, respectivamente = flujo
q de entrada lateral por unidad de volumen

Como se muestra arriba, las ecuaciones de aguas poco profundas 1D y 2D son muy similares, con las
principales diferencias: las ecuaciones 1D solo se derivan para fuerzas que actúan en la dirección x (Figura
6­1), mientras que las ecuaciones 2D tienen en cuenta las fuerzas que actúan en ya sea la dirección xey, o
la dirección xeyz (Figura 6­1); las ecuaciones 1D utilizan un término de fuerza adicional para describir
fuerzas adicionales debido a contracciones y expansiones severas (basadas en un coeficiente empírico),
mientras que las ecuaciones 2D resuelven esto directamente con la inclusión de los términos de viscosidad
de Foucault/esfuerzo turbulento; y las ecuaciones 2D tienen un término agregado para tener en
cuenta la rotación de la Tierra (Coriolis), que no se puede incluir en el enfoque 1D. En la práctica general,
los efectos de Coriolis se ignoran en las ecuaciones 2D promediadas lateralmente, a menos que esté
modelando áreas más cercanas a los polos norte y sur, donde el efecto de Coriolis es más predominante.

Los términos de las ecuaciones presentadas se pueden modificar para tener en cuenta las diferencias de
densidad, o se pueden agregar términos adicionales para otras fuerzas, como el viento y los flujos de lodo/
escombros. Tenga en cuenta que la presencia de diferencias de densidad no viola la aproximación de Bousinnesq.

Forma de onda de difusión de la ecuación del momento

Se pueden realizar simplificaciones de las formas 1D y 2D de las ecuaciones, siendo la más común la
forma de onda de difusión de las ecuaciones. La aproximación de onda de difusión se utiliza a menudo
en formas de solución tanto 1D como 2D. Para llegar a la forma de difusión de las ecuaciones, simplemente

6­4
Machine Translated by Google

TD­41 Capítulo 6

elimine los términos de aceleración (cambios en la velocidad con respecto al tiempo y el espacio) en las ecuaciones
de momento. Esta forma de ecuación de momento solo contiene fuerzas de gravedad, fricción y presión hidrostática.

La forma de difusión de la ecuación del momento se puede combinar con la ecuación de continuidad y escribirse en términos
de resolver únicamente la elevación de la superficie del agua. Esto hace que la ecuación sea más fácil de resolver,
generalmente más estable para una amplia gama de problemas y requerirá mucho menos tiempo de cálculo. Sin embargo, sin
los términos de aceleración en las ecuaciones, las ecuaciones de ondas de difusión son menos precisas que las ecuaciones
completas y son menos aplicables a toda la gama de problemas que el modelador puede necesitar resolver. Esto se discutirá
más adelante en este documento.

De ahora en adelante, se asumirá que 2D implica ecuaciones de aguas poco profundas promediadas verticalmente en 2D, a
menos que se especifique lo contrario.

Diferencias computacionales

Para comprender mejor las diferencias entre el modelado de flujo inestable 1D y 2D, el modelador debe conocer todas las
diferencias computacionales entre los dos enfoques.
La siguiente es una descripción de las principales diferencias computacionales entre el enrutamiento de flujo inestable 1D
y 2D.

Superficie del agua y velocidades

La diferencia más obvia entre los modelos 1D y 2D es que los modelos 1D solo calculan una única elevación de la superficie
del agua en cada sección transversal, mientras que los modelos 2D calculan una superficie de agua única para cada celda/
cara del modelo (ver ejemplo en la Figura 6­2). . Además, las velocidades en un modelo 1D se calculan como velocidades
promedio (vertical y horizontalmente) para el canal principal, el banco izquierdo y el banco derecho dentro de cada sección
transversal (ver ejemplo en la Figura 6­
3). Algunos modelos 1D, como HEC­RAS, calculan velocidades 1D entre los puntos de ruptura del valor n de Manning para el
área overbank y, además, pueden publicar los resultados del proceso en velocidades más espaciales. Sin embargo,
esas velocidades se basan en el supuesto de que el flujo es perpendicular a la sección transversal y que la
distribución del flujo es sólo una función de la distribución de transporte de la sección transversal.

Los modelos 2D calculan velocidades promedio (vertical y horizontalmente) para cada cara de celda.
Sin embargo, el detalle de las superficies del agua y las velocidades en un modelo 2D depende del número de celdas que
utiliza un modelador a lo largo del canal principal y las llanuras aluviales (ver ejemplo en la Figura 6­4).

6­5
Machine Translated by Google

Capítulo 6 TD­41

Figura 6­2. Ejemplo de gráfico de superficie de agua para una solución de modelo 1D y 2D.

Figura 6­3. Ejemplo de salida de velocidad de un modelo 1D.

6­6
Machine Translated by Google

TD­41 Capítulo 6

Figura 6­4. Ejemplo de salida de velocidad en una sección transversal de un modelo 2D.

Pérdidas por fricción

Las pérdidas por fricción se calculan para modelos 1D multiplicando una pendiente de fricción promedio (Sf )
por la longitud entre las secciones transversales. La pendiente de fricción se calcula en cada sección
transversal con la ecuación de Manning. Se pueden utilizar diferentes técnicas de promedio para
calcular la pendiente de fricción promedio entre secciones transversales. Además, la distancia entre las
secciones transversales se pondera según el flujo en la margen izquierda, el canal principal y la margen
derecha, y sus correspondientes longitudes de alcance. Los modelos 2D también calculan una pendiente de
fricción en cada cara de las celdas; sin embargo, la pendiente de fricción generalmente no se promedia sobre
la celda, ya que la dirección del flujo es bidimensional. La pendiente de fricción se utiliza para calcular una
fuerza de fricción en cada cara y se incluye como un componente en la solución de la ecuación de momento bidimensional.

Cálculos de transporte

Para los modelos 1D, el transporte se calcula para el canal principal de la sección transversal como un área
de flujo separada. Esto significa que el área completa, el perímetro mojado y la rugosidad se utilizan para
calcular un transporte único para el canal principal en cada elevación de la superficie del agua. El transporte
en los márgenes también se divide en áreas de flujo separadas en función de las irregularidades. Por
lo tanto, un área de sobrebanco se tratará como un área de transporte separada si tiene un único valor n de
Manning, o puede tratarse como varias áreas de transporte si tiene múltiples valores n de Manning. En
la Figura 6­5 se muestra un ejemplo de una sección transversal 1D dividida en áreas de transporte separadas.
Como se muestra, el margen izquierdo se divide en dos áreas de transporte, aunque tienen la

6­7
Machine Translated by Google

Capítulo 6 TD­41

mismo valor n de Manning. El transporte se divide de esta manera porque la ecuación de Manning se desarrolló
para áreas de flujo separadas (originalmente solo canales). Dada la variabilidad de la profundidad y la
velocidad, y la relación entre el área y el perímetro mojado en la ecuación de Manning, es más apropiado
tratar la margen izquierda como dos áreas de flujo separadas.

Figura 6­5. Ejemplo de áreas de transporte separadas para una sección transversal 1D.

Los modelos 2D calculan el transporte por separado para cada cara de cada celda. Esto puede generar
diferencias en los resultados entre los modelos 1D y 2D, incluso solo para el flujo del canal principal. Por
ejemplo, digamos que tenemos un canal trapezoidal simple. Un modelo 1D tratará el canal como un área de
flujo única y calculará el transporte utilizando el área total y el perímetro mojado para cada elevación de la
superficie del agua. Para un modelo 2D dependerá de cuántas celdas se estén utilizando para modelar el
canal principal. Si se utiliza una sola celda, entonces se utilizará toda el área y el perímetro mojado para calcular
el transporte, como en un modelo 1D. Sin embargo, si se utilizan varias celdas para modelar el canal
principal, entonces el transporte se calcula por separado para cada cara que cruza el canal. Este enfoque
producirá una cantidad diferente de transporte para una superficie de agua determinada. En general, dividir
el canal en pedazos producirá una mayor cantidad de transporte para una superficie de agua determinada.
El efecto neto de esto es que cuando un canal se divide en varios pedazos para calcular el transporte, la superficie
del agua calculada será menor debido al aumento percibido en la eficiencia del transporte del canal más alto.

En la Figura 6­6 se muestra un ejemplo de canal trapezoidal. Este canal fue modelado con secciones
transversales 1D y también con celdas 2D. Para el modelado 2D, se crearon tres modelos 2D diferentes
(usando HEC­RAS). Se creó un modelo 2D como un modelo de una sola celda en todo el canal, otro con 4
celdas en todo el canal (como se muestra en la Figura 6­6) y un tercero con 10 celdas en todo el canal. Todos
los modelos utilizaron la misma condición límite aguas abajo, que era una curva de calificación de valor único.
Los resultados de este experimento se muestran en la Figura 6­7. El modelo 1D y el modelo 2D, que utilizan
una sola celda para el canal, dan básicamente los mismos resultados (las elevaciones más altas de la superficie
del agua calculadas). Los modelos 2D que utilizaron 4 y 10 celdas a lo largo del canal también dieron
aproximadamente los mismos resultados, pero inferiores que el modelo 1D y el enfoque de una sola celda.
Esto se debe al hecho de que si se divide el canal en pedazos y se calcula por separado el transporte
para esas piezas, no se obtiene el mismo transporte que cuando se calcula el transporte como un área de flujo
única. El transporte es altamente no lineal con respecto a la relación entre el área y el perímetro mojado. Las
mayores diferencias ocurrirán cuando haya bancos empinados que se computan por separado de la
porción principal del canal.
En el ejemplo que se muestra en la Figura 6­6 (4 celdas), el transporte que se calcula para las celdas que
representan los bancos es muy bajo, debido a la gran cantidad de relación entre el perímetro mojado y el área de flujo.

6­8
Machine Translated by Google

TD­41 Capítulo 6

Las dos celdas en el medio tienen una relación entre el perímetro húmedo y el área de flujo mucho menor
y, por lo tanto, un transporte mucho mayor. El transporte general para todo el canal termina siendo mayor
que si lo calculara como un área de flujo única, y el resultado neto es una superficie de agua más baja.

Figura 6­6. Ejemplo de canal trapezoidal para cálculos de transporte 1D y 2D.

Figura 6­7. Resultados de la computación de transporte con modelos 1D y 2D de diferentes tamaños de celda.

6­9
Machine Translated by Google

Capítulo 6 TD­41

Físicamente, lo que realmente sucederá es que el agua más lenta en el área de las orillas ralentizará el agua
que va hacia el medio del canal, y el agua más rápida en el medio del canal afectará la velocidad del agua sobre
las orillas del canal. . Este fenómeno se maneja en modelos 2D mediante el uso de modelado de turbulencia. La
fuerza de corte agua­agua entre dos celdas adyacentes se puede explicar directamente en modelos 2D con
modelado de turbulencia. En la Figura 6­8 se muestra el mismo experimento; sin embargo, el modelado
de turbulencia se activó para todos los enfoques de modelado 2D. Como se muestra, ahora todos los modelos,
1D, 2D de una sola celda, 2D de cuatro celdas y 2D de diez celdas, obtienen aproximadamente las mismas
elevaciones de la superficie del agua.
Sin embargo, se tuvo que ingresar un coeficiente de turbulencia especificado por el usuario para los modelos
2D. Este coeficiente fue calibrado para dar los mismos resultados que el modelo 1D.

Figura 6­8. Resultados de ejemplo para un canal trapezoidal modelado en 1D y 2D con el modelado
de turbulencia activado.

Pérdidas por contracción y expansión

Las contracciones y expansiones del flujo son inherentemente un fenómeno de flujo tridimensional. Los modelos 1D
manejan las pérdidas de fuerza/energía en contracciones y expansiones utilizando coeficientes empíricos
multiplicados por un cambio en la altura de velocidad (modelos 1D de flujo estacionario e inestable).
Además, los modelos de flujo inestable 1D calculan fuerzas de aceleración espacial y fuerzas diferenciales
de presión a través de la contracción y expansión del flujo. Sin embargo, las fuerzas de aceleración y presión
espacial sólo se calculan en una dirección.

6­10
Machine Translated by Google

TD­41 Capítulo 6

Los modelos 2D capturan las fuerzas debidas a contracciones y expansiones con fuerzas de presión
bidimensionales y términos de aceleración espacial. Se capturan fuerzas adicionales con modelos de
turbulencia para tener en cuenta todas las pérdidas de fuerza/energía asociadas con una
contracción/expansión de flujo. Como se mencionó anteriormente, el modelado de turbulencia requiere
que el usuario ingrese y calibre coeficientes de turbulencia.

Áreas de flujo ineficaces

Los modelos 1D utilizarán toda la sección transversal para representar el agua en movimiento activo.
Sin embargo, partes de una sección transversal a menudo tendrán agua que no se transporta (sin
velocidad) en la dirección aguas abajo. Debido a las áreas sin transporte, el modelado 1D requiere el
uso de áreas de "flujo ineficaz" para modelar las partes de la sección transversal que están húmedas pero que
no producen ningún transporte aguas abajo. Esto requiere mucho criterio por parte del usuario para poder
colocar correctamente estas áreas de flujo ineficaz dentro de la sección transversal. Además, es posible
que sea necesario cerrar las áreas de flujo ineficaz cuando la superficie del agua suba lo suficiente, debido a
que el área ahora transportará agua en dirección aguas abajo para lograr una mayor elevación de la
superficie del agua. Por lo tanto, también se debe ingresar un parámetro adicional, llamado “elevación de
activación” para cada área de flujo ineficaz. Esta elevación del disparador también requiere el criterio del usuario.

Los modelos 2D no requieren que el usuario defina áreas de flujo ineficaces. Las áreas de flujo ineficaz, o
zonas de recirculación, se calculan automáticamente en función de las ecuaciones 2D. Sin embargo,
los modelos de turbulencia y los coeficientes de turbulencia pueden afectar el tamaño de las zonas de
recirculación y la velocidad del agua cerca de estas zonas.

6­11
Machine Translated by Google

Capítulo 6 TD­41

Ejemplos de aplicación

Como se mencionó anteriormente, definitivamente existen áreas de aplicación del modelado en las que el modelado 2D
puede producir mejores resultados que el modelado 1D. A continuación se muestran algunos ejemplos de aplicaciones en
las que el modelado 2D siempre dará mejores resultados que el modelado 1D:

• Diques que se desbordan y/o rompen donde el agua puede ir en muchas direcciones. Si el área interior de un
dique tiene una pendiente, el agua viajará por tierra en potencialmente muchas direcciones antes de encontrar
su camino hacia el punto más bajo del área protegida, y luego el agua comenzará a acumularse y
potencialmente sobrepasará y/o romperá el dique. en el extremo inferior del sistema. Sin embargo, si un área
protegida es pequeña y, en última instancia, toda el área se llenará hasta alcanzar el nivel de una piscina,
entonces el modelado 1D es adecuado para predecir la superficie final del agua y la extensión de la inundación.
En la Figura 6­9 se muestra un ejemplo de un sistema de diques con una brecha.
Como se muestra, el agua va en muchas direcciones, por calles, a través de caminos y hacia áreas bajas.

• Áreas y/o eventos en los que no se conoce completamente el recorrido del flujo del agua. Por ejemplo, si se va a
realizar un análisis de rotura de presa, pero el área aguas abajo no está confinada o es muy plana,
entonces puede que no sea posible tener un conocimiento exacto de hacia dónde irá el agua. Los modelos
1D requieren conocimiento de la ruta del flujo para todos los eventos antes de poder definir la geometría del
modelo. Los modelos 2D no tienen este requisito, ya que el flujo puede ir en cualquier dirección dentro del
dominio del modelado 2D. Además, debido a que una rotura de presa es altamente no lineal, se debe
utilizar la forma completa de las ecuaciones (incluidos los términos de aceleración), ya que una
rotura de presa tiene cambios extremadamente rápidos en profundidad y velocidad con respecto al espacio y
el tiempo. En la Figura 6­10 se muestra un ejemplo de rotura de presa, en el que el flujo sale a un área muy
plana y luego viaja en muchas direcciones.

• Llanuras de inundación muy amplias y planas, de modo que cuando los caudales salen al margen
En esta área, el agua tomará múltiples rutas de flujo y tendrá diferentes elevaciones y velocidades de la
superficie del agua en múltiples direcciones.

6­12
Machine Translated by Google

TD­41 Capítulo 6

Figura 6­9. Ejemplo de rotura de un sistema de diques con agua yendo en muchas direcciones.

Figura 6­10. Ejemplo Rompe presa que sale a una zona extremadamente plana y se extiende. Las profundidades del
agua se muestran en tonos de azul (el azul oscuro indica una mayor profundidad del agua).

6­13
Machine Translated by Google

Capítulo 6 TD­41

• Bahías y estuarios: dentro de las bahías y estuarios, el flujo irá continuamente en múltiples direcciones debido
a las fluctuaciones de las mareas y los flujos de los ríos que ingresan a la bahía/estuario en múltiples
lugares y momentos. En la Figura 6­11 se muestra el estuario del río Columbia.
Como se muestra, el agua se moverá en muchas direcciones dentro del estuario. Para obtener
velocidades y direcciones de flujo precisas, sería necesario un modelado 2D o incluso 3D.
Los ríos impulsados por las mareas también requieren el uso de la forma completa de las ecuaciones,
incluidos los términos de aceleración. La propagación de ondas en general no se puede resolver con
precisión sin los términos de aceleración en la ecuación del momento.

Figura 6­11. Estuario inferior del río Columbia con profundidades de agua mostradas en tonos de azul.

• Arroyos muy trenzados: en un arroyo muy trenzado, el canal se divide y

combina numerosas veces. Este tipo de sistema es muy difícil de modelar con precisión con un enfoque
de modelado 1D.

• Abanicos aluviales: es discutible que cualquier modelo numérico de límite rígido pueda capturar con precisión
un evento de inundación en un abanico aluvial, debido a la naturaleza episódica de las evoluciones del flujo
que pueden cambiar la superficie del terreno, el transporte y la dirección completa de los canales durante el
proceso. evento. En general, el agua se dividirá en muchas direcciones y se esparcirá sobre el abanico
aluvial. Esto sería extremadamente difícil de modelar en 1D.

• Flujo alrededor de curvas abruptas, en las que se produce una cantidad significativa de superelevación.
Por ejemplo, la Figura 6­12 muestra un modelo 2D de un canal rectangular de concreto que rodea una
curva de 180 grados. Se utilizó un modelo 2D para capturar la superelevación.

6­14
Machine Translated by Google

TD­41 Capítulo 6

que ocurre alrededor de la curva. Las velocidades que ocurren en este tipo de ubicación son muy
complejas en las tres dimensiones y muy probablemente tendrían patrones de flujo en espiral. Si se
necesitara una distribución de velocidad detallada a través de la curva, entonces un modelo 3D o un
modelo físico sería más adecuado para este tipo de problema. En otras palabras, esta es otra área
donde se requieren las ecuaciones 2D y 3D completas, ya que los términos de aceleración se utilizan
para modelar la elevación alrededor de la curva exterior. El avance de las olas en general sólo se
puede calcular con la inclusión de los términos de aceleración, y esto incluye el ascenso de la
superficie del agua alrededor de un objeto o una pared.

Figura 6­12. Ejemplo de superelevación de la superficie del agua en una curva cerrada.

• Aplicaciones donde existen gradientes de velocidad vertical. Estos problemas pueden requerir la
aplicación de un modelo 2D o un modelo 3D. En la Figura 6­13 se muestra un ejemplo de un modelo
2D promediado lateralmente que representa el gradiente de velocidad vertical.

Figura 6­13. Modelo detallado 2D promediado lateralmente de gradientes de velocidad vertical. Las direcciones
se muestran x y z.

• Aplicaciones en las que es muy importante obtener velocidades detalladas para la hidráulica del flujo sobre
o alrededor de un objeto, como el estribo de un puente o los pilares de un puente. Este tipo de
aplicación también requiere la forma completa de las ecuaciones, ya que la forma de onda de
difusión no puede explicar el avance de la onda frente al objeto, ni la rápida desaceleración y
aceleración frente al objeto y alrededor de él. En la Figura 6­14 se muestra un ejemplo de un
modelo 2D que representa velocidades detalladas del flujo que circula alrededor de los pilares. Sin embargo, el flujo

6­15
Machine Translated by Google

Capítulo 6 TD­41

alrededor de obstrucciones es implícitamente de naturaleza 3D; por lo tanto, si se requiere un

comportamiento detallado de la velocidad, la socavación o los sedimentos, entonces es necesario un modelo 3D no hidrostático.

Figura 6­14. Gráfico detallado de velocidad del modelo 2D del flujo que circula alrededor de los muelles desde la plataforma de una estación
de ferrocarril.

Las siguientes son áreas en las que el modelado 1D puede potencialmente producir resultados tan buenos como el modelado
2D (desde la perspectiva de las elevaciones calculadas de la superficie del agua y los hidrogramas de flujo/etapa),
con menos esfuerzo (desde el punto de vista del desarrollo, calibración y aplicación del modelo, como así como un
punto de vista del tiempo computacional):

Ríos y llanuras aluviales en los que las direcciones y fuerzas de flujo dominantes siguen la trayectoria general del
flujo del río, y donde el canal y las llanuras aluviales están directamente conectados y la superficie del agua
a lo largo de las secciones transversales es relativamente constante. Esta situación cubre potencialmente muchos
sistemas fluviales, pero es obviamente discutible en cuanto a la importancia de que las velocidades y fuerzas
laterales y verticales impacten las elevaciones calculadas de la superficie del agua y el límite de inundación
resultante. Un ejemplo de este tipo de río se muestra en la Figura 6­15, que muestra la confluencia de los ríos
Allegheny y Monongahela en Pittsburgh, PA. En esta situación de ejemplo, los dos ríos fluyen a través de colinas y
pequeñas cadenas montañosas. Los ríos están profundamente arraigados en la topografía y hay muy poca zona
de ribera para que el agua se esparza cuando el flujo sale del canal. Las líneas de flujo siguen muy de cerca
la línea central de la corriente y, por lo tanto, el agua es altamente unidimensional.

6­16
Machine Translated by Google

TD­41 Capítulo 6

Figura 6­15. Ejemplo de un sistema fluvial de flujo altamente unidimensional (Allegheny ­ Monongahela
Ríos, confluencia en Pittsburgh, PA). El color verde a rojo indica el terreno (elevación baja a alta)
y los azules indican la profundidad del agua (el azul oscuro indica mayor profundidad).

Arroyos empinados que son impulsados altamente por gravedad y tienen pequeñas áreas de desbordamiento.
En consecuencia, esta situación es también un área donde la forma de onda de difusión de las ecuaciones
producirá resultados similares a la forma completa de las ecuaciones. La razón por la que las ecuaciones producen
resultados similares se debe al hecho de que la magnitud de los términos de aceleración es muy pequeña en comparación
con las fuerzas de gravedad y de fricción que dominan el movimiento del agua en corrientes empinadas.

Sistemas fluviales que contienen muchos puentes/cruces de alcantarillas, presas, presas y otras estructuras cerradas, diques,
estaciones de bombeo, etc., que impactan las etapas y flujos calculados dentro del sistema fluvial. Para este tipo de sistemas
fluviales el estado actual del arte en modelos 1D todavía está por delante de los modelos 2D; sin embargo, esto
cambiará con el tiempo a medida que mejore el software de modelado 2D. Además, este también es un buen
lugar para realizar modelados combinados 1D y 2D para el sistema fluvial. Modelado 1D para el canal principal y
estructuras hidráulicas, y luego modelado 2D para las áreas de sobreterramiento y llanuras de inundación.

En la Figura 6­16 se muestra un ejemplo de un modelo combinado 1D/2D del sistema del río Truckee que fluye a través de
Reno, NV y luego pasa por el área de Truckee Meadows.
El río principal es empinado y tiene muchos puentes que cruzan el arroyo. Mientras que el flujo en el río principal es
altamente unidimensional, el flujo en las áreas sobre la ribera es altamente bidimensional. La modelización del río principal
de Truckee es complicada debido a los numerosos puentes y presas que cruzan el arroyo, diques intermitentes
que se

6­17
Machine Translated by Google

Capítulo 6 TD­41

rebasado en ambas direcciones y diferencias significativas en la superficie del agua entre el canal y los
márgenes. Posteriormente, esta es una situación de ejemplo ideal para utilizar el modelado 1D para
el río principal y el modelado 2D para todas las áreas de ribera.

Figura 6­16. Ejemplo de modelo 1D/2D del río Truckee cerca de Reno, NV. El color verde a rojo indica el
terreno (elevación baja a alta) y los azules indican la profundidad del agua (el azul oscuro indica
mayor profundidad).

Aplicaciones de pronóstico en tiempo real que modelan sistemas fluviales medianos a grandes (100 millas
o más del sistema fluvial) y ejecutan pronósticos de períodos de tiempo más largos (es decir, semanas
a meses). Incluso con los tremendos avances en la computación multiprocesador y en la computación
GPU (Unidades de Procesador de Gráficos), todavía existen importantes limitaciones espaciales y
de tiempo de simulación sobre qué y cuándo podemos usar efectivamente los modelos 2D en el
pronóstico en tiempo real. Los tiempos de ejecución de las simulaciones para el modelado 2D seguirán
mejorando con el tiempo.

Áreas en las que los datos básicos no respaldan el beneficio potencial de utilizar un modelo 2D.
Si no existe un terreno de ribera detallado y una batimetría de canal, o si solo hay secciones
transversales detalladas en los lugares estudiados, muchos de los beneficios del modelo 2D no se
obtendrán debido a la poca precisión de los datos del terreno.

6­18
Machine Translated by Google

TD­41 Capítulo 6

Modelo de desarrollo
Si bien el desarrollo de un modelo hidráulico es similar, independientemente de si es 1D o 2D, también existen
muchas diferencias significativas en el proceso de desarrollo del modelo. Esta sección analiza el proceso de
desarrollo del modelo y destaca las diferencias entre los dos enfoques. Además, se analizan las ventajas y
desventajas de los enfoques de modelado 1D y 2D.

Datos del terreno

Como se analizó anteriormente, todo desarrollo de modelos hidráulicos requiere datos del terreno. Específicamente,
los modelos 2D requieren que se construya un modelo de terreno completo que incluya la topografía de todo el dominio
de modelado, incluido el canal. Por otro lado, para los modelos 1D es común utilizar secciones transversales
topográficas para el canal y las áreas de llanura aluvial, o un modelo de terreno para las áreas de sobrebanco y
secciones transversales topográficas del área del canal. Para obtener más detalles sobre las necesidades del terreno
para los enfoques de modelado 1D y 2D, consulte el Capítulo 4.

Discretización del modelo

Antes de diseñar un modelo 1D, se debe conocer la trayectoria del flujo del agua para todos los eventos que se
van a modelar. Además, si las rutas/dirección del flujo cambian con diferentes eventos, esto requeriría más de una
representación geométrica del sistema. Por otro lado, los modelos 2D no requieren conocimiento de las trayectorias
exactas del flujo del agua. Requieren conocimiento del alcance de las inundaciones (lo que también es válido para los
modelos 1D).

Una vez definidas las extensiones del modelo, el canal y el área de la llanura aluvial se definen con secciones
transversales (modelos 1D) o una malla computacional (modelos 2D). Para el enfoque de modelado 1D, el nivel de
detalle se define por la cantidad de secciones transversales que se utilizan y por garantizar que las secciones
transversales se coloquen en las ubicaciones correctas. En general, para el modelado 1D, las secciones
transversales deben ubicarse espacialmente para capturar con precisión el terreno, las pérdidas por fricción y las
pérdidas por contracciones y expansión. Sin embargo, el número de secciones transversales necesarias para el
modelado 1D depende en gran medida de la pendiente del arroyo (es decir, los arroyos más empinados requieren más
secciones transversales a distancias más cortas entre sí).

Además, se deberán colocar secciones transversales en las roturas de rasante (cambios de pendiente);
alrededor de estructuras hidráulicas; a cambios significativos en la rugosidad; uniones y divisiones de flujo; y
lugares donde el flujo cambiará. En el modelado 1D, para obtener una representación precisa del área de flujo real, los
modeladores deben diseñar secciones transversales perpendiculares al flujo, así como definir áreas de flujo
ineficaces y obstrucciones bloqueadas al flujo dentro de cada sección transversal.
Esto requiere conocimiento de las rutas de flujo y qué parte de la sección transversal estará activa para la variedad de
eventos.

El desarrollo de un modelo 2D consiste en definir un tamaño de celda/elemento de malla base para el problema y luego
refinar/engrosar la malla en las áreas apropiadas. En general, se requieren tamaños de celda más pequeños para
ubicaciones donde el terreno cambia significativamente y donde la elevación de la superficie del agua y/o la velocidad
están cambiando. Este requisito se debe al hecho de que, dependiendo del modelo utilizado, se calcula una única
superficie de agua para cada celda y una velocidad promedio es

6­19
Machine Translated by Google

Capítulo 6 TD­41

calculado para cada cara de celda. Dependiendo del nivel deseado de detalle de las superficies y velocidades del agua, se
necesitarán más celdas en áreas donde se están produciendo más cambios.

El desarrollo de una buena malla 2D tendrá requisitos adicionales además de los enumerados anteriormente. A menudo es
necesario refinar y alinear las celdas contenidas con el canal principal para obtener mejores resultados. Esto incluirá
alinear las caras de la celda a lo largo del terreno elevado de las estaciones principales del banco del canal, para
garantizar que el flujo permanezca en el canal hasta que se exceda la elevación del banco. Alinear las celdas con el flujo
también reducirá cualquier difusión numérica que pueda estar presente en el esquema de solución que se utiliza para
resolver las ecuaciones. Y finalmente, el número de celdas a lo largo del canal afectará los detalles de la distribución de
velocidad calculada. Además, la mayoría de los modelos 2D representan computacionalmente cada celda como una única
elevación y cada cara como una sola línea. Algunos modelos utilizan triángulos para cada celda, que representan la
superficie terrestre con 3 puntos. Estos tipos de modelos requerirán más celdas a lo largo del canal para poder
representar correctamente el área del canal y el perímetro mojado.

Algunos modelos 2D utilizan tecnología de subcuadrícula (HEC­RAS), que representa la celda como una relación detallada
de elevación­volumen y cada cara como una sección transversal detallada, según el terreno subyacente. Los
modelos de subred requerirán muchas menos celdas para representar correctamente el área del canal y el perímetro mojado.
Este tipo de modelos pueden modelar un canal con precisión con solo unas pocas celdas para todo el canal (3 a 4
celdas). Sin embargo, si el usuario desea una distribución de velocidad detallada dentro del canal, entonces se necesitarán
más celdas, ya que cada celda/cara solo calcula una única velocidad promedio. En general, para obtener una buena
distribución de la velocidad en todo el canal, necesitará utilizar al menos de 5 a 10 celdas en todo el canal, dependiendo
de qué tan detallada desee que se discretice la velocidad.

Los requisitos adicionales para crear buenas mallas 2D incluirán alinear celdas a lo largo de la parte superior de las carreteras,
diques, muros contra inundaciones y cualquier otra barrera para el flujo, incluso si se trata de una barrera natural del terreno.

En el modelado 2D no es necesario definir áreas de flujo ineficaces, ya que el enfoque 2D puede identificar automáticamente
áreas de flujo activo y zonas de remolinos. Algunos modelos 2D no requieren que las celdas de la malla estén alineadas con
la dirección del flujo. Sin embargo, alinear las celdas de la malla para el flujo del canal principal permitirá que el modelo
produzca menos error numérico y una menor difusión de la elevación de la superficie del agua (atenuación del pico). En la
Figura 6­17 se muestra un ejemplo de una malla 2D con las celdas del canal principal más pequeñas que las celdas del
banco superior y las celdas/elementos del canal principal alineados con el flujo.

Si se requieren descripciones detalladas del comportamiento de la velocidad alrededor de estructuras como pilares de
puentes, vertederos en curva u otras estructuras, se debe proporcionar suficiente resolución alrededor de ellas para
permitir que el modelo 2D capture con precisión el comportamiento de interés. En la Figura 6­18 se muestra un ejemplo de
una malla detallada 2D de mayor resolución alrededor de una estructura.

6­20
Machine Translated by Google

TD­41 Capítulo 6

Figura 6­17. Ejemplo de malla 2D con malla detallada del canal principal, con rejillas alineadas al
flujo. El color verde a rojo indica el terreno (de baja a alta elevación).

6­21
Machine Translated by Google

Capítulo 6 TD­41

Figura 6­18. Ejemplo de malla 2D con una malla detallada de mayor resolución alrededor de las estructuras. Verde
El color rojo indica el terreno (de baja a alta elevación).

Definición de coeficientes de rugosidad

La definición espacial de los valores de fricción (coeficientes de rugosidad) es muy diferente entre los enfoques de
modelado 1D y 2D. Ambos enfoques pueden basarse en la estimación de valores iniciales a partir de datos de uso/
cobertura del suelo y fotografías aéreas. Sin embargo, este enfoque de estimación no es necesario para los modelos
1D. En cambio, para los modelos 1D, se ingresan coeficientes de rugosidad para cada sección transversal,
generalmente con un mínimo de valores de sobrebanco izquierdo, canal y sobrebanco derecho. Además, los
coeficientes de rugosidad se pueden discretizar aún más horizontalmente tanto para las áreas sobre el banco como
para el canal principal. También un enfoque común para los modelos 1D es tener opciones para variar la
rugosidad verticalmente en las secciones transversales, o poder variar la rugosidad con cambios en el flujo.
tarifas.

Además, algunos modelos 1D incluso tienen opciones para cambiar la rugosidad debido a cambios estacionales
de vegetación o temperatura. Si bien obtener estimaciones detalladas de la rugosidad inicial en un modelo 1D puede
llevar algo de tiempo, cambiar la rugosidad para la calibración del modelo es muy fácil. Con el uso de tablas para
mostrar y modificar valores de rugosidad, se pueden cambiar rápidamente porciones más grandes de un modelo y se
puede calibrar un modelo de manera iterativa. Además, los factores de flujo versus rugosidad y rugosidad estacional
también facilitan la calibración de un modelo, así como la prueba de la sensibilidad del modelo a los valores de rugosidad
calibrados.

Desarrollar estimaciones iniciales de coeficientes de rugosidad para modelos 2D es similar a 1D, en el sentido de que
el modelador puede utilizar información de fotografía aérea y uso del suelo para crear capas espaciales que pueden
relacionarse con estimaciones iniciales o coeficientes de rugosidad. Además, la mayoría de los modelos 2D tienen
una forma de anular los coeficientes de rugosidad iniciales con polígonos definidos por el usuario. Estos polígonos
se pueden utilizar para establecer mejores estimaciones iniciales de los coeficientes de rugosidad, por ejemplo

6­22
Machine Translated by Google

TD­41 Capítulo 6

definir la rugosidad para las áreas principales del canal, o los polígonos se pueden usar para cambiar las
estimaciones iniciales durante el proceso de calibración.

Cambiar los coeficientes de rugosidad para el proceso de calibración generalmente requiere mucho más tiempo
y es más difícil para los modelos 2D que para los modelos 1D. La razón de esta diferencia es que los modelos 2D
requieren que se defina la rugosidad para cada cara de la malla 2D. Además, el proceso de calibración es más difícil
para los modelos 2D porque se deben tomar decisiones sobre qué extensión del área espacial requiere que se cambie
la rugosidad para poder calibrar un modelo.

Estructuras Hidráulicas

La definición y el modelado de estructuras hidráulicas pueden ser muy diferentes entre los enfoques de modelado 1D y
2D. En general, el modelado 1D suele basarse en ecuaciones y coeficientes semiempíricos para modelar
estructuras hidráulicas. La hidráulica de la estructura se puede definir completamente con una única o familia de
curvas de potencia. Sin embargo, también hay aplicaciones en las que se utiliza la ecuación de energía o de
momento (usando secciones transversales) para resolver la superficie del agua y la velocidad hacia arriba, a través
y fuera de la estructura. Para este tipo de modelado, se requiere la ecuación de momento completa (no la forma de
onda de difusión) para capturar con precisión la aceleración y desaceleración a través de la estructura.

Para el modelado 2D de estructuras hidráulicas, se puede adoptar un enfoque similar al modelado 1D, y lo más habitual
es que se aplique a la gran mayoría del modelado de estructuras. Sin embargo, existe la opción de modelar el flujo a
través de la estructura en un enfoque completamente bidimensional, que requiere mucho más trabajo y tiempo pero
produce información mucho más detallada. Un ejemplo de aplicación de este enfoque es el flujo a través de
un puente. Si el agua no va a llegar a la plataforma del puente y va a estar presurizada, entonces se puede
utilizar un enfoque de modelado 2D detallado (usando muchas celdas/elementos pequeños). En este caso, el propósito
del enfoque de modelado 2D es capturar más detalles de cómo la superficie del agua y las velocidades cambian a
medida que el flujo se acerca al puente, atraviesa el puente y se expande a medida que el agua sale del puente, así
como capturando separaciones de flujo y zonas de remolinos. Un ejemplo del enfoque de modelado 2D detallado para
un puente se muestra en la Figura 6­19 (la figura muestra las velocidades resultantes a través de un puente a partir de
un modelo 2D detallado del puente).

Este mismo enfoque se puede utilizar para muchos tipos de estructuras hidráulicas, según el nivel de detalle deseado.
Sin embargo, para estructuras donde el flujo pasa a través de una profundidad crítica y entra en caída libre, actuando
como una cascada, este enfoque a menudo producirá resultados incorrectos. Las ecuaciones de aguas poco
profundas, tal como se definen y utilizan en casi todos los programas 2D, se definen para un flujo que varía
gradualmente y bajo el supuesto de una distribución de presión hidrostática consistente a través del campo
de flujo. Estos dos supuestos no son ciertos para el flujo que pasa por una profundidad crítica y entra en caída
libre.

Incluso cuando el flujo general de un río es constante (flujo de entrada = flujo de salida), el campo de flujo en las
proximidades de una estructura hidráulica puede ser inestable. Este fenómeno se puede capturar con modelos 2D y
3D, pero no con modelos 1D.

6­23
Machine Translated by Google

Capítulo 6 TD­41

Figura 6­19. Ejemplo de salida de velocidad para un modelo 2D detallado de un puente (la velocidad se superpone al terreno)
donde el color verde a rojo indica una elevación baja a alta).

Condiciones de borde

La definición de condiciones de contorno es similar entre los enfoques de modelado 1D y 2D. En los enfoques
de modelado 1D, todas las condiciones de contorno se adjuntan directamente a una sección transversal, un
área de almacenamiento o una estructura hidráulica interna. En los enfoques de modelado 2D, las
condiciones de contorno se adjuntan a celdas/elementos de límites externos, celdas/elementos/caras internas y
también a estructuras hidráulicas internas. Se pueden utilizar todos los mismos tipos de condiciones de contorno:
flujo, etapa, profundidad normal (ecuación de Manning); curvas de calificación; precipitación; aguas subterráneas;
y formas de controlar estructuras hidráulicas durante la simulación. La única diferencia es que en el
modelado 2D, algunos de los tipos de condiciones de contorno permiten variaciones en las elevaciones de la
superficie del agua a lo largo de la ubicación de la condición de contorno (es decir, flujo de salida en profundidad
normal). Además, algunos tipos de condiciones de contorno de modelos 2D también pueden requerir información
de velocidad para capturar el impulso del flujo con mayor precisión a medida que ingresa al sistema.

Una consideración adicional para elegir la ubicación de la condición de contorno en un modelo 2D es la necesidad
de alejar esa condición de contorno lo suficiente del área de interés para estar seguro de que la ubicación de
la condición de contorno no influye en los resultados. Un ejemplo podría ser un arroyo con un gran desbordamiento
en un lado donde se aplica la condición de límite aguas arriba. Aplicar el flujo a esa condición de contorno que se
está tratando como un elemento 1D podría sesgar la distribución del flujo en el banco superior a cierta
distancia aguas abajo del

6­24
Machine Translated by Google

TD­41 Capítulo 6

condición de contorno ya que la distribución de flujo 1D supuesta se propaga aguas abajo dentro del dominio
2D.

Calibración del modelo

El proceso de calibración del modelo es uno de los pasos más importantes en el desarrollo de un modelo
hidráulico. Se requiere la calibración de cualquier modelo hidráulico para comprender si el modelo es capaz de
reproducir inundaciones pasadas y si puede predecir inundaciones futuras con confianza. El proceso de
calibración también permite una mayor comprensión de la sensibilidad del modelo a los datos, las fuerzas
de fricción y otros coeficientes empíricos. Un modelo que no está calibrado es sólo un experimento numérico.
Además, un modelo no calibrado puede o no estar siquiera cerca de reproducir elevaciones y flujos realistas de
la superficie del agua en todo el sistema.
Sin embargo, a veces, en situaciones de emergencia, donde no hay datos para calibrar el modelo o no hay
tiempo suficiente para calibrarlo, los resultados de un modelo no calibrado son mejores que ningún resultado.

El proceso de calibración de un modelo es muy similar entre los enfoques de modelado 1D y 2D; sin
embargo, la calibración del modelo 2D puede ser más difícil y llevar más tiempo. Los datos necesarios para
realizar una calibración de modelo son prácticamente los mismos, excepto que una calibración más
detallada de modelos 2D requiere mediciones detalladas de la distribución de velocidad para garantizar la
precisión de los perfiles de velocidad que pueden estar disponibles en ubicaciones calibradas donde se
realizan mediciones directas para un rango de flujos por agencias como el Servicio Geológico de Estados
Unidos. La cantidad y calidad de los datos disponibles variarán de un lugar a otro. En algunos casos, es posible
que no haya datos observados suficientes para realizar una calibración completa, sino que se puede caracterizar
mejor como una calibración limitada de ciertos aspectos de la simulación o en partes limitadas del dominio del
modelo.

Los datos generales para la calibración del modelo son: cobertura/uso del suelo observado (posiblemente para
diferentes años o épocas del año); registros de lluvia (o datos de lluvia cuadriculados); elevaciones de la
superficie del agua observadas en medidores de corriente; flujos calculados en medidores de corriente (basados
en la superficie del agua observada y una curva de clasificación establecida); marcas altas de agua por
escombros y manchas de agua; extensiones/límites de las inundaciones a partir de fotografías aéreas e
inspecciones de campo; relatos anecdóticos/elevaciones de agua de investigaciones de campo de hogares y
negocios que fueron impactados durante inundaciones históricas; y cualquier fotografía o video que se tomó durante la inundación o poc
Toda esta información se utiliza en el proceso de calibración del modelo 1D o 2D. Además de estos datos, cualquier
medición detallada de la velocidad que se tome en las secciones transversales a menudo puede ser muy útil
para una mayor calibración de un modelo 2D. Este tipo de información de velocidad generalmente se recopila en
estructuras de interés o lugares de interés, como: puentes, alcantarillas, presas, aliviaderos, curvas
cerradas, alrededor de diques y también cerca de estructuras de formación de ríos. Este proceso a menudo
requiere una mirada crítica a todos los datos en los que se basa. Los medidores de caudal registran la etapa,
que se convierte en un caudal mediante una curva de clasificación que puede tener sesgo. Las mediciones de
nivel de agua pueden estar sujetas a problemas mecánicos o tener influencias hidráulicas localizadas en la lectura
creadas por la estructura a la que están adheridas, etc. Las marcas de agua altas pueden ser sospechosas si se
recogen en áreas donde influyen la acción del oleaje, el peralte u otros factores. los niveles de agua en ese lugar
y las elevaciones obtenidas con una unidad GPS portátil que tiene una precisión vertical y horizontal limitada.

6­25
Machine Translated by Google

Capítulo 6 TD­41

El proceso de calibración de un modelo hidráulico comienza con la hidrología. Un modelo hidráulico es tan bueno
como las condiciones límite de flujo proporcionadas. Tener buenas condiciones de contorno es imperativo en
cualquier tipo de modelado de superficies terrestres y riberas. Generalmente, se utilizará un modelo hidrológico
para realizar el modelado de toda la cuenca y proporcionar datos de flujo al modelo hidráulico para las
condiciones de contorno. Entonces, el primer paso es garantizar que el modelo hidrológico esté bien calibrado.
Durante la calibración del modelo hidráulico, puede ser necesario volver al modelo hidrológico para calibraciones/
cambios adicionales basados en los resultados del enrutamiento del modelo hidráulico y los cálculos de
elevaciones de la superficie del agua, flujos y tiempos del hidrograma. Esto debe considerarse como un paso
común, en el que el modelado hidrológico y las calibraciones del modelado hidráulico se realizan juntos, no como
ejercicios separados.

La calibración de los alcances del modelado 1D se logra cambiando/ajustando lo siguiente:

• Parámetros de rugosidad hidráulica. •

Coeficientes de contracción y expansión. • Extensiones
de área de flujo ineficaces y elevaciones desencadenantes de altura. •
Coeficientes de estructura hidráulica. •
Coeficientes de pérdida por flexión (a veces llamados pérdidas menores).
• Información sobre las condiciones de los límites, como pendientes de energía o incluso curvas de clasificación potenciales.
valores.
• Información de bloqueo de escombros en estructuras.
• Dimensiones de ruptura del dique y valores de tiempo.

6­26
Machine Translated by Google

TD­41 Capítulo 6

Los ajustes del modelo de los coeficientes de rugosidad se realizan según el alcance, en lugar de secciones
transversales individuales. Los cambios deben realizarse gradualmente entre ejecuciones del modelo, para comprender
mejor los efectos de esos cambios tanto en la elevación de la superficie del agua como en la atenuación del
hidrograma y el tiempo de viaje con el objetivo final de hacer que el modelo imite más fielmente el entorno
físico. y características del flujo. Los cambios realizados aguas arriba en un modelo afectarán localmente las elevaciones
de la superficie del agua, pero también afectarán la magnitud y el momento de los flujos que van aguas abajo. Los
cambios realizados aguas abajo en un modelo no sólo afectarán las elevaciones de la superficie del agua localmente,
sino que también pueden afectar las elevaciones de la superficie del agua aguas arriba a cierta distancia. Cuanto
más plano sea el arroyo, mayor será la distancia que se verá afectada por los cambios en la elevación de la
superficie del agua aguas abajo. En la Figura 6­20 se muestra un ejemplo de un tramo de río modelo 1D calibrado para
el Bajo Río Columbia. En la Figura 6­20 se muestran tanto los datos observados de los medidores como las marcas
de agua altas recopiladas después de la inundación. Hay una dispersión notable en algunos de los datos
de la marca de límite superior, lo cual es común, ya que este tipo de datos observados pueden tener una
dispersión e incertidumbre considerables en los valores.

Figura 6­20. Ejemplo de un modelo 1D calibrado para el Bajo Río Columbia.

Calibrar un modelo 2D es muy similar, pero con algunas de las siguientes diferencias. En el modelado 1D, los
cambios en la rugosidad se realizan en secciones transversales y/o grupos de secciones transversales dentro
de un tramo. En general, ajustar los coeficientes de rugosidad, ya sea hacia arriba o hacia abajo según sea necesario,
para una variedad de secciones transversales es fácil mediante el uso de tablas de valores y factores de rugosidad
para aumentar o disminuir los valores en un porcentaje. Además, los factores de rugosidad que varían según el flujo y/
o la estación generalmente están disponibles y son fáciles de usar.

Ajustar la rugosidad en un modelo 2D es más difícil, porque el modelador no sólo tiene que elegir una distancia a lo largo
del río para hacer ajustes, sino también decidir qué partes de esa área deben ajustarse. Esta mayor complejidad se
debe a la variabilidad espacial mucho mayor que existe

6­27
Machine Translated by Google

Capítulo 6 TD­41

definido para coeficientes de rugosidad de modelado 2D (es decir, cada cara de celda/elemento tiene
potencialmente un valor de rugosidad único). Además de la extensión espacial para realizar ajustes a los
coeficientes de rugosidad, otra decisión necesaria es cuánto ajuste se debe hacer a cada uno de los tipos de
uso del suelo de rugosidad base.

Si bien los modelos 2D no utilizan coeficientes empíricos de contracción y expansión, y no hay necesidad de
áreas de flujo ineficaces y obstrucciones bloqueadas en las secciones transversales, existen coeficientes
empíricos de turbulencia que pueden necesitar ser ajustados para tener en cuenta mejor el alcance de la
turbulencia. zonas de remolinos, así como la difusión natural del campo de flujo a medida que el hidrograma
viaja aguas abajo. Es discutible qué parte de la difusión de un hidrograma se debe a fuerzas de fricción
versus turbulencia. La calibración de modelos 2D también requerirá los mismos tipos de cambios que los
modelos 1D para los siguientes aspectos: cambios en los coeficientes de la estructura hidráulica que están
afectando los resultados; cambios en las condiciones de contorno; bloqueos de escombros e información
sobre rupturas de diques.

Como se mencionó anteriormente, se puede lograr un mayor refinamiento de los modelos 2D con mediciones de
velocidad observadas. Si se está desarrollando un modelo 2D para obtener un conocimiento detallado de la
distribución de la velocidad en áreas específicas, entonces haber observado la información de la velocidad,
incluso si la información es para eventos de menor magnitud, puede ser muy útil para garantizar que el
modelo produzca resultados razonables. La calibración de la información de velocidad generalmente requerirá
cambios en la rugosidad, posibles coeficientes de viscosidad de remolinos y tal vez incluso ajustes del
terreno, si es apropiado y justificado.

Problemas de tiempo y costos

El tiempo para desarrollar un modelo con el enfoque de modelado 1D o 2D puede variar, dependiendo del tipo y
propósito del modelo. Por ejemplo, si es necesario desarrollar un modelo “rápido y sucio” durante una
emergencia por inundación, entonces es mucho más rápido diseñar un polígono de área de flujo 2D, establecer
un tamaño básico de celda/elemento, adjuntar algunas condiciones de contorno y comenzar a crear simulaciones.
Hacer un modelo rápido con el enfoque 1D definitivamente llevará más tiempo. Sin embargo, si se está
desarrollando un modelo detallado, ambos enfoques finalmente tomarán aproximadamente la misma cantidad
de tiempo para desarrollar el modelo inicial. Sí, se requiere una cantidad considerable de trabajo para
diseñar secciones transversales para un modelo hidráulico detallado. Sin embargo, también se requiere
mucho trabajo para refinar un área de flujo 2D con una resolución de malla adicional. Además, algunos modelos
pueden requerir la delineación de líneas de ruptura para definir áreas de flujo 2D. Estimar los parámetros
base, como: coeficientes de rugosidad, coeficientes de estructura hidráulica y otros coeficientes
empíricos necesarios para los enfoques de modelado, lleva aproximadamente la misma cantidad de tiempo.

Además, otro factor que varía el tiempo de desarrollo del modelo es el proceso de calibración del modelo.
Durante este proceso se realizan ajustes al modelo y luego se calcula el modelo para todos los eventos
y escenarios deseados. La principal discrepancia de tiempo se debe al hecho de que los modelos 2D requieren
mucho más tiempo de cálculo que los modelos 1D. Por ejemplo, es común que un modelo 2D requiera varias
horas de cálculo, posiblemente incluso un día o más. Por otro lado, los modelos 1D generalmente se ejecutan
en minutos u horas. Debido a esta importante diferencia de tiempo computacional, el proceso de calibrar
y realizar análisis alternativos para una variedad de eventos consumirá mucho más tiempo utilizando el
enfoque de modelado 2D.

6­28
Machine Translated by Google

TD­41 Capítulo 6

Además de los requisitos de tiempo de cálculo, los requisitos de hardware para ejecutar modelos 1D y 2D son
generalmente diferentes. En general, los modelos 1D son de un solo subproceso, lo que significa que solo requieren un
procesador de un solo núcleo. Por otro lado, los modelos 2D generalmente pueden aprovechar tantos
procesadores centrales como estén disponibles en la máquina. Para un modelado 2D serio, el modelador debe
obtener una máquina computacional de alto nivel para dicho trabajo, o incluso investigar las opciones de HPC
(Computación de alto rendimiento) y computación en la nube.

Conocimientos, habilidades y habilidades del modelador

Ya sea que realice modelado 1D o 2D, el modelador debe tener una buena experiencia en hidrología, hidráulica,
hidráulica de estructuras y soluciones numéricas de ecuaciones no lineales. Además, no importa qué modelo
de computadora y enfoque de modelado se seleccione, el modelador debe tomar clases sobre cómo usar esa pieza
de software específica de manera efectiva para el modelado 1D y/o 2D. El autoestudio es fantástico para empezar a
utilizar un software específico, pero el conocimiento adicional adquirido al tomar una clase de los desarrolladores/expertos
en el uso del software a menudo resultará invaluable.

Adquirir conocimientos, habilidades y capacidades en el modelado con enfoques de modelado 1D o 2D es una cuestión
de tiempo y oportunidad. El modelador necesita realizar estudios detallados de los sistemas fluviales en los que
tenga buen terreno, información estructural y datos históricos para realizar calibraciones del modelo. Si el modelador
es nuevo en el modelaje, siempre debe buscar orientación y ayuda de modeladores experimentados. Además, todos
los modelos deben ser revisados por expertos independientes para garantizar que el trabajo realizado se haya
realizado con una elección de modelado adecuada (1D o 2D); aplicado correctamente; contiene la cantidad adecuada
de detalles para la ubicación y el tipo de estudio determinados; fue calibrado efectivamente y está reproduciendo
resultados hidráulicos razonables para los eventos y alternativas que se están modelando.

Resumen de ventajas y desventajas del modelado 1D y 2D

Existen claras ventajas y desventajas tanto para el enfoque de modelado 1D como para el 2D. Esta sección
del documento resume las ventajas y desventajas que se han descrito anteriormente en este documento. La discusión
se basa en las ventajas y desventajas desde la perspectiva del modelado 1D versus el modelado 2D.

Ventajas del modelado 1D:

En general, los modelos 1D requieren menos datos del terreno, ya que la parte del canal del modelo puede
provenir de levantamientos transversales detallados separados.

A menudo (pero no siempre) los modelos 1D son más fáciles de calibrar, debido a la simplicidad de
cambiar parámetros como los coeficientes de rugosidad y otras variables.

El modelado de estructuras hidráulicas suele ser más fácil y requiere menos datos y recursos computacionales.
requisitos.

Se requieren significativamente menos tiempo y recursos computacionales para los modelos 1D.

6­29
Machine Translated by Google

Capítulo 6 TD­41

El enfoque de modelado 1D puede ser más apropiado para áreas en las que los datos básicos no respaldan
la capacidad de desarrollar un modelo 2D razonable.

Desventajas del modelo 1D:

Se debe conocer la trayectoria del flujo del agua, para todos los eventos, antes de desarrollar el modelo.
No siempre es posible conocer la ruta del flujo, especialmente en áreas planas, lo que hace que esas
decisiones a menudo sean subjetivas y menos precisas.

Para obtener una representación precisa del área de flujo real, las secciones transversales 1D deben
disponerse perpendiculares al flujo. Sin embargo, cumplir este requisito no siempre es posible para toda la
gama de eventos y, por lo tanto, puede requerir más de una representación geométrica del sistema para
eventos de flujo bajo y alto que tienen diferentes patrones de flujo.

Las pérdidas de energía y/o fuerza debido a contracciones y expansiones requieren que el modelador defina
coeficientes empíricos y áreas de flujo ineficaces.

La dirección del flujo durante el evento se limita a la ruta de flujo definida del 1D
modelo.

La salida de velocidad se limita a los valores promedio para el canal principal y los overbanks. Si bien se puede
calcular una mayor discretización de la velocidad a partir de la salida 1D, se limita a la suposición de que el
flujo es perpendicular a la sección transversal y que la distribución del flujo es solo una función de la forma de
la sección transversal y los valores de rugosidad.

El mapeo del área inundada se basa en la suposición de que la superficie del agua
cambia linealmente entre dos secciones transversales cualesquiera, y que la superficie del agua sea plana dentro de las
áreas de almacenamiento.

Ventajas del modelado 2D:

No es necesario conocer la trayectoria del flujo del agua, para todos los eventos, para desarrollar el modelo.
Sin embargo, es necesario definir correctamente el alcance de las inundaciones.

La dirección del flujo puede cambiar durante el evento. El agua puede moverse en cualquier
dirección, basada en la energía y el impulso del flujo.

La velocidad, el impulso y la dirección del flujo se tienen en cuenta con mayor precisión con el modelado 2D.
Esta responsabilidad es especialmente cierta para el flujo que pasa por carreteras, diques, barreras,
estructuras, alrededor de curvas y en cruces/divisiones de flujo. Además, se pueden utilizar modelos 2D para
analizar zonas de remolinos dentro del campo de flujo. Alrededor de las curvas, los modelos 2D producen
elevaciones precisas de la superficie del agua, pero las distribuciones de velocidad pueden ser erróneas
debido a la existencia de flujo helicoidal.

Las pérdidas de energía y fuerza debidas a contracciones y expansiones, etc., se contabilizan directamente y
no requieren coeficientes empíricos, mayor rugosidad ni áreas de flujo ineficaces definidas por el usuario.

El mapeo del área inundada, así como las velocidades y los riesgos de inundación (profundidad x velocidad)
es más preciso.

6­30
Machine Translated by Google

TD­41 Capítulo 6

El modelado detallado de estructuras hidráulicas, en un enfoque de modelado 2D completo, puede

proporcionar más información sobre la distribución del flujo que se acerca, atraviesa y sale de una
estructura.

Desventajas del modelo 2D:

Se requieren modelos de terreno más precisos y detallados para desarrollar un modelo 2D preciso. El
terreno debe incluir los detalles de los canales en todas las ubicaciones dentro del modelo, así como
capturar correctamente características como carreteras, bermas y diques.
Es posible que los conjuntos de datos LiDAR demasiado filtrados o los conjuntos de datos que se hayan procesado con un tamaño
de cuadrícula demasiado grande no resuelvan adecuadamente estas características clave del terreno que influyen en los
comportamientos y patrones del flujo.

Definir y modificar valores de rugosidad requiere una mayor definición espacial y puede ser
más difícil y requiere más tiempo durante el proceso de calibración.

Requiere significativamente más tiempo computacional y/o recursos computacionales. Puede requerir
la compra de una computadora de muy alto nivel (muchos núcleos, CPU rápidas, mucha RAM y
disco duro rápido) o utilizar HPC y soluciones de computación en la nube. Puede

requerir el uso de tamaños de cuadrícula más grandes de lo deseable para el problema, a fin de reducir los
tiempos de ejecución a una cantidad de tiempo manejable.

Es posible que realmente no se produzcan mejores resultados si los datos utilizados para realizar el
modelado (terreno, datos del canal y rugosidad) no soportan el nivel requerido para un
modelado 2D preciso.

La decisión entre las opciones de modelos 1D o 2D debe mantener la fidelidad a los objetivos del proyecto que
se investiga y el principio del "camino más simple y técnicamente sólido para lograr los objetivos". La falta de
datos para un modelo 2D o la disponibilidad de un modelo 1D no justifica, en sí misma, la elección de un modelo
a menos que esté respaldado por los objetivos del proyecto. Por ejemplo, si la extensión de la inundación es
de interés para un río que incluye curvas u otras características 3D, un modelo 1D o 2D podría ser suficiente
aunque la distribución de velocidad sea incorrecta; por otro lado, si el comportamiento de los sedimentos es
de interés, podría ser más apropiado un modelo 2D modificado para flujo helicoidal o un modelo 3D. El ingeniero
debe aprovechar el conocimiento adecuado del sistema para alcanzar con éxito los objetivos del proyecto.

6­31
Machine Translated by Google
Machine Translated by Google

TD­41 Capítulo 7

Capítulo 7

Bidimensional versus tridimensional

Modelado
La cuestión del modelado hidráulico 2D versus 3D es una cuestión muy compleja. La razón principal de
esta complejidad es que la mayoría de los modelos en los que el 3D es apropiado incluyen amplias áreas
donde domina el comportamiento 2D. Desafortunadamente, las áreas 3D impactan las áreas 2D y viceversa.
Este capítulo aborda principalmente únicamente el modelado hidrostático 3D; sin embargo, tres ejemplos de
aplicaciones utilizan modelos 3D no hidrostáticos. De ahora en adelante, 3D implicará 3D­hidrostático, a menos
que se indique lo contrario.

Ecuaciones bidimensionales y tridimensionales

Las leyes físicas que gobiernan el flujo del agua son: (1) el principio de conservación de la masa
(continuidad) y (2) el principio de conservación del momento. Estas leyes se expresan
matemáticamente en forma de ecuaciones diferenciales parciales, a las que en adelante nos
referiremos como ecuaciones de continuidad y de momento. En la derivación de las
ecuaciones de flujo inestable bidimensionales y tridimensionales utilizadas en este capítulo,
se hacen varios supuestos sobre el flujo:

1. Se considera que el flujo es de un fluido incompresible.

2. La distribución de presión se considera hidrostática.
3. La aceleración vertical del agua se considera insignificante.
4. La pendiente del lecho se considera suave.
5. Los efectos de la fricción límite se pueden tener en cuenta con leyes de resistencia al flujo
derivadas para flujos estacionarios (p. ej., ecuación de Manning).
6. La aproximación de Bousinnesq es válida.

Las ecuaciones 2D se presentaron anteriormente en el Capítulo 6 y, por lo tanto, no se repetirán aquí.

Las ecuaciones de flujo 3D se presentan a continuación.

Ecuaciones tridimensionales (3D)

La ecuación tridimensional de continuidad y momento se puede escribir en forma de ecuación diferencial

parcial, con respecto a la presión (P), la velocidad (u, v, w), y comúnmente se muestran en los libros de
texto de hidráulica de la siguiente manera:

Ecuación de continuidad:

uvw
 0
xyz

7­1
Machine Translated by Google

Capítulo 7 TD­41

Ecuaciones de momento:

uuu
t ( x
) (
y
ultravioleta
) uuu
()
z
 
fv
1 P1
X
Tx 
 
 . S
X

0 0

v ( )vu  vv
() VW )

z
 fu

1 P1
 
 
.t y
Sy
t X y ( 0 y 0

, P
P xyz

, xygxyz
a , 
dz
  , ,

z

dónde:
tu = velocidad en la dirección x =
v velocidad en la dirección y w =
velocidad en la dirección z
Tx = Esfuerzos turbulentos en la dirección x

Ty = Esfuerzos turbulentos en la dirección y f =

Parámetro de Coriolis
PAG = Presión


 = Densidad de referencia
= Densidad que varía espacialmente
Pa = Presión atmosférica que varía espacialmente =
g gravedad
S X= Fuentes/sumideros de impulso adicionales en la dirección x

S y= Fuentes/sumideros de impulso adicionales en la dirección y

La densidad que varía espacialmente, , se calcula utilizando ecuaciones de estado para temperatura (T),
salinidad (S) y sedimento (Se), ya sea individualmente o como combinación. La ecuación de estado se representa como:

  
ST , 

 SS 

TT 
0 S 0 t 0

 
en 
Ci SG i
 1
i 1 SGi C
 SG
i i
 1
dónde:
s = 0,78 kg m­3/ppt =
0,15 kg m­3/ oC =
diferencia de densidad debida al sedimento
i = clase de sedimento, 1 a n
SG = Gravedad específica =
Concentración del sedimento clase C

La Figura 7­1 proporciona un diagrama que define los símbolos utilizados en las ecuaciones 3D.

7­2
Machine Translated by Google

TD­41 Capítulo 7

Figura 7­1. Definición de símbolos utilizados en ecuaciones 3D.

Ejemplos de aplicación
Hay un conjunto distinto de problemas/aplicaciones de ingeniería que gravitan hacia la aplicación de un modelo
numérico 3D, generalmente cuando el modelador está interesado en una descripción tridimensional de la
velocidad (batimetría compleja, curvaturas o estructuras hidráulicas cercanas), gradiente de densidad flujos
impulsados o flujos impulsados por gradiente de temperatura. A continuación se muestran algunos ejemplos
específicos de estos problemas/aplicaciones.

• Al modelar la estratificación vertical debido a la densidad variable (agua dulce/agua salada): Este es un
problema común que se encuentra en ambientes estuarinos donde el estuario es profundo o contiene
un canal de navegación profundo a través de áreas poco profundas (por ejemplo, Bahía Mobile, Bahía
Galveston), Figura 7­2. En tales aplicaciones, el agua salina empujada por las mareas hacia el
estuario provoca la estratificación del estuario (Figura 7­3), con agua dulce de fuentes
terrestres en la superficie y agua salina debajo. Esta estratificación también da como resultado que la
dirección del flujo entre la superficie y el fondo sea opuesta en ciertas condiciones de marea. La
Figura 7­3 también ilustra la variación de velocidad entre la superficie y el fondo. Esta variación en
la velocidad es una condición que los modelos 2D no pueden simular; los modelos 2D promediados
lateralmente pueden replicar esta estratificación pero no la distribución del flujo entre el canal
profundo y el poco profundo.

7­3
Machine Translated by Google

Capítulo 7 TD­41

Figura 7­2. Canal profundo rodeado de áreas poco profundas, el rojo indica más profundo y el verde menos profundo.

Figura 7­3. Estratificación de la salinidad y diferencias de velocidad en un estuario.

• Al modelar la estratificación debido a la temperatura: de manera similar a la salinidad,

La temperatura provoca diferencias de densidad entre la superficie y el fondo. Las diferencias de
densidad impulsadas por la temperatura se observan con mayor frecuencia en embalses o lagos.
El agua es particularmente receptiva a este fenómeno, donde durante el caluroso verano la parte
superior de la columna de agua es más cálida que la parte inferior, y durante el invierno la superficie es
más fría que el fondo.

7­4
Machine Translated by Google

TD­41 Capítulo 7

El segundo escenario en particular es problemático, ya que causa una estratificación inestable y

eventualmente resulta en un “vuelco” del agua donde el agua más fría se hunde y el agua más caliente
sube. El “vuelco” tiene consecuencias graves y negativas para la calidad del agua y la ecología. Además,
los modelos numéricos promediados lateralmente en 2D pueden simular esta estratificación de temperatura
de los yacimientos, pero son intrínsecamente incapaces de capturar la hidrodinámica asociada con el
"vuelco". La Figura 7­4 proporciona un ejemplo de estratificación de temperatura.

Figura 7­4. Estratificación de temperatura en un lago.

• Al modelar el flujo alrededor de una curva: el movimiento de un fluido incompresible, como el agua, alrededor
de una curva se rige por un equilibrio entre las fuerzas de inercia y de presión.
Las fuerzas de inercia hacen que el agua fluya hacia el exterior de la curva y da como resultado la creación
de un gradiente de elevación. El exterior de la curva alcanza una elevación mayor que el interior de la curva.
Esta "superelevación" provoca una diferencia de presión entre el exterior y el interior, lo que da como
resultado que el flujo en la parte inferior sea hacia el interior y el flujo en la parte superior sea hacia el exterior.
Esta condición de flujo se llama "helicoidal". La generación de “superelevación” se puede simular mediante
modelos 2D, y el comportamiento amplio de la distribución del momento debido al flujo helicoidal se puede
simular agregando términos de corrección de la distribución del momento a los modelos promediados 2D.

La Figura 7­5 presenta un campo de velocidad promedio en profundidad 2D corregido utilizando el

método de transporte de vorticidad (Bernard y Schneider, 1992; Finnie et al., 1999). La magnitud de las
correcciones aplicadas al campo de velocidad se calcula únicamente con el análisis 2D de aguas poco profundas.

7­5
Machine Translated by Google

Capítulo 7 TD­41

Las ecuaciones se presentan en la Figura 7­6. Las variaciones verticales en la velocidad pueden estimarse
a partir de estas correcciones para su uso en modelos de transporte cuasi­3D (Brown. 2012). Sin embargo,
la determinación de la distribución exacta de la velocidad o el comportamiento adecuado del transporte
de sedimentos a lo largo de la profundidad y a lo largo de la curva requiere un modelo numérico 3D. La Figura
7­7 y la Figura 7­8 proporcionan un ejemplo de esta distribución de flujo helicoidal.

Figura 7­5. Campo de velocidad promediado en profundidad calculado en el río Mississippi con transporte de vorticidad.

7­6
Machine Translated by Google

TD­41 Capítulo 7

Figura 7­6. Cambio en la velocidad calculada producida por el método de transporte de vorticidad.

Figura 7­7. Flujo helicoidal en el recodo de un río.

7­7
Machine Translated by Google

Capítulo 7 TD­41

Figura 7­8. Diferencia de velocidad entre la superficie y el fondo en el recodo de un río.

• Modelado de Retiro Selectivo: Estrategias de retiro selectivo, para lavado de reservorios,

desvíos de ríos, etc., se utilizan para optimizar las estrategias de gestión. La extracción selectiva
implica la eliminación de agua en ciertos lugares de la columna de agua. Los procesos
involucrados en este retiro son inherentemente de naturaleza 3D y no pueden replicarse ni
simularse utilizando modelos 2D o posiblemente 2D promediados lateralmente. Sin embargo, los
modelos 2D se pueden utilizar para reducir la cantidad de diseños alternativos que se simularán utilizando modelos 3D.

Una de las preguntas frecuentes de los ingenieros, planificadores y diseñadores hidráulicos es "¿se puede
utilizar un modelo 2D para diseñar compuertas y/o estructuras de aliviaderos?" La respuesta a esta pregunta
es específica del proyecto, compleja y dependiente de la experiencia. Por ejemplo, no se deben utilizar modelos
2D para diseñar compuertas y/o aliviaderos presurizados; sin embargo, se recomiendan encarecidamente los
modelos 2D para reducir las alternativas de diseño para las que se requiere modelado físico y/o 3D. En
general, para estas estructuras se recomiendan 3D­Non Hidrostático (3D­NH) y 3D­Non
Hidrostático Multifásico (3D­NHMP). Se recomienda el modelado físico de estas estructuras en
circunstancias en las que la estructura es excepcionalmente única y no existe conocimiento para medir el
rendimiento de los modelos 3D­NH/3D­NHMP.

7­8
Machine Translated by Google

TD­41 Capítulo 7

Modelado 3D de Estructuras

A continuación se presentan tres ejemplos para ilustrar la comparación de modelos 2D y 3D para tales esfuerzos de diseño. En
los siguientes tres ejemplos, la terminología 3D indica 3D­NH y 3D­NHMP. Los ejemplos siguientes suponen que los resultados
del modelo 3D son precisos.

• Modelado de flujo a través de compuertas: Modelado de flujo a través de compuertas para regulación de piscinas,
la gestión del flujo aguas abajo, etc. es una parte relativamente rutinaria de la gestión del agua y la hidráulica. Sin
embargo, modelar la relación entre la cantidad de flujo y la etapa aguas arriba no lo es. Este tipo específico de problema
se maneja mediante modelado 3D­NHMP y modelado físico cuando la complejidad del proyecto lo requiere.
Una pregunta común es "¿se puede utilizar el modelado 2D para guiar adecuadamente los criterios de selección o para
reducir la cantidad de diseños a simular usando 3D?" Considere un tramo recto de un río con 9 compuertas abiertas de
par en par, que se requieren para conducir el flujo bajo ciertas condiciones. A modo de comparación, la Figura 7­9
proporciona los resultados del modelado para este ejemplo mostrando los resultados de velocidad de 2D
(profundidad promediada) y 3D­NHMP (sin presurizar) en un plano vertical a través de los centros de las puertas.
Los resultados de velocidad 2D no imitan suficientemente los resultados de 3D­NHMP y la distribución del flujo a
través de las compuertas y, sin embargo, las zonas de recirculación aguas abajo de las compuertas se reproducen
adecuadamente (Figura 7­9).

Consideremos ahora el mismo ejemplo, excepto que la puerta central penetra 0,5 metros en la columna de agua. La
Figura 7­10 ilustra la comparación de velocidades entre los resultados 2D y las simulaciones presurizadas 3D­NHMP para
esta nueva situación en la que la puerta central penetra 0,5 metros en la columna de agua. Nuevamente, la
comparación indica que la velocidad promedio en profundidad del modelo 2D no imita suficientemente las velocidades
del modelo 3D­NHMP (Figura 7­10). Hay diferencias obvias en la zona presurizada, sin embargo, los resultados
fuera de esta área tampoco coinciden con la simulación 3D­NHMP. Este ejercicio de modelado, de un diseño
de puerta simple, indica que "sí, las simulaciones 2D se pueden usar para guiar los diseños que deben
simularse usando 3D­NHMP y los diseños que deben descartarse de la consideración". Sin embargo, los resultados del
modelo 2D “NO” deben usarse para los diseños finales de aberturas de puertas.

7­9
Machine Translated by Google

Capítulo 7 TD­41

Figura 7­9. Comparación de velocidad entre modelos 2D (profundidad promediada) y 3D­NHMP (sin presión) para
un tramo recto de un río con nueve compuertas del mismo tamaño. Colores escalados para
ilustrar patrones, no valores exactos.

7­10
Machine Translated by Google

TD­41 Capítulo 7

Figura 7­10. Comparación de velocidad promedio en profundidad entre modelos 2D y 3D­NHMP (presurizado) para
un tramo recto de un río con nueve compuertas donde la compuerta central penetra 0,5 metros en la
columna de agua. Colores escalados para ilustrar patrones, no valores exactos.

7­11
Machine Translated by Google

Capítulo 7 TD­41

• Diseño de aliviadero inclinado simple: Los aliviaderos son esenciales para la gestión del embalse.
estructuras y conducir agua aguas abajo del embalse. Un componente esencial del diseño de un aliviadero
es la determinación de las cargas hidráulicas y la duración del endurecimiento necesario para evitar la
socavación de la estructura. Las mejores prácticas estándar de la industria para el diseño de aliviaderos
se basan en la utilización de modelos físicos y/o modelos 3D­NHMP.
La disponibilidad del software numérico 3D­NHMP ha llevado a los ingenieros a utilizar modelos físicos
para validar el modelo numérico y luego utilizar los modelos numéricos para iterar sobre el diseño óptimo.
Los modelos 3D tienen un tiempo de respuesta más rápido, así como costos más bajos que los modelos
físicos, pero son varias veces más exigentes desde el punto de vista computacional que los modelos 2D.
Por lo tanto, sigue existiendo la tentación de utilizar modelos 2D para reducir aún más el tiempo dedicado al
diseño mediante el uso de modelos 2D para eliminar diseños obviamente ineficaces.

La Figura 7­11 muestra el diseño de un embalse junto con la ubicación del aliviadero y el canal de salida.
El aliviadero debe conducir una Inundación Máxima Probable (PMF) de 6.179 m3 /seg, y el caudal de
descarga de diseño es de 158 metros. La Figura 7­12 muestra los resultados de velocidad de un
modelo 3D, y la Figura 7­11 y la Figura 7­13 muestran los resultados de un modelo 2D. Tanto el modelo
3D como el 2D muestran un comportamiento de velocidad similar. Tanto el modelo 3D como el 2D simulan
ubicaciones similares para el salto hidráulico; sin embargo, el modelo 2D no disipa la misma energía en
comparación con el modelo 3D.

Además, el modelo 3D simula velocidades entre 1 y 5 m/s después del salto (Figura 7­12), mientras
que el modelo 2D simula entre 8 y 10 m/s (Figura 7­14). Esto indica que si se utiliza el modelo 2D para
diseñar el aliviadero, se requeriría una plataforma endurecida más larga, cuando, como lo muestran los
resultados 3D, la energía se habría disipado dentro de una plataforma significativamente más corta. Este
resultado es esperado debido al hecho de que los modelos 2D no tienen todos los mecanismos a través de
los cuales se disipa la energía, es decir, la mezcla vertical y la interacción entre aire/agua/estructura. Por
lo tanto, dependiendo del régimen de flujo, el control aguas abajo, los parámetros de diseño necesarios
y otras características, el modelo 2D podría ser capaz de replicar suficientemente la física del sistema,
pero no es probable. Las mejores prácticas indican que el modelo 2D debe validarse a un modelo 3D y
solo entonces usarse para guiar la selección de diseños que deben simularse usando 3D.

Figura 7­11. Batimetría y velocidad 2D para un aliviadero inclinado que debe transmitir un Máximo Probable
Inundación (PMF) de 6.179 m3 /seg.

7­12
Machine Translated by Google

TD­41 Capítulo 7

Figura 7­12. Aliviadero inclinado (que debe transmitir un PMF de 6.179 m3 /seg), velocidad 3D (los colores representan
velocidad) y la superficie del agua (el espesor representa la superficie del agua).

Figura 7­13. Batimetría, velocidad 2D y perfil de la superficie del agua para el aliviadero inclinado (que debe
transportar un PMF de 6.179 m3 /seg).

• Diseño de aliviadero con pendiente escalonada: Otro diseño común de un aliviadero implica la
creación de escalones para la disipación de energía antes del estanque amortiguador. La Figura
7­14 muestra un ejemplo de aliviadero escalonado. Este ejemplo de aliviadero, construido de hormigón, debe pasar

7­13
Machine Translated by Google

Capítulo 7 TD­41

una Inundación Máxima Probable de ~3,145 m3 /s, y el control de colada está en ~437 metros.
La Figura 7­15 muestra el comportamiento hidráulico simulado del modelo 3D­NHMP del aliviadero y el
estanque amortiguador, y la Figura 7­14 muestra lo mismo para un modelo 2D.

Figura 7­14. Batimetría para un aliviadero escalonado que deberá pasar un PMF de aproximadamente 3.145 m3 /s.

Un vistazo rápido indica que ambos modelos proporcionan una hidráulica similar, pero en una
inspección más cercana se observa que la velocidad máxima simulada por el modelo 3D­NHMP es de 10
m/s (Figura 7­15), en comparación con ~15 m/s calculada por el modelo 2D. modelo (Figura 7­16).
Se espera esta diferencia en la velocidad máxima, ya que el modelo 3D­NHMP incurre en pérdida de
energía debido a que el agua cae pendiente abajo sobre los escalones, pero el modelo 2D no puede
explicar estas pérdidas.

Figura 7­15. Resultados de la hidráulica simulada 3D­NHMP para un aliviadero escalonado (que debe pasar un
PMF de aproximadamente 3.145 m3 /s).

Cabe señalar que si los pasos son menos severos, los resultados 2D pueden proporcionar resultados
similares a los del 3D­NHMP. La simulación 3D­NHMP muestra que existe un chorro de agua bajo la
superficie en el control del agua de cola, y el chorro se extiende hacia la superficie aguas abajo, el modelo 2D

7­14
Machine Translated by Google

TD­41 Capítulo 7

El modelo es incapaz de simular este comportamiento. El modelo 2D muestra la existencia

de un salto (Figura 7­16), que no se muestra en la simulación 3D­NHMP. En esta ubicación
el modelo 3D­NHMP muestra la existencia de interacción aire­agua y arrastre de aire, el
modelo 2D es incapaz de simular este comportamiento.

Figura 7­16. Resultados de la simulación hidráulica 2D para un aliviadero escalonado (que debe pasar un PMF de
aproximadamente 3.145 m3 /s) para diferentes rugosidades de Manning (n). El espesor indica profundidad.

7­15
Machine Translated by Google

Capítulo 7 TD­41

Un buen modelador numérico sabrá que el modelo 2D es incapaz de contabilizar todas las pérdidas en este
sistema, y podría emprenderse una modificación de los valores de rugosidad. Esta decisión es aceptable
y está justificada dentro de lo razonable, entendiendo que el modelado 2D se está realizando como un medio
para reducir la cantidad de modelado físico y/o 3D­NHMP, y no para el diseño final. La Figura 7­16, panel
inferior, muestra los resultados con la rugosidad de Manning aumentada de 0,015 a 0,025. Tenga en cuenta
que las velocidades coinciden con las reportadas por el modelo 3D­NHMP, el comportamiento en el pie es
erróneo, pero el comportamiento promedio en el aliviadero, no en la cuenca amortiguadora, es correcto.

Cabe señalar que este aumento en la rugosidad de Manning no está físicamente justificado y no debe
alentarse.

Los tres ejemplos presentados anteriormente ilustran que los modelos 2D simulan un comportamiento en general
inconsistente con los modelos 3D­NHMP. Un buen modelador numérico reconocerá que el modelo 2D es incapaz de dar
cuenta, implícitamente, de todas las pérdidas en este sistema. En este caso es aceptable aumentar la rugosidad del
modelo, dentro de lo razonable, para obtener una mejor comprensión del sistema. Por lo tanto, se recomiendan modelos
3D­NHMP para dicho diseño, complementados con modelado físico si es necesario. Sin embargo, se recomienda que
los modelos 2D se puedan utilizar como herramienta de eliminación de diseño preliminar para reducir potencialmente
la cantidad de modelos físicos y/o simulaciones 3D­NHMP. Sin embargo, los modelos 2D no deben usarse para diseños
finales de este tipo de estructuras (como se muestra arriba).

En las secciones siguientes, 3D indica 3D­Hidrostático a menos que se especifique lo contrario.

Modelo de desarrollo
Si bien el desarrollo de un modelo hidráulico es similar, independientemente de si el modelo es 2D o 3D, también existen
muchas diferencias significativas en el proceso de desarrollo del modelo. También existen diferencias significativas en
el desarrollo de modelos entre 3D hidrostático y 3D no hidrostático. Esta sección analiza el proceso de
desarrollo del modelo y destaca las diferencias entre los dos enfoques: hidrostático 2D y 3D. Además, se discuten las
ventajas y desventajas de los enfoques de modelado 2D y 3D.

Datos del terreno/batimétricos

Como se analizó anteriormente, todo desarrollo de modelos hidráulicos e hidrodinámicos requiere datos del terreno. Los
modelos 3D requieren que se construya un modelo batimétrico o de terreno completo que incluya datos para todo el
dominio de modelado. Dependiendo de la estrategia de transformación de malla 3D utilizada, es posible que sea
necesaria una modificación que incluya el suavizado de la batimetría. La siguiente sección proporciona ejemplos de
varias transformaciones de malla utilizadas en el modelado 3D. Para obtener detalles sobre las necesidades del terreno
para enfoques de modelado 2D, consulte el Capítulo 4.

Transformación de malla

Los modelos numéricos 3D requieren un método adecuado para rastrear la superficie libre transitoria.
Históricamente, se ha utilizado una variedad de descripciones de coordenadas, siendo la más popular:

• Transformada sigma ( ) (Figura 7­17A),

7­16
Machine Translated by Google

TD­41 Capítulo 7

• Transformada del plano Z (Figura 7­17B) y •

Euleriana lagrangiana arbitraria (ALE) (Figura 7­17C).

Cada uno de estos tiene ventajas y desventajas. La cuadrícula de transformación Z es ideal para simulaciones en
3D en las que un fluido más denso debe permanecer atrapado dentro del canal profundo, la cuadrícula de
transformación es ideal para simulaciones en las que un fluido más denso debe caer en un canal profundo y la
cuadrícula de ALE es propicia para ambos mencionados. condiciones. Aunque la cuadrícula ALE es ventajosa
para la mayoría de las condiciones de modelado, en general requiere una resolución más alta que la
cuadrícula con transformada z y .

Figura 7­17. Estrategias de mallado para 3D.

Una breve nota sobre el mallado 3D­no hidrostático (3D­NH)

El mallado para modelos 3D­NH es significativamente diferente del mallado para modelos 3D­Hidrostáticos.
Los modelos 3D­NH en general no se limitan a columnas verticales y pueden tener un diseño de malla completamente
desestructurado. Esto se debe al hecho de que los modelos 3D­NH resuelven la presión en lugar de asumir una
distribución hidrostática.

Discretización del modelo

En general, los modelos 2D y 3D no requieren conocimiento de las trayectorias exactas del flujo del agua antes
de diseñar el modelo. Sin embargo, el conocimiento de las rutas de flujo puede resultar muy útil para
desarrollar una malla más detallada para el modelado 2D y 3D. Por ejemplo, para modelar con precisión

7­17
Machine Translated by Google

Capítulo 7 TD­41

flujo de salinidad dentro y fuera de un estuario, la información sobre la ubicación de la estratificación debe
conocerse a priori para proporcionar suficiente resolución vertical para capturar esta estratificación y para evitar la
resolución vertical en áreas que no experimentan esta estratificación. Los modelos 3D son
computacionalmente costosos y la eficiencia en la resolución es esencial.

Una vez definidas las extensiones del modelo, los conceptos básicos del desarrollo de un modelo 3D son similares
al modelo 2D en resolución horizontal. Se requiere una mayor resolución en términos de tamaño más pequeño
de celda/elemento para ubicaciones de interés o ubicaciones donde hay un cambio rápido en la hidrodinámica o el
transporte. Estas regiones incluyen a menudo estructuras, curvas, canales de navegación, ubicaciones de
afluencia, contracciones o expansiones y similares.

Independientemente de la rapidez del cambio, los modelos 3D requieren que la resolución horizontal sea suficiente
para capturar características batimétricas esenciales para simular el campo de flujo/transporte 3D. La
Figura 7­18 proporciona un ejemplo de este aumento de resolución en una ubicación donde un canal de
navegación se divide en dos y se requiere un aumento de resolución para proporcionar una resolución
vertical adecuada cuando se genera información en capas.

Figura 7­18. Mallado horizontal para modelos 3D, el color verde a rojo indica las características batimétricas (de más
profundo a menos profundo).

Además de las características de flujo de interés, los recursos computacionales disponibles también influyen en
la resolución horizontal. Esta influencia es especialmente cierta para los códigos hidrostáticos 3D; en
consecuencia, estos modelos 3D requieren que se mantenga una estructura de columnas, y cada celda/
elemento horizontal requerirá una descripción de capas verticales. Sin embargo, una resolución horizontal
innecesariamente alta aumentará los costos de cálculo y los tiempos de ejecución de la simulación 3D.

7­18
Machine Translated by Google

TD­41 Capítulo 7

La información de capas verticales 3D y la manera en que se disponen las capas depende en gran medida del tipo
de descripción de coordenadas utilizada y de la naturaleza de la simulación. Por ejemplo, si se utiliza la
transformación , entonces el número de capas en la vertical se llevará a cabo en toda la extensión de la simulación
y será necesario suavizar cualquier cambio abrupto de pendiente o discontinuidad en la batimetría. Sin
embargo, si se utiliza la transformación z, entonces se debe determinar el espaciado vertical de las capas
(teniendo en cuenta que esta transformación modifica la batimetría a escalón). De lo contrario, si se utiliza una
malla ALE 3D, entonces las capas verticales deben mantener los nodos en niveles de contorno constantes. En
la Figura 7­19 se presenta un ejemplo de capas ALE. Tenga en cuenta que en una capa los nodos se disponen
aproximadamente al mismo nivel de contorno.

Figura 7­19. Ejemplo de mallado vertical Arbitrario­Lagrangiano Euleriano (ALE) para modelos 3D.

Definición de coeficientes de rugosidad

La definición de valores de fricción (coeficientes de rugosidad, altura de rugosidad, etc.) es muy similar entre 2D y
3D. Los coeficientes de rugosidad, en general, se definen en función de los datos disponibles, como el uso/cobertura
del suelo, el tipo de vegetación, el comportamiento de los sedimentos, etc. La malla/cuadrícula 3D inicial debe
construirse teniendo en cuenta los distintos tipos de dominios/materiales/materiales de rugosidad. regiones que podrían
estar en el dominio modelado, esto es esencial ya que crear dominios/materiales/regiones de rugosidad adicionales en
una cuadrícula/malla 3D es un proceso tedioso y que requiere mucho tiempo. Una vez que se completan las
simulaciones 3D iniciales, el proceso de calibración y validación se puede utilizar para guiar la modificación de los
valores de rugosidad iniciales cuando sea necesario.

Definición de cierre de turbulencia, viscosidad de remolino de fondo y valores de difusión

Los modelos 3D, en comparación con los modelos 2D, proporcionan medios sólidos para calcular los valores de cierre
de turbulencia. Estos valores de cierre de turbulencia se utilizan para tener en cuenta los procesos que surgen como
resultado de cambios caóticos y/o gradientes de presión y velocidad. La turbulencia aumenta la presión hidráulica o

7­19
Machine Translated by Google

Capítulo 7 TD­41

Pérdidas hidrodinámicas en el campo de flujo, más a menudo a través de la generación de vórtices que pueden
interactuar o no entre sí. En el modelado numérico 2D y 3D, los modelos de cierre de turbulencia se utilizan para
tener en cuenta las pérdidas en el campo de flujo a escalas menores que la resolución de malla/celda.
La descripción matemática de los esquemas de cierre de turbulencias está involucrada, y el lector interesado
puede consultar a Savant (2015), Blumberg et al., (1992), Mellor y Yamada (1982) y Henderson­Sellers (1982)
para una descripción detallada. discusión. La Tabla 7­1 presenta pautas generales para la utilización del cierre de
turbulencia en modelos 3D.

Tabla 7­1. Esquemas sugeridos para el cierre de turbulencias

Turbulencia vertical Turbulencia horizontal
Condición de flujo Cierre Cierre
Flujo estuarino estratificado para
Mellor­Yamada Nivel 2 Smagorinski
Sistemas de micro mareas
Flujo estuarino estratificado para Mellor­Yamada Nivel 2.5, k­ Mellor­Yamada Nivel 2.5, k­
Macrosistemas de mareas Smagorinski
Caudales fluviales sin
Mellor­Yamada Nivel 2 Smagorinski
estructuras hidráulicas
Flujos fluviales con Mellor­Yamada Nivel 2.5, k­ Mellor­Yamada Nivel 2.5, k­
Estructuras Hidráulicas Smagorinski
Embalses con
k­ k­
Retiro selectivo
Embalses/Lagos estratificados Mellor­Yamada Nivel 2 Smagorinski

Modelos de cierre de turbulencia de nivel superior o de 2 ecuaciones, como Mellor­Yamada Nivel 2.5 y k­
implican la solución de ecuaciones de transporte para la generación y disipación de turbulencias. Como
consecuencia, el costo de cálculo aumenta en comparación con los modelos de 1 ecuación como el nivel 2 de
Mellor­Yamada y otros modelos más simples.

La mayoría de los modelos numéricos 3D ofrecen la opción de proporcionar valores de fondo para el tensor
de viscosidad turbulenta. Estos valores de fondo suelen depender de la simulación y del sistema, y son difíciles
de determinar. Como consecuencia, estos valores proporcionan un método conveniente para estabilizar la
simulación numérica y como herramientas de calibración. Se recomienda que estos valores se mantengan en
niveles que estabilicen la simulación, pero no más altos a menos que las observaciones lo justifiquen.

Para simulaciones 3D que involucran el transporte de constituyentes, se requiere la especificación de valores de

difusión de fondo y turbulencia. Además, las sugerencias proporcionadas en la Tabla 7­1 también son válidas
para la difusión turbulenta. Las ecuaciones de transporte tanto en 2D como en 3D incorporan difusión,
turbulenta y de fondo, como términos que redistribuyen la masa. Como consecuencia, se debe tener extremo
cuidado al especificar los valores de difusión de fondo. Los valores de difusión de fondo son extremadamente
tediosos de determinar y son especificados por el usuario y, por lo tanto, proporcionan una manera conveniente de
hacer coincidir las observaciones. La utilización de los valores de difusión de fondo de esta manera podría
permitir que la simulación coincida con los valores observados bajo una condición, pero con toda probabilidad
fallará la validación bajo un conjunto diferente de condiciones de flujo. Por lo tanto, para un modelo de pronóstico,
los valores de fondo no deben usarse como herramientas de calibración y deben establecerse en un valor que
estabilice el transporte, pero no en un valor superior, a menos que las observaciones lo justifiquen.

7­20
Machine Translated by Google

TD­41 Capítulo 7

Estructuras hidráulicas y otras

La definición y el modelado de estructuras hidráulicas es similar en los enfoques de modelado 2D y 3D, con la
única diferencia en la estructura de capas verticales y la presencia de cálculos de turbulencia que
inherentemente tienen en cuenta las pérdidas de expansión y contracción. Las estructuras no presurizadas,
como puentes no sumergidos, alcantarillas, etc., se pueden modelar directamente utilizando las
ecuaciones de movimientos; para estructuras presurizadas se deben tener en cuenta consideraciones de
presión adicionales y divisiones de flujo.

Al modelar estructuras hidráulicas, los modelos 3D pueden proporcionar más información hidráulica sobre la
socavación de pilares de puentes y otros fenómenos en comparación con las simulaciones 2D. Sin embargo,
debido a la presencia de aceleraciones verticales, los modelos 2D y 3D hidrostáticos son inapropiados para
información de flujo o transporte de campo cercano, como sedimentos. En cambio, para estas simulaciones,
se deben utilizar modelos 3D no hidrostáticos.

Condiciones de borde

La especificación de límites es similar en modelos 2D y 3D. Los modelos 2D en general especifican límites en
bordes/células externos y bordes/células internos como entradas, salidas, etapas, precipitaciones, evaporación
e infiltración. Los modelos 3D permiten la especificación de los mismos límites, pero estas especificaciones
pueden variar en profundidad si se desea. Los límites de flujo supercrítico, comunes en los modelos
2D, no se utilizan en los modelos 3D debido a las deformaciones extremas de la superficie. El transporte baroclínico,
como la salinidad y la temperatura, puede provocar que los límites 3D tengan flujos en direcciones opuestas.
Este tipo de transporte se observa comúnmente en lugares con límites de elevación de la superficie del agua.

Calibración del modelo

El proceso de calibración del modelo es uno de los pasos más importantes en el desarrollo de un modelo
hidráulico. Se requiere la calibración de cualquier modelo hidráulico para comprender si el modelo es capaz
de reproducir flujos/transporte pasados y si tiene alguna posibilidad de predecir eventos de flujo/transporte futuros.
El proceso de calibración también permite una mayor comprensión de la sensibilidad del modelo a los
datos, las fuerzas de fricción y otros coeficientes empíricos. Un modelo que no está calibrado es en realidad sólo
un experimento numérico. De hecho, un modelo no calibrado puede o no estar siquiera cerca de reproducir
elevaciones, flujos y transporte realistas de la superficie del agua en todo el sistema. Nunca se deben
utilizar modelos 3D no calibrados para tomar decisiones de ingeniería, ni siquiera en situaciones de
emergencia. Esta afirmación se debe a que el esfuerzo necesario para crear un modelo 3D requerirá generalmente
escalas de tiempo mayores a las que permitirá la emergencia, y la complejidad del modelado 3D se presta a
cometer errores.

El proceso de calibrar un modelo es muy similar entre los enfoques de modelado 2D y 3D; sin embargo, la
calibración del modelo 3D puede resultar más difícil y llevar más tiempo. Los datos necesarios para realizar la
calibración de un modelo son prácticamente los mismos, excepto que una calibración más detallada de los
modelos 3D requiere mediciones detalladas de la distribución de velocidad para garantizar la precisión de
los perfiles de velocidad vertical. Para flujos 3D dependientes de la densidad, también se requiere la calibración
del transporte y puede ser muy compleja. Esta complejidad de calibración se debe a las interacciones de retroalimentación entre

7­21
Machine Translated by Google

Capítulo 7 TD­41

hidrodinámica y transporte (es decir, la hidrodinámica influye en el transporte y el transporte, a cambio, influye en
la hidrodinámica a través de impactos de densidad).

Los datos necesarios para calibrar y validar un modelo 3D son:

• Elevaciones de la superficie del

agua, • Perfiles de velocidad horizontal y vertical, •
Divisiones de flujo, y •
Perfiles verticales de los componentes transportados.

Problemas de tiempo y costos

El tiempo para desarrollar un modelo con el enfoque de modelado 2D o 3D puede variar, dependiendo del tipo y
propósito del modelo. Por ejemplo, si es necesario desarrollar un modelo “rápido y sucio” durante una
emergencia por inundación, entonces es mucho más rápido diseñar un polígono de área de flujo 2D, establecer un
tamaño básico de celda/elemento, adjuntar algunas condiciones de contorno y calcular. Por otro lado, crear modelos
3D puede llevar mucho tiempo y ser tedioso. La creación de un modelo 3D requiere la creación de un dominio
2D, así como información de capas verticales y, por lo tanto, su creación siempre requiere más tiempo. Estimar
los parámetros base, tales como: coeficientes de rugosidad, parámetros de fondo y otros coeficientes
empíricos necesarios para cada uno de los enfoques de modelado, requiere aproximadamente la misma
cantidad de tiempo.

En general, el costo computacional de los modelos 2D, que resuelven las ecuaciones de continuidad, x
e y de momento para u, v y h (es decir, 3 grados de libertad por celda/elemento computacional) es menor que el
de los modelos 3D. Esta diferencia de costo computacional se debe a que los modelos 3D resuelven la
continuidad, el momento x, y y el supuesto hidrostático para resolver u, v, w y h (es decir, 4 grados de libertad
por celda/elemento computacional); y, por lo tanto, siempre son más costosos desde el punto de
vista computacional de resolver en comparación con 2D. La adición de capas verticales aumenta aún más el gasto
computacional.

Los avances recientes en modelos 2D han permitido la simulación de hidrodinámica a escala decenal y a
escala continental; sin embargo, incluso los modelos 3D más avanzados están limitados a escalas anuales y
escalas espaciales significativamente más pequeñas, del orden de decenas de millas.

Por otro lado, los modelos 3D requieren un hardware computacionalmente más potente que el de los modelos 2D.
Por lo tanto, la mayoría de las veces, los modelos 3D se ejecutan en sistemas multiprocesador, como computadoras
de alto rendimiento (HPC) o computación en la nube. Las máquinas de escritorio se pueden utilizar para modelos
3D relativamente pequeños dependiendo de la resolución del modelo, el recuento de celdas/elementos y las escalas
espaciales/temporales de simulación.

Dependiendo del tipo de modelo, los modelos 3D pueden proporcionar resultados de profundidad, presión,
velocidades, concentraciones, etc. En consecuencia, estos resultados pueden ser del orden de gigabytes, si no
mayores, y son complejos de procesar. Sin embargo, las herramientas de visualización para modelos 3D
son primitivas en comparación con las herramientas de visualización para modelos 2D. Como consecuencia, los
ingenieros crean formas ad hoc de extraer datos generados por modelos 3D y/o extraer información de capas y
visualizar cada capa o columna como 2D. Estas prácticas ad hoc añaden un coste de tiempo adicional en el
posprocesamiento de los resultados.

7­22
Machine Translated by Google

TD­41 Capítulo 7

Alternativamente, el costo computacional y de tiempo de ejecutar un modelo 3D se puede reducir si se utiliza un

modelo 2D para "iniciar en caliente" el modelo 3D, pero esto requiere la existencia de resultados 2D o la creación
de una simulación 2D.

Conocimientos, habilidades y habilidades del modelador

Ya sea que realice modelado 2D o 3D, el modelador debe tener una buena experiencia en hidráulica/
hidrodinámica, hidráulica alrededor de estructuras, transporte y soluciones numéricas de ecuaciones no lineales.
Además, no importa qué modelo de computadora y enfoque de modelado se seleccione, el modelador debe tomar
clases sobre cómo usar ese software específico de manera efectiva para el modelado 2D y/o 3D. El autoestudio es
excelente para comenzar a usar modelos 2D, pero el conocimiento adicional adquirido al tomar una clase de
los desarrolladores/expertos en el uso del software a menudo resultará invaluable. Por el contrario, los modelos 3D
son extremadamente complejos y el autoestudio es de menor utilidad. Esta discrepancia se debe a que aunque
todos los modelos 3D resuelven las mismas ecuaciones de movimiento y transporte, la discretización en horizontal y
vertical y las buenas prácticas son específicas del modelo. Además, la aplicación y utilización de las opciones de
turbulencia varían entre modelos y dentro del mismo modelo dependiendo de ciertas condiciones.

Por lo tanto, el modelador debe familiarizarse con las técnicas para preprocesar y posprocesar la malla/cuadrícula 3D
para satisfacer las buenas prácticas de estratificación para el modelo que se utiliza y para visualizar los resultados de la
simulación.

Además, el modelador necesita realizar estudios detallados de sistemas en los que tenga buena información
batimétrica, de estructura y datos históricos para realizar calibraciones del modelo. Si el modelador es nuevo
en el modelaje, siempre debe buscar orientación y ayuda de modeladores experimentados. Además, todos
los modelos deben ser revisados por expertos independientes para garantizar que el trabajo realizado se haya
realizado con una opción de modelado adecuada (2D o 3D); aplicado correctamente; contiene la cantidad adecuada
de detalles para la ubicación y el tipo de estudio determinados; fue calibrado efectivamente y está reproduciendo
resultados hidráulicos razonables para los eventos y alternativas que se están modelando.

Resumen de ventajas y desventajas del modelado 2D y 3D

Existen claras ventajas y desventajas tanto para el enfoque de modelado 2D como para el 3D. Esta
sección del documento resume las ventajas y desventajas que se han descrito anteriormente en este documento.
La discusión se basa en las ventajas y desventajas desde la perspectiva del modelado 2D versus el modelado 3D.
Aunque las ventajas y desventajas del modelado 2D se enumeraron anteriormente, aquí se presentarán desde
la perspectiva de comparación con el modelado 3D.

Ventajas del modelado 2D:

No es necesario conocer la trayectoria del flujo del agua, para todos los eventos, para desarrollar el
modelo. Sin embargo, es necesario definir correctamente la extensión del dominio.

La dirección del flujo puede cambiar durante el evento. El agua puede moverse en cualquier dirección,
basado en la energía y el impulso del flujo.

7­23
Machine Translated by Google

Capítulo 7 TD­41

La velocidad, el momento y la dirección del flujo se tienen en cuenta con mayor precisión. De hecho, esta precisión es
especialmente cierta para el flujo que pasa por carreteras, diques, barreras, estructuras, alrededor de curvas y en cruces/
divisiones de flujo. Además, los modelos 2D se pueden utilizar para analizar zonas de remolinos 2D dentro del
campo de flujo. Alrededor de las curvas, los modelos 2D producen elevaciones precisas de la superficie del agua, pero las
distribuciones de velocidad pueden ser erróneas debido a la existencia de flujo helicoidal.

Las pérdidas de energía y fuerza debidas a contracciones y expansiones, etc., se contabilizan directamente y no
requieren coeficientes empíricos, mayor rugosidad ni áreas de flujo ineficaces definidas por el usuario.

El mapeo del área inundada, así como las velocidades y los peligros de inundación (profundidad x
velocidad) es exacta.

El modelado detallado de estructuras hidráulicas, en un enfoque de modelado 2D completo, puede proporcionar más
información sobre la distribución del flujo que se acerca, atraviesa y sale de una estructura.

Tiempos de cómputo rápidos, así como tiempo de obtención de resultados.

Desventajas del modelo 2D:

Es posible que se produzca una falsa sensación de confianza en los resultados del modelo 2D en escenarios en los que el modelo 3D
Los efectos no son obvios, como el flujo alrededor de curvas, el flujo alrededor de estructuras, la simulación de embalses/
lagos, etc.

No se puede obtener información sobre la estructura vertical del flujo a partir de modelos 2D.

Ventajas del modelado 3D:

Los modelos 3D proporcionan una mejor representación de la hidráulica/hidrodinámica de un sistema.

Los modelos 3D proporcionan la única manera de modelar de manera realista sistemas hidráulicos que se estratifican y el
transporte resultante. Estos incluyen sistemas como embalses, lagos y estuarios.

Los modelos 3D pueden proporcionar un comportamiento de velocidad preciso para el flujo alrededor de las curvas y
replicar con precisión el sedimento resultante u otro comportamiento de transporte.

Los modelos 3D proporcionan una mejor representación de la hidráulica alrededor de las estructuras. Sin embargo, para
representar con precisión la hidráulica alrededor de las estructuras, se deben utilizar modelos 3D no hidrostáticos.

Los modelos 3D proporcionan un comportamiento hidráulico preciso para sistemas hidráulicos que incluyen
retirada selectiva.

Las pérdidas de energía en forma de remolinos, las pérdidas por expansión/contracción, así como las pérdidas debidas a
obstrucciones sumergidas, como presas, se pueden replicar con precisión mediante modelos 3D.

En cierto modo, los modelos 3D tienen una base más física que los modelos 2D. Los parámetros como los coeficientes de
viscosidad de remolino que el usuario especifica en los modelos 2D se calculan internamente mediante modelos 3D.

Desventajas del modelado 3D:

7­24
Machine Translated by Google

TD­41 Capítulo 7

Se debe conocer la trayectoria del flujo del agua, para todos los eventos, para desarrollar la dirección vertical.
Estructura de capas en el modelo para modelos hidrostáticos 3D. El

proceso para adquirir datos adecuados, como viscosidad de fondo, etc., para el modelado 3D requiere pruebas de campo/
recopilación de datos.

El desarrollo de una malla/cuadrícula 3D generalmente requiere la descripción inicial del dominio.

en 2D.

El desarrollo de mallas en un modelo 3D puede ser complicado y puede requerir iteraciones. La naturaleza

compleja de los modelos 3D aumenta la probabilidad de inestabilidades. Esta inestabilidad requiere un conocimiento
profundo del código y puede requerir reinicios frecuentes del modelo, ya sea mediante arranque en caliente o
modificando los parámetros.

El tiempo y el coste de cálculo necesarios para el 3D pueden ser extremos y prohibitivos.

Se requiere un posprocesamiento de los resultados para obtener datos utilizables. La visualización de mallas/
cuadrículas y resultados en 3D puede requerir capacitación y conocimientos adicionales de software de
posprocesamiento.

La decisión entre la elección de modelos 2D o 3D debe mantener la fidelidad a los objetivos del proyecto que se
investiga y el principio del "camino más simple y técnicamente sólido para lograr los objetivos". Por ejemplo, si la
extensión de las inundaciones es de interés para un río que incluye curvas u otras características 3D, un modelo 2D
podría ser suficiente aunque la distribución de velocidad sea incorrecta; por otro lado, si el comportamiento de los
sedimentos es de interés, un modelo 3D podría ser más apropiado. Otra posibilidad es que el mismo sistema requiera
un modelo 2D o un modelo 3D dependiendo del sistema hidráulico que se esté investigando. El ingeniero debe
aprovechar el conocimiento adecuado del sistema para alcanzar con éxito los objetivos del proyecto.

7­25
Machine Translated by Google
Machine Translated by Google

TD­41 Capítulo 8

Capítulo 8

Modelos físicos hidráulicos

Las constantes mejoras en las técnicas numéricas y la tecnología informática durante el último medio siglo han
ampliado enormemente la gama de aplicaciones prácticas de modelado numérico.
No obstante, existen muchos problemas prácticos de ingeniería hidráulica para los cuales los modelos hidráulicos
físicos a escala siguen siendo el único método de solución o el más rentable.
La revolución digital que permitió avances en el modelado numérico también ha ampliado las opciones para controlar y
medir con precisión los flujos y fenómenos relacionados en modelos a escala. Por tanto, también se está ampliando la
gama de aplicaciones potenciales y la fidelidad de los modelos físicos hidráulicos.
Además, existen problemas para los cuales puede ser apropiada una combinación de modelado numérico y físico (es
decir, modelado híbrido).

Este documento no presenta orientación detallada sobre el desarrollo, uso y limitaciones de los modelos hidráulicos
físicos a escala. Se recomienda consultar con un experto en la materia apropiado al principio del proceso de estudio para
cualquier proyecto donde la aplicación de un modelo físico a escala pueda ser apropiada.

Las áreas generales de aplicabilidad incluyen:

• Flujos internos a estructuras hidráulicas, incluyendo: • Sistemas

de llenado y vaciado de esclusas.

• Obras de salida

• Unidades hidroeléctricas •

Determinación de fuerzas que actúan sobre las estructuras •

Canales de alta velocidad (particularmente para geometría compleja de canales/estructuras o flujos cerca de
profundidad crítica) •

Navegación de poco calado en las proximidades de estructuras •

Protección local contra la socavación (particularmente diseño de

aliviaderos) • Socavación alrededor de estructuras

atípicas pilares de

puentes • Paso de peces • Entradas de bombas (según lo especificado por normas formales)

Los modelos físicos a escala pueden clasificarse como modelos de procesos, utilizados para investigar procesos
físicos genéricos, o como modelos de diseño, utilizados para investigar aspectos específicos de una alternativa
propuesta o una falla experimentada. Los modelos de procesos son una fuente crítica de información para los avances
científicos en los campos de la mecánica de fluidos ambientales, sedimentación, ingeniería fluvial, etc. y para el
desarrollo y verificación de modelos numéricos. Los modelos de diseño son una herramienta ampliamente aceptada para
explorar alternativas de diseño y optimizar costos y rendimiento.

8­1
Machine Translated by Google

Capítulo 8 TD­41

al mismo tiempo que reduce la incertidumbre. Los modelos de diseño también pueden proporcionar datos
de validación para modelos numéricos a escala y así reducir la incertidumbre en los modelos numéricos a
escala prototipo. Para estructuras hidráulicas y otros sistemas hidráulicos complejos, los ahorros en costos de
construcción, operaciones y mantenimiento y los beneficios de las mejoras en el desempeño pueden exceder el
costo de un modelo de diseño por márgenes significativos (ASCE, 2000).

8­2
Machine Translated by Google

TD­41 Capítulo 9

Capítulo 9

Resumen

El modelado numérico de sistemas hidráulicos naturales y artificiales ha sido una parte integral de la gestión
desde la llegada de la informática moderna. Los primeros días del modelado numérico vieron la creación y amplia
utilización de modelos numéricos simplificados que se restringían a calcular propiedades del flujo a
escala gruesa, como los tiempos de residencia basados en el flujo y el volumen.
A medida que los recursos computacionales estuvieron disponibles, los modelos numéricos aumentaron en
complejidad para incluir simulaciones multidimensionales. Estas simulaciones quedaron relegadas a pequeñas
escalas temporales y espaciales, del orden de horas y millas respectivamente. No fue hasta la disponibilidad
de modelos 1D robustos que el modelado numérico de sistemas hidráulicos se volvió rutinario y lo
suficientemente extendido como para ser utilizado en la planificación y la gestión.

Los recursos computacionales sufrieron una revolución en la última década del siglo XX y el modelado
numérico multidimensional se convirtió en un facilitador realista para diseñadores, planificadores y gerentes
junto con el modelado 1D. Los avances recientes en computación de alto rendimiento (HPC), unidades de
procesamiento gráfico (GPU), computación en la nube y computación de escritorio han permitido simulaciones
numéricas multidimensionales de sistemas hidráulicos en detalle. La disponibilidad de 1D y multidimensionales
(2D y 3D) ha llevado a la pregunta de "¿cuáles usar?" o más apropiadamente "¿cuáles no usar?"

La elección del tipo de modelo numérico, en la mayoría de los casos, depende de la habilidad, experiencia,
conocimiento y preferencia del modelador. Un modelador experto puede obtener información útil de todo
tipo de modelos numéricos; sin embargo, la selección juiciosa del modelo numérico apropiado promete importantes
ahorros de tiempo y costos. Además, la elección de un modelo numérico depende del propósito para el cual se
investiga el sistema.
El mismo sistema podría requerir un modelo numérico diferente basado en los objetivos del proyecto, el plazo,
los recursos computacionales y los requisitos de datos. Por ejemplo, la determinación de la extensión de las
inundaciones se puede realizar con modelos 1D, 2D y 3D, pero el tiempo y los costos computacionales
pueden indicar que se debe utilizar un modelo 1D porque los modelos 2D/3D serán relativamente prohibitivos en
tiempo y cálculo. El mismo sistema, si se investiga en busca de velocidades detalladas, requerirá la utilización
de un modelo 2D/3D en lugar del modelo 1D.

La falta de disponibilidad de datos suele ser una mala excusa para disuadir el uso de un modelo si los objetivos
y la física del sistema indican que se requiere un modelo 2D o 3D. En tales casos, se deben adquirir datos
adicionales para simular con precisión el sistema y lograr los objetivos del proyecto.
Este documento ha presentado pautas y sugerencias sobre dónde son válidas las diversas técnicas numéricas
y en qué medida. La Tabla 9­1 presenta varios sistemas que comúnmente simulan los ingenieros hidráulicos y el
tipo de modelo más simple de usar.

9­1
Machine Translated by Google

Capítulo 9 TD­41

Tabla 9­1. Modelado recomendado para varios sistemas comúnmente modelados.

3D­NH /
Sistema Propósito 1D 2D 3D 3D­NHMP Usar
Río Extensiones de las inundaciones
(SW2, DW) P, E, D
Río Velocidades (SW2) P, E, D
Río Comportamiento General Transporte* (SW2) P, E, D
Río Comportamiento detallado del transporte* (SW2) P, E, D
Río
Extensiones de las inundaciones
(SW2, DW) P, E, D
(con curvas)
Río (SW2 con 2D (P,E)
Velocidades/Transporte*
(con curvas) correcciones) 3D (P, E, D)
Flujos supercríticos Extensiones de las inundaciones
(SW2) P, E, D
2D (P,E)
Flujos supercríticos Inundación (SW2 con 3D (P,E)
(con curvas)** Extensiones/Velocidades/Transporte* correcciones) 3D­NH/NHMP
(P,E,D)

Rotura/incumplimiento de presa Extensiones de las inundaciones

(SW2) P, E, D

Rotura/incumplimiento del dique Extensiones de las inundaciones

(SW2, DW) P, E, D

Río
Extensiones de las inundaciones
(SW2) P, E, D
(con estructuras)
Río 2D (P,E)
Velocidades/Transporte* (SW2)
(con estructuras) 3D (P,E)
Embalses/Lagos Tiempo de residencia P, E, D
Embalses/Lagos Transporte* P, E, D
estuarios Agua
(SW2) P, E, D
(Bien mezclada) Superficie/Velocidades/Transporte*
Agua
Estuarios (Estratificados) P, E, D
Superficie/Velocidades/Transporte*
2D/3D (P,E)
Compuertas y aliviaderos** Superficie/velocidades del agua (SW2) 3D­NHMP
(P,E,D)
Nota: DW = Onda difusiva, SW2 = 2D Agua poco profunda, 1D = Unidimensional, 2D = Bidimensional, 3D = 3D
Hidrostático, 3D­NH = 3D­No Hidrostático, 3D­NHMP = 3D­Multifásico No Hidrostático, P = Planificación,
E = Evaluación, D = Diseño, Transporte* = Incluyendo sedimento, ** = Modelo físico recomendado

9­2
Machine Translated by Google

TD­41 Referencias

Referencias
Sociedad Estadounidense de Ingenieros Civiles. 2000. Modelado hidráulico: conceptos y práctica,
Manuales e informes de la ASCE sobre prácticas de ingeniería No. 97. Instituto de Recursos Hídricos
y Ambientales, ASCE, Reston, VA. https://
[Link]/10.1061/9780784404157

Bernard, RS y Schneider, ML 1992. Modelado numérico de profundidad promediada para canales curvos.
Informe técnico HL­92­9, Estación experimental de vías navegables de ingenieros del ejército de EE.
UU., Laboratorio de hidráulica, Vicksburg,
MS. [Link]

Blumberg, AF, B. Galperin y DJ O'Connor. 1992. Modelado de la estructura vertical de flujos en canales
abiertos. Revista de Ingeniería Hidráulica vol. 118, Número 8, 1119­1134. [Link]
10.1061/(ASCE)0733­9429(1992)118:8(1119)

Brown, GL 2012. Un modelo de sedimento suspendido cuasi­3D que utiliza un conjunto de factores de corrección
aplicados a una ecuación de difusión de advección de profundidad promediada. Tercer Simposio
Internacional del IIDH sobre Flujos Superficiales, Iowa City, IA, 4 al 6 de junio de 2012.

Chow, VT 1959. Hidráulica de canal abierto, McGraw­Hill, Nueva York, NY.

Ellet, C. (1853). Los ríos Mississippi y Ohio: contiene planes para la protección del delta contra inundaciones y
mejora de la navegación del Ohio y otros ríos mediante embalses. Informe al Congreso, Lippincot, Grambo, and
Co., Filadelfia, PA.

Finne, J., B. Donnell, J. Letter y RS Bernard. 1999. Corrección de flujo secundario para cálculos de flujo
promediado en profundidad. Revista de Ingeniería Mecánica. vol. 125, Número 7, 848­863. [Link]
10.1061/(ASCE)0733­9399(1999)125:7(848)

Centro de Ingeniería Hidrológica (HEC). 2016. Sistema de análisis de ríos HEC­RAS, Manual del usuario de
modelado 2D, versión 5.0, CPD­68A. Cuerpo de Ingenieros del Ejército de EE. UU., HEC, Davis, CA.

Centro de Ingeniería Hidrológica (HEC). 2016. Sistema de análisis de ríos HEC­RAS, Guía de aplicaciones, Versión
5.0.4, CPD­70. Cuerpo de Ingenieros del Ejército de EE. UU., HEC, Davis, CA.

Centro de Ingeniería Hidrológica (HEC). 2016. Sistema de análisis de ríos HEC­RAS, Manual de referencia
hidráulica, versión 5.0, CPD­69. Cuerpo de Ingenieros del Ejército de EE. UU., HEC, Davis, CA.

Centro de Ingeniería Hidrológica (HEC). 2016. Sistema de análisis de ríos HEC­RAS, Manual del usuario,
Versión 5.0, CPD­68. Cuerpo de Ingenieros del Ejército de EE. UU., HEC, Davis, CA.

Henderson­Sellers, B. 1982. Una fórmula simple para los coeficientes de difusión vertical en remolinos
en condiciones de estabilidad no neutral. Revista de investigación geofísica: océanos, vol. 87, Número
C8, 5860–5864.
[Link]

1
Machine Translated by Google

Referencias TD­41

Mellor, GL y T. Yamada. 1982. Desarrollo de un modelo de cierre de turbulencia para problemas

de fluidos geofísicos. Reseñas de geofísica y ciencia espacial, vol. 20, núm. 4, 851–875. https://
[Link]/
10.1029/RG020i004p00851

Savant, G. 2015. Hidráulica adaptativa tridimensional para aguas poco profundas (ADH­SW3): cierre
de turbulencia. CR­15­1. Cuerpo de Ingenieros del Ejército de EE. UU., Investigación y desarrollo
de ingenieros ­ Laboratorio de hidráulica y costas, ERDC/CHL, Vicksburg, MS. http://
[Link]/11681/1801

Jeremy A. Sharp, Tate. O. McAlpin, Ronald E. Heath, Gary C. Lynch y Howard E. Park.
2013. Investigación hidrodinámica 2D de Olmsted Cofferdams. TR­13­6. Cuerpo de Ingenieros del
Ejército de EE. UU., Investigación y desarrollo de ingenieros ­ Laboratorio de hidráulica y costas,
ERDC/CHL, Vicksburg, MS.
[Link]