CURSO DE INICIACIÓN IBER

UNIDAD 1 PRIMEROS PASOS

pág. 1
CURSO DE INICIACIÓN IBER

1.1. INSTALACIÓN DEL PROGRAMA

1.1.1. Iberaula. Familiarizándonos con el entorno
El programa Iber ha de descargarse gratuitamente desde el portal web de Iberaula, siendo
necesario además registrarse como usuario en el Iberaula.

[Link]

Una vez registrado como usuario, es posible la descarga del programa.

Como bien se explica en el apartado de descargas, los miembros del Aula Iber no
proporcionarán asistencia ni soporte técnico a nivel de usuario. Sin embargo, sí existe un
foro IBER donde los distintos participantes plantean cuestiones.

1.1.2. Registrarse como usuario

Para registrarse como usuario se deberá hacer doble click sobre el link “Inicio Sesión”:

Aparecerá un pequeño cuestionario con datos básicos, que hemos de rellenar.

pág. 2
CURSO DE INICIACIÓN IBER

Damos al botón “Guardar” y nuestro registro habrá finalizado.

1.1.3. Descargar el programa. Documentación

Para descargar el programa, deslizamos el ratón y lo situamos en la pestaña “Modelo Iber”,
apareciendo un desplegable con el menú Descarga. Pincharemos sobre él.

Nos situamos entonces en el entorno Modelo > descarga, en el cual tenemos la posibilidad
de descargar no sólo el programa IBER v 2.6, sino también documentación relativa al
programa. Conviene realizar la lectura del texto del apartado descarga, así como del
contrato de licencia.

pág. 3
CURSO DE INICIACIÓN IBER

Tenemos por tanto a nuestra disposición:

• Iber 2.6
• Manual de referencia hidráulico
• Manual de usuario
• Tutoriales

Al descargar el programa aparece la siguiente ventana:

Tendremos que elegir en función del sistema operativo que tengamos en nuestro PC.

pág. 4
CURSO DE INICIACIÓN IBER

Es importante señalar, que como está escrito en el contrato de licencia, ésta sólo es válida
para un puesto de instalación o una instalación en red. Por tanto, ha de instalarse allí
donde vayamos a hacer uso del programa.

1.1.4. Instalación
Pulsamos doble click sobre la aplicación correspondiente, y aparece el asistente de
instalación de Iber.

pág. 5
CURSO DE INICIACIÓN IBER

Pulsamos “Siguiente”, y a continuación seleccionamos el idioma a utilizar:

Aparece entonces el Acuerdo de Licencia que una vez leído hemos de seleccionar “Acepto
el acuerdo” y pulsar “Siguiente”.

A continuación, nos pedirá el programa la carpeta de destino, pudiendo elegir la que viene
por defecto u otra cualquiera. Continuamos una vez elegida la ubicación dando a “Siguiente”.

Continuamos pulsando “Siguiente”,

pág. 6
CURSO DE INICIACIÓN IBER

Y decimos si queremos o no crear un Icono en el escritorio, lo cual resulta muy útil.

Una vez más pulsamos “Siguiente”, y por último, “Instalar”.

Apareciendo la barra verde que se irá completando según se instale el programa.

Cuando el proceso haya finalizado, aparecerá la siguiente ventana:

pág. 7
CURSO DE INICIACIÓN IBER

Pulsando “Finalizar”, el proceso de instalación habrá terminado con éxito.

1.1.5. Arranque del programa

En el escritorio, hacemos doble clic sobre el icono de Iber 2.6

Y el programa arrancará, apareciendo inicialmente la referencia al pre y postprocesador,

Y finalmente aparecerá el en nuestra pantalla el entorno de Iber, listo para trabajar.

pág. 8
CURSO DE INICIACIÓN IBER

1.2. INTRODUCCIÓN. MODELOS HIDRÁULICOS

1.2.1. Modelos hidráulicos

Hidráulica: se denomina hidráulica a la rama de la física y la ingeniería que se encarga del
estudio de las propiedades mecánicas de los fluidos, dependientes de las fuerzas que se
interponen con la masa y el empuje de la misma.

Para estudiar la mecánica de fluidos, y en particular, el agua, se han utilizado los modelos
hidráulicos, los cuales simulan de manera cercana a la realidad el comportamiento del agua.
Son los modelos físicos y modelos numéricos o computacionales.

[Link]. Modelos físicos

Son reproducciones a escala de distintos escenarios, que simulan los distintos eventos que
puedan materializarse en dicho terreno.

[Link]. Modelos numéricos

Realizan el estudio del comportamiento del agua usando soluciones numéricas de

ecuaciones no lineales que describen mediante modelos los procesos físicos que tienen
lugar.

1.2.2. Tipos de modelos numéricos

Entre todos los modelos numéricos, podemos encontrar distintas clasificaciones:

• En función del enfoque físico del problema.

• Del número de dimensiones con las que se trabaja
• De la variación temporal o no del fenómeno
• De método matemático de resolución de las ecuaciones

pág. 9
CURSO DE INICIACIÓN IBER

[Link]. Enfoque físico

En función del enfoque físico del problema podemos distinguir entre modelos Eulerianos y
modelos Lagrangianos.

MODELOS EULERIANOS: EL cálculo se realiza discretizando la masa del fluido en

elementos (por ejemplo, el tramo de cauce que se encuentra entre dos secciones
transversales), en función de las fuerzas interiores y exteriores que se aplican al elemento
se obtiene las variables hidráulicas representativas del elemento en cada momento.

MODELOS LAGRANGIANOS: Se descompone el fluido en partículas y se asocia a cada

partícula unas coordenadas ξi, el movimiento del fluido estará definido si en todo momento
se conoce la posición ocupada por cada partícula material.

[Link]. Número de dimensiones

En función del número de dimensiones que trata el modelo podemos distinguir entre
modelos de orden 0 (régimen uniforme en una sección), una dimensión (normalmente la X),
dos dimensiones (pueden ser vertical XZ, horizontal XY) y tres dimensiones (XYZ).

MODELOS UNIDIMENSIONALES: Los procesos fluviales se pueden simplificar en flujos de

agua unidimensionales en algunos escenarios.

La profundidad y velocidad sólo varían en la dirección longitudinal del canal, cuyo eje se
supone aproximadamente una línea recta, y la velocidad es constante en cualquier punto de
la sección transversal al citado eje. Es este supuesto, el movimiento del agua se denomina
‘fase fluvial’, esto es, concentrado en colectores prismáticos irregulares.

En este caso el movimiento del agua está gobernado por las ecuaciones de Saint-Venant:
continuidad y conservación de la cantidad de movimiento.

Estas formulan, respectivamente, la relación de caudal entrante y saliente en una superficie

cerrada al ritmo que varía su densidad, y la igualdad entre variación de la cantidad de
movimiento y resultante de las fuerzas exteriores que actúan sobre el fluido:

siendo A, área para la sección transversal al flujo; v, velocidad media; x, abcisa de posición
(situada en el conducto de la sección considerada); B, anchura de la lámina de agua para
un calado y; y, calado (espesor de la lámina de agua); t, tiempo; g, aceleración de la
gravedad; z, cota para la base del conducto; St, pendiente de fricción, que es función de la
velocidad, rugosidad y radio hidráulico.

MODELOS BIDIMENSIONALES: Durante los eventos de avenidas es habitual que el flujo

de agua no quede restringido a la línea central unidireccional de un canal principal, sino que
el caudal desborde ocupando las márgenes del canal principal y las llanuras de inundación.

pág. 10
CURSO DE INICIACIÓN IBER

En estas condiciones, simplificar el flujo a una dimensión resulta insuficiente para modelizar
los flujos secundarios perpendiculares a la dirección principal de avance del agua en el
canal.

La simplificación a 2D lleva asociadas una serie de hipótesis que definen el modelo de aguas
poco profundas (shallow water). En estos modelos, el grosor de la capa de fluido se
considera pequeño comparado con la escala longitudinal horizontal típica. Las hipótesis
fundamentales para estos modelos matemáticos son las siguientes:

• Las ondas que se producen en la superficie varían gradualmente, lo cual es

equivalente a decir que la distribución presiones en la vertical es hidrostática,
o que la aceleración en la vertical es pequeña.
• Las pérdidas por fricción de un flujo transitorio no difieren mucho de las
mismas pérdidas en un flujo estacionario.
• La pendiente promedio del fondo del cauce es tan pequeña que la tangente
puede ser aproximada por el ángulo y las medidas en el fondo son
equivalentes a las medidas en el plano horizontal.
•
Las ecuaciones de aguas someras pueden expresarse de la siguiente forma:

Ecuación de la dinámica

Las ecuaciones de aguas someras son no lineales y el sistema de ecuaciones es hiperbólico,

lo que puede dar lugar a que aparezcan discontinuidades en la solución (saltos hidráulicos
o frentes de onda). Por ese motivo se emplean esquemas numéricos de aproximación, entre
los que destaca el método de los volúmenes finitos, dada su estabilidad de cálculo,
convergencia numérica y velocidad de cómputo.

Los modelos bidimensionales están orientados habitualmente a su computación con mallas

de distintos tipos, cuadrangulares, triangulares, regulares o irregulares, siendo la más
extendida en los modelos la computación sobre mallas triangulares.

MODELOS TRIDIMENSIONALES: La modelización 3D completa de flujos turbulentos se

basa en las ecuaciones de Navier-Stokes promediadas según Reynolds (RANS).

Pueden expresarse de la siguiente manera para un fluido no newtoniano incomprensible:

pág. 11
CURSO DE INICIACIÓN IBER

Este tipo de modelos precisa mucho tiempo de cálculo, con lo cual no es de utilidad en la
gran mayoría de las aplicaciones de la dinámica fluvial.

[Link]. Variación temporal del caudal

Desde el punto de vista de la variación temporal los modelos se pueden dividir en:

MODELOS EN RÉGIMEN PERMANENTE: el caudal y las condiciones de contorno

permanecen constantes a lo largo del tiempo.

El flujo uniforme tiene lugar cuando la pendiente del lecho (S0) es igual a la pendiente de
fricción (Sf) y todos los demás términos de la ecuación de conservación de la cantidad de
movimiento se consideran despreciables.

Al ser el flujo es permanente, la línea de energía, la superficie de la lámina de agua y la

solera son líneas paralelas. Considerando esta simplificación metodológica, existen varias
fórmulas empíricas que permiten obtener las relaciones caudales-calados de forma
biunívoca, en función de parámetros tan simples como la rugosidad del lecho.

La más conocida es la ecuación de Manning, que relaciona la velocidad con parámetros

morfométricos y con un índice de rugosidad tabulado para los distintos materiales del lecho:

Fórmula de Manning

siendo: v, velocidad media de la corriente; R, radio hidráulico (A/P); S, pendiente para la

línea de energía (coincidente con la del lecho); n, coeficiente de rugosidad.

El número adimensional n, básico en la ecuación de los caudales (Q=v·A) o de calados a

partir de estos, se encuentra tabulado para las diferentes tipologías de lechos fluviales,
según su rugosidad y la presencia de obstáculos al lecho.

MODELOS EN RÉGIMEN GRADUALMENTE VARIADO: El caudal es variable en el

espacio.

El flujo es permanente en el tiempo, con un caudal constante en el tiempo. Sin embargo,

existe una variación de la velocidad en el espacio, y por tanto del calado al no modificarse
el caudal. En este caso el régimen recibe el nombre de gradualmente variado, y en él se
produce una distribución hidrostática de las presiones. Los perfiles pueden analizarse
considerando régimen supercrítico y subcrítico.

Cuando el régimen pasa de lento a rápido y viceversa, la profundidad debe pasar por valores
críticos (Yc), conociéndose en el primer caso estos puntos como secciones de control, y

pág. 12
CURSO DE INICIACIÓN IBER

produciéndose un característico resalto hidráulico en el segundo de los casos. En tales

situaciones puede estimarse el número de Froude a partir de la ecuación de Froude:

Fórmula de Froude

siendo: v, velocidad del flujo; g, aceleración de la gravedad; hm calado de la sección.

MODELOS EN RÉGIMEN VARIABLE: tanto el caudal como las condiciones de contorno

pueden variar a lo largo del tiempo.

Cuando se trata de estimar calados modelizando un flujo no permanente, esto es, con
variaciones de la velocidad de la corriente en el tiempo, este tipo de modelos persiguen la
resolución de las ecuaciones de Saint-Venant empleando distintos esquemas numéricos.

Fórmula de Saint-Venant

Siendo U el vector de las variables del flujo, F es el tensor del flujo y H es término
independiente o término fuente.

[Link]. Método matemático de resolución

Atendiendo al método matemático de resolución de las ecuaciones podemos encontrar las

siguientes metodologías:

Al ser el agua un medio continúo constituido por un número casi infinito de partículas.
Determinar la velocidad de cada partícula es virtualmente imposible y poco práctico; por ello
este continuo se divide en elementos discretos de tamaño finito, que la memoria de una
computadora es capaz de manejar. Estos métodos resuelven las ecuaciones que tienen
dominio continuo mediante la solución en un número finitos de puntos discretos en dicho
dominio. Cuando los valores buscados (tirante, velocidad, etc.) en dichos puntos discretos
son encontrados, la solución en cualquier otro punto puede ser aproximada mediante
métodos de interpolación. Estos modelos numéricos proporcionan mucho detalle y precisión,
siendo capaces de manejar condiciones de borde e iniciales complejas, para los cuales no
existen en la mayoría de casos soluciones analíticas. Sin embargo, la mayoría de estos
métodos requieren la solución simultánea de todos los nodos del dominio, lo que incrementa
la carga computacional haciéndolos más lentos.

MODELOS EN DIFERENCIAS FINITAS (MDF): discretización del fluido en diferencias

finitas.

MODELOS EN ELEMENTOS FINITOS (MEF): Se discretiza el fluido en elementos finitos.

MODELOS EN VOLÚMENES FINITOS (MVF): Por último y más complejo, el método de

volúmenes finitos.

pág. 13
CURSO DE INICIACIÓN IBER

MODELOS LAGRANGIANOS (SPH): Entre los métodos lagrangianos el más destacado es

el método de Hidrodinámica de Partículas Suavizadas (Smoothed Particle Hydrodynamics,
SPH).

Los modelos SPH se caracterizan, en primer lugar, por la discretización Lagrangiana del
medio continuo, mediante un conjunto de partículas ó parcelas de fluido y en segundo lugar,
a la reducción del problema al estudio de la interacción entre las mismas, en base a un
esquema de promediado que permite pasar de los valores puntuales de las variables de
interés (presión, densidad, velocidad, posición) a los valores correspondientes de campo
continuo. En cierto modo esta discretización permite asimilar el problema de la mecánica de
fluidos al clásico de n cuerpos de la mecánica celeste, lo que hace que no sea casual el
desarrollo inicial de esta técnica en la comunidad astrofísica.

1.3. MODELO IBER

1.3.1. Modelo Iber. Aplicación y capacidades

La interfaz de Iber está basada en GID ([Link]), el software de Pre-proceso y

Postproceso desarrollado por el CIMNE (Centro Internacional de Métodos Numéricos en
Ingeniería. El programa surge para poder disponer de una herramienta con la cual realizar
simulaciones hidráulicas que faciliten la aplicación de la legislación. Ante la existencia de los
programas Turbillón y Carpa, Iber surge como unión de ambos.

Los campos de aplicación de Iber son:

• Simulación del flujo en lámina libre en cauces naturales

• Evaluación de zonas inundables. Cálculo de las zonas de flujo preferente
• Cálculo hidráulico de encauzamientos
• Cálculo hidráulico de redes de canales en lámina libre
• Cálculo de corrientes de marea en estuarios
• Estabilidad de los sedimentos del lecho
• Procesos de erosión y sedimentación por transporte de material granular

Las capacidades y características más destacadas del modelo Iber son las siguientes:

• Resolución integrada de las ecuaciones de Saint Venant 2D

• Esquemas explícitos en volúmenes finitos con mallas no estructuradas
• Capacidad de resolver flujo subcrítico y supercrítico, incluyendo resaltos
hidráulicos móviles
• Mojado y secado del dominio con la conservación exacta del volumen de
agua
• Modelización de la turbulencia mediante modelos de diferente complejidad
• Cálculo de la infiltración
• Tensión superficial por viento
• Estructuras internas: puentes, compuertas y vertederos
• Delimitación de la zona de flujo preferente según RDPH (Reglamento de
Dominio Público Hidráulico). Delimitación de la vía de intenso desagüe y
zonas de grave riesgo para personas y bienes.
• Evolución del lecho debido a transporte de sedimentos por carga de fondo y
en suspensión
• Interfaz amigable de pre y postproceso
• Integración en GIS
• Verificado y contrastado con soluciones analíticas, con otros modelos, con
ensayos de laboratorio y con medidas de campo

El modelo Iber es un modelo que se encuentra en continuo desarrollo, y está abierto a

mejoras o adaptaciones. Entre las líneas prioritarias de desarrollo del modelo se encuentran:

pág. 14
CURSO DE INICIACIÓN IBER

• Evaluación del hábitat fluvial

• Definición de caudales ecológicos con métodos hidrobiológicos
• Nuevas condiciones de contorno internas
• Transporte de mezclas de sedimento
• Modelos de calidad de aguas

El modelo Iber, para poder realizar todas estas funciones, consta de tres módulos distintos:
se incluye un módulo hidrodinámico, un módulo de turbulencia, y un módulo de transporte
de sedimentos por carga de fondo y por carga en suspensión. Tanto las ecuaciones
hidrodinámicas, como los modelos de turbulencia y el modelo de transporte de sedimentos,
se resuelven en forma integral por el método de volúmenes finitos en una malla no-
estructurada. El método de volúmenes finitos es uno de los más extendidos y comúnmente
utilizados en dinámica de fluidos por ser muy eficiente para la resolución de leyes de
conservación.

1.3.2. Módulo hidrodinámico

El módulo hidrodinámico resuelve las ecuaciones de aguas someras promediadas en

profundidad bidimensionales, también conocidas como ecuaciones de [Link] 2D.
Dichas ecuaciones asumen las siguientes hipótesis:

• Distribución de presión hidrostática

• Distribución uniforme de velocidad en profundidad
La hipótesis de presión hidrostática se cumple razonablemente en el flujo en ríos, así como
en las corrientes generadas por la marea en estuarios y zonas costeras.

La hipótesis de distribución uniforme de velocidad en profundidad se cumple de forma

habitual en ríos y estuarios, siempre y cuando no existan procesos relevantes de
estratificación debido a diferencias de salinidad, de temperatura o al viento. En la actualidad,
los modelos numéricos basados en las ecuaciones de aguas someras promediadas en
profundidad bidimensionales son los más utilizados en estudios de evaluación de zonas
inundables, dinámica fluvial, corrientes de marea, transporte de sedimentos y calidad de
aguas. El hecho esencial en los modelos llamados de aguas poco profundas, es que el
grosor h de la capa de fluido es pequeño no comparado con la escala longitudinal horizontal
típica L.

Para entender el alcance de las hipótesis asociadas, hay que deducir las ecuaciones de
aguas poco profundas a partir de la ecuación de conservación de la masa, lo cual se
desarrollará posteriormente en el curso.

El módulo hidrodinámico de Iber tiene la capacidad de considerar los siguientes procesos:

• Flujo no estacionario en régimen rápido y en régimen lento

• Formación de resaltos hidráulicos no estacionarios
• Fricción de fondo según formulación de Manning
• Tensiones turbulentas calculadas según diversos modelos de turbulencia
• Viscosidad turbulenta constante
• Modelo parabólico
• Modelo de longitud de mezcla
• Modelo k-e
• Variación temporal de la cota del fondo debido a transporte de sedimentos
• Condiciones de contorno abierto tipo:
• Hidrograma (entrada / salida)
• Marea
• Vertido crítico
• Vertedero
• Curva de gasto

pág. 15
CURSO DE INICIACIÓN IBER

• Condiciones de contorno tipo pared

• Deslizamiento libre
• Fricción de pared según ley logarítmica
• Condiciones de contorno internas:
• Puentes (coeficiente de pérdida de carga / coeficiente de arrastre)
• Vertedero
• Compuerta
• Muro / Dique encauzamiento
• Infiltración según las formulaciones de:
• Green-Ampt
• Horton
• Lineal
• Rozamiento superficial por viento según formulación de Van Dorn
• Salida de resultados de Riesgo según RDPH
• Utilidades para el cálculo de la zona de flujo preferente según RDPH
• Frentes de inundación no estacionarios

1.3.3. Módulo de turbulencia

En gran parte de procesos hidráulicos, la turbulencia tiene un papel importante y es

necesario disponer de una herramienta con la que simular estos procesos de viscosidad
turbulenta, energía turbulenta y disipación.

El modelo Iber incorpora diversos modelos de turbulencia tipo Boussinesq para el cálculo de
las tensiones tangenciales turbulentas.

Se resuelven en el módulo de turbulencia e intervienen tanto en el cálculo de la

hidrodinámica como en el del transporte en suspensión. Se incluyen los siguientes modelos
de turbulencia tipo Boussinesq para aguas someras:

• Viscosidad turbulenta constante

• Modelo parabólico
• Modelo de longitud de mezcla
• Modelo k-ε de Rastogi y Rodi
Esto permite seleccionar el más adecuado en cada caso, teniendo en cuenta la complejidad
del modelo, tiempo de cálculo, complejidad del fujo y necesidades de precisión del cálculo.

1.3.4. Módulo de transporte de sedimentos

En el módulo de transporte de sedimentos se resuelven las ecuaciones de transporte por

carga de fondo y por carga en suspensión. Teniendo en cuenta ambos modos de transporte
se calcula la evolución de la cota del fondo debido a procesos de sedimentación y erosión
mediante la ecuación de Exner. En la versión actual únicamente se consideran
granulometrías uniformes.

[Link]. Transporte de sedimento por carga de fondo

El caudal sólido de fondo se calcula mediante formulaciones empíricas en función de la

tensión de fondo. El módulo de transporte de sedimentos por carga de fondo incluye las
siguientes características:

• Umbral de movimiento de Shields

• Formulaciones para caudal sólido de fondo
• Meyer Peter-Muller con corrección de Wong-Parker (D=2-
30mm)
• Van Rijn (D=0.2-2mm)

pág. 16
CURSO DE INICIACIÓN IBER

• Corrección por pendiente de fondo en inicio del arrastre

(tensión crítica en talud)
• Corrección por pendiente de fondo en transporte sólido
(magnitud y dirección)
• Marcelo García / Mhin-Duc
• Separación de tensiones de Einstein por
formas de fondo y grano
• Condiciones de contorno tipo
sedimentograma (caudal sólido de fondo
variable en tiempo)
• Condición de cota de fondo no erosionable
(puntos fijos)
[Link]. Transporte de sedimento en suspensión

El transporte en suspensión se calcula resolviendo la ecuación de convección-difusión

promediada en profundidad para la concentración de sedimento, incluyendo un término de
sedimentación / erosión que modela el intercambio de sedimento entre el lecho y la carga
en suspensión.

Las principales características de este módulo son:

• Incorporación de transporte por difusión turbulenta

• Término de deposición / Re suspensión
• Cálculo de la concentración de sedimento en suspensión según
formulaciones de:
• Van Rijn
• Smith
• García
• Cálculo de la velocidad de sedimentación de las partículas
según formulación de van Rijn
• Condición de contorno de concentración SS variable en
tiempo
1.3.5. El entorno de Iber

El entorno Iber consta de una interfaz amigable que resulta sencilla de manejar cuando se
conocen las capacidades del programa.

Incluye módulos de preproceso de datos y postproceso de resultados basados en el sistema

GiD, y ambos módulos son compatibles con entornos SIG, con la facilidad de importado y
exportado de datos que ello conlleva.

En el preproceso se introducen los datos, manejando dos conceptos; el de la introducción

de geometría y mallado con las condiciones de contorno oportunas.

[Link]. Preproceso

En el preproceso se introducen los datos, manejando dos conceptos; el de la introducción

de geometría y mallado con las condiciones de contorno oportunas.

La introducción de la geometría consiste en introducir en el modelo una representación más

o menos idealizada del problema a estudiar. Se crea y se edita el modelo y se le asignan las
condiciones de contorno y propiedades. Una vez realizado esto, ha de mallarse la superficie.
Se realiza por tanto una discretización del modelo generado a partir de la geometría. Se crea
y edita una malla que recoge las características de la geometría. Sobre esta malla se realiza
el cálculo.

pág. 17
CURSO DE INICIACIÓN IBER

A la hora de realizar estas operaciones, hemos de tener en cuenta que el preprocesador

tiene las siguientes capacidades:

• Aplicaciones CAD para diseño de geometrías de cálculo

• Generación de mallas estructuradas y no-estructuradas mediante diferentes
algoritmos de mallado
• Capacidad de importar imágenes georreferenciadas
• Definición de propiedades y parámetros de cálculo variables espacialmente
• Posibilidad de importar la rugosidad de fondo de forma automática a partir
de datos geo-referenciados de usos del suelo
• Mensajes de aviso frente a errores comunes
• Posibilidad de reemprender un cálculo interrumpido

[Link]. Postproceso

El postproceso es el proceso posterior al cálculo numérico, y que nos permite volcar los
resultados que necesitemos de manera gráfica y muy visual, además de permitirnos obtener
resultados numéricos en puntos concretos para poder analizar que está ocurriendo
realmente.

El postprocesador de Iber permite realizar todas las operaciones siguientes:

• Visualización de campos 2D escalares y vectoriales

• Extracción de secciones y perfiles longitudinales
• Visualización de resultados temporales mientras se ejecuta la simulación
• Visualización de zonas de riesgo según Reglamento del Dominio Público
Hidráulico
• Visualización de zonas inundables y frentes seco-mojado no estacionarios
• Herramientas para la visualización de los resultados en formato SIG
• Exportación de resultados en formato raster para ser visualizados
directamente en entornos SIG
• Visualización de resultados sobre imágenes geo-referenciadas
• Creación de videos con animación de resultados no estacionarios
•
[Link]. Mallador de Iber

Uno de los procesos que requieren mayor tiempo y esfuerzo a la hora de desarrollar un
estudio de simulación numérica del flujo en ríos es la generación de la malla de cálculo. Un
río tiene una geometría irregular y la construcción de una malla eficiente no es evidente. Una
buena malla debe de ser irregular, con el fin de minimizar el número de elementos con
transiciones suaves. Para ello son muy adecuados los métodos de mallado basados en el
error cordal (máxima diferencia entre el terreno original y la malla).

En Iber se incluyen los siguientes métodos de generación de mallas:

• Mallas estructuradas
• Malla de triángulos irregulares generada con error cordal
• Malla de triángulos rectángulos equiláteros o RTIN (Rectangular
Triangulated Irregular Network), generada con error cordal
• Malla regular con todos los puntos del MDT (Modelo Digital del Terreno)
• Importación de mallas generadas con SMS (formato .2dm)
Para estas metodologías se requiere como información de partida un Modelo Digital del
Terreno en formato ASCII de arc/Info.

pág. 18
CURSO DE INICIACIÓN IBER

1.4. BASES CONCEPTUALES HIDRÁULICAS

1.4.1. Ecuaciones gobernantes
[Link]. Ecuaciones de partida

El módulo Iber, como se ha comentado, resuelve las ecuaciones de aguas someras

promediadas en profundidad bidimensionales, también conocidas como ecuaciones de
[Link] 2D.

Si bien el flujo del agua es 3D, se considera la simplificación a 2D. Esta aproximación lleva
asociadas una serie de hipótesis que definen el modelo de aguas poco profundas (shallow
water). En el modelo de flujo llamado de aguas poco profundas el hecho esencial es que el
grosor h de la capa de fluido es pequeño no comparado con la escala longitudinal horizontal
típica L.

Para entender el alcance de las hipótesis asociadas, hay que deducir las ecuaciones de
aguas poco profundas a partir de la ecuación de conservación de la masa.

(1.4.1)

donde u, v, w son las tres componentes de la velocidad del flujo en las direcciones x, y, z
respectivamente, y de las ecuaciones de Navier-Stokes para flujo incompresible, que en
forma diferencial conservativa se escriben como:

(1.4.2)

con ρ la densidad del fluido, g el valor de la aceleración en la vertical, τij las componentes
del tensor de esfuerzos y p la presión.

Para fijar las soluciones de las ecuaciones diferenciales (1.4.1) y (1.4.2) es preciso definir
las condiciones de contorno en las fronteras. En este caso existen dos zonas: la entre fase
fluido-sólido (fondo) que es fija y la entre fase fluido- fluido (superficie libre) que puede

pág. 19
CURSO DE INICIACIÓN IBER

cambiar continuamente. Las condiciones de frontera en ambas superficies se presentan en

dos formas: cinemáticas y dinámicas.

Las condiciones cinemáticas están relacionadas con la velocidad y nos dicen que las
partículas de agua en su movimiento no pueden cruzar ninguna frontera. Para el fondo
significa, que la componente de la velocidad normal a la superficie sólida debe ser cero
(fondo sólido, impermeable y fijo).

Las condiciones de contorno dinámicas nos dan información sobre las fuerzas que actúan
en los contornos. Si el flujo es viscoso y el fondo fijo, la fuerza que actúa es la de la
viscosidad y por lo tanto las partículas que se encuentran en contacto con el fondo están
pegadas a él, por lo cual puede imponerse la condición de no deslizamiento. En la superficie
libre se supone continuidad de esfuerzos.

[Link]. Ecuaciones de aguas poco profundas 3D

Para especificar un modelo simplificado de aguas poco profundas se tienen que considerar
escalas típicas que permitan realizar un análisis de órdenes de magnitud.
Las escalas longitudinales en la vertical incluyen:

• la profundidad de agua h;
• el grosor δ de la capa límite en el fondo o en la superficie libre. En ríos, lagos
y costas, las capas límite suelen ser más gruesas que la profundidad de
agua, y entonces la capa de agua suele encontrarse incluida dentro de la
capa límite. En este caso, la profundidad h es la escala vertical relevante.
• la variación del nivel de agua (es decir, amplitud de ondas generadas en el
agua), que vuelve a ser una escala interna si no se impone directamente
desde el exterior.
Las escalas longitudinales en la horizontal también pueden ser:

• las dimensiones físicas del receptáculo que se considera: anchura de un río,

longitud de un puerto, embalse;
• distancias en las cuales las fuerzas externas varían de forma significativa:
variaciones del viento y de la presión atmosférica en una tormenta;
• la longitud de onda, que es una escala interna determinada por otros factores
como la frecuencia de las mareas y las dimensiones físicas.
Con las hipótesis de ciertas escalas también se pueden simplificar las ecuaciones y
tendríamos por lo tanto las ecuaciones de aguas poco profundas 3D:

(1.4.3)

(1.4.4)

(1.4.5)

(1.4.6)

pág. 20
CURSO DE INICIACIÓN IBER

[Link]. Ecuaciones de aguas poco profundas 2D (promediadas en profundidad)

Una vez obtenidas las ecuaciones de aguas poco profundas en 3D, pasamos a la forma
bidimensional las ecuaciones. Las hipótesis fundamentales para estos modelos
matemáticos son las siguientes:

• Las ondas que se producen en la superficie varían gradualmente, es decir,

la distribución presiones en la vertical es hidrostática, o que la aceleración
en la vertical es pequeña.
• Las pérdidas por fricción de un flujo transitorio no difieren mucho de las
mismas pérdidas en un flujo estacionario.
• La pendiente promedio del fondo del cauce es tan pequeña que la tangente
puede ser aproximada por el ángulo y las medidas en el fondo son
equivalentes a las medidas en el plano horizontal.
Para pasar a la forma bidimensional de las ecuaciones, dejando la profundidad como
variable dependiente, hay que realizar más simplificaciones. Para eliminar de las ecuaciones
el movimiento en la dirección vertical z se promedian las ecuaciones en esta dirección
haciendo uso de las definiciones de los promedios de las variables.

(1.4.7)

(1.4.8)

Las ecuaciones estándar de aguas poco profundas, se obtienen haciendo algunas

aproximaciones más como despreciar los esfuerzos laterales, el efecto de Coriolis las
fuerzas externas (viento, gradientes de presión atmosféricos).

Los términos de fricción del agua con el fondo del cauce se representan por Sf, pendiente
de la línea de energía en cada dirección

(1.4.9)

(1.4.10)

dando lugar al siguiente sistema de ecuaciones que es la forma más conocida de

representación del sistema de ecuaciones de aguas poco profundas

pág. 21
CURSO DE INICIACIÓN IBER

(1.4.11)

(1.4.12)

(1.4.13)

[Link]. Forma vectorial de las ecuaciones de aguas poco profundas 2D (promediadas

en profundidad)

Este sistema de ecuaciones en su forma conservativa, es decir, escritas las ecuaciones de

la forma más cercana posible a un sistema de leyes de conservación de masa y cantidad de
movimiento, es

(1.4.14)

(1.4.15)

U representa el vector de variables conservadas, h profundidad del agua, hu y hv caudales

unitarios a lo largo de las direcciones coordenadas x, y respectivamente, F y G son los flujos
de las variables conservadas a través de los lados de un volumen de control, y contienen el
flujo convectivo y los gradientes de presión hidrostática. La parte derecha de la igualdad en
el sistema de ecuaciones, H, contiene las fuentes y sumideros de la cantidad de movimiento
a lo largo de las dos direcciones coordenadas, provenientes de las variaciones del fondo del
cauce y de las pérdidas por fricción que deben estar relacionadas con el campo de
velocidades.

Son ecuaciones no lineales y el sistema de ecuaciones es hiperbólico, lo que puede dar

lugar a que aparezcan discontinuidades en la solución (saltos hidráulicos o frentes de onda).
Por ese motivo se emplean esquemas numéricos de aproximación, entre los que destaca
el método de los volúmenes finitos, dada su estabilidad de cálculo, convergencia
numérica y velocidad de cómputo.

1.4.2. Propiedades de los sistemas de ecuaciones hiperbólicas

[Link]. Método de las características

pág. 22
CURSO DE INICIACIÓN IBER

El método de las características, puede ser descrito para un problema 1D como una técnica
en la que el problema de resolver dos ecuaciones diferenciales en derivadas parciales puede
ser sustituido por el de resolver cuatro ecuaciones diferenciales ordinarias.

Esta descripción implica una situación continua en la que las derivadas están definidas en
todas partes y en casi todas las direcciones.

Para comprender su utilidad, que reside principalmente en darnos una visión intuitiva del
comportamiento de las soluciones de las ecuaciones diferenciales hiperbólicas, nos
situaremos en el caso más simple y seguiremos la deducción realizada por Abbott [Abbott
1992]. Consideremos un canal prismático horizontal, de sección rectangular y anchura
unidad, y despreciemos todos los términos de pérdida de energía. Tal y como hemos visto,
las ecuaciones de conservación de masa y cantidad de movimiento pueden ser escritas:

(1.4.16)

o equivalentemente

(1.4.17)

Utilizando que c2 = g h, y mediante algunas transformaciones simples del sistema escrito

bajo la última forma, obtendremos:

(1.4.18)

pág. 23
CURSO DE INICIACIÓN IBER

Este resultado es análogo a (1.3.17), obtenido de realizar un cambio de variables en el

sistema original. Si las comparamos con la expresión general para la variación de una
cantidad cualquiera s = s (x, t)

(1.4.19)

identificaremos su significado. Podemos interpretar las ecuaciones (1.3.18) como que la

cantidad u + 2c se mantiene constante a lo largo de una línea en el plano (x, t) dada por la
pendiente u + c. Y, análogamente para la cantidad u − 2c.

(1.4.20)

Las direcciones definidas en el plano por estas ecuaciones se denominan direcciones

características (positiva, C +, y negativa, C −, respectivamente). A las cantidades
conservadas a lo largo de ellas se les llama invariantes de Riemann (J + y J −,
respectivamente). Estas son las cuatro ecuaciones diferenciales ordinarias que reemplazan
a las ecuaciones diferenciales en derivadas parciales en el método de las características.

Podemos decir que en un problema descrito por unas ecuaciones como las de aguas poco
profundas, existen dos direcciones características en cada punto del plano (x, t) y, cuando
se dibujan las tangentes, aparecen dos familias de líneas características a lo largo de las
cuales las cantidades correspondientes (u + 2c) o (u − 2c) son invariantes.

Las características Proporcionan información importante del comportamiento de la solución.

Los siguientes gráficos son muy representativos:

pág. 24
CURSO DE INICIACIÓN IBER

1.5. ESQUEMAS NUMÉRICOS. VOLÚMENES

FINITOS
1.5.1. Resolución por volúmenes finitoS
En un esquema de resolución en volúmenes finitos se discretiza el fluido en volúmenes
finitos de modo que:

• El flujo pasa a través de las aristas de las celdas

• Lo que sale de una celda entra en la celda de al lado
• Se realiza un balance de entrada y salida para cualquier variable

pág. 25
CURSO DE INICIACIÓN IBER

Por ello tiene las siguientes ventajas respecto a los modelos de elementos finitos y
diferencias finitas:

• Impone conservación de forma natural

• Existe una gran flexibilidad geométrica
• Se resuelven las ecuaciones en forma integral
• La discretización es muy intuitiva
Los desarrollos de los esquemas de resolución de las ecuaciones son bastantes complejos
y a continuación daremos unas pinceladas simplemente por recordar o conocer las
diferentes posibilidades.

Como hemos comentado anteriormente, se resuelven las ecuaciones en su forma

conservativa, y por tanto:

Por tanto, queda a resolver la siguiente ecuación con un término de convección y otro de
difusión:

1.5.2. Esquemas centrados

Para resolver el término convectivo, se puede discretizar el flujo convectivo mediante un
esquema centrado como el que sigue, sin considerarse la dirección en que se transmite la
información:

pág. 26
CURSO DE INICIACIÓN IBER

1.5.3. Esquemas descentrados

Otra opción es discretizar el flujo convectivo mediante un esquema centrado considerando
la dirección en que se transmite la información, pueden ser como los centrados, de orden 1
o de orden 2.

[Link]. Esquemas descentrados de primer orden

El siguiente esquema descentrado es de orden 1.

El término Gn de difusión numérica estabiliza la solución, aunque introduce errores.

Este esquema de resolución tiene las características siguientes:

• Es numéricamente estable
• No genera oscilaciones de Φ
• Muy difusivo
•
[Link]. Esquemas descentrados de segundo orden

Entre ellos están los esquemas TVD (Total Variation Diminishing) de alta resolución

• Reducen el valor de la difusión numérica manteniendo la solución estable

• Son más precisos, pero menos estables que esquemas orden 1

Estos esquemas se resuelven siguiendo los siguientes pasos:

• Se realiza una extrapolación lineal de las variables en cada celda a partir del
valor medio en la celda y del gradiente.
• Posteriormente se extrapolan los valores de las variables de nodos a aristas

pág. 27
CURSO DE INICIACIÓN IBER

• Estos valores extrapolados se utilizan en vez de los valores nodales en el

esquema numérico
El modelo Iber utiliza un esquema conservativo descentrado de alta resolución.

1.5.4. Discretización temporal

Para la resolución de estos esquemas numéricos se realiza una discretización temporal,

Cumpliendo la condición de CFL (Courant-Friedrichs-Levy)

Con lo que se obtiene la condición de un paso de tiempo de cálculo de valor:

Esta condición será necesario cumplirla para conseguir estabilidad.

pág. 28