PRE TRATAMIENTO DE DATOS

El modelo bidimensional restringe el ingreso de información procesada a ciertos formatos,

entre los formatos que se pueden leer en el modelo se encuentra el sistema de
codificación ASCII. Así, a través de un SIG se transforman los datos a este tipo de código.
Se ocupa una imagen en formato JPG de la zona de interés creada a partir de un SIG, la
cual debe estar referenciada, con la finalidad de elaborar con mayor detalle la zona de
modelación.
METODOLOGÍA IBER
La simulación bidimensional toma en cuenta la siguiente metodología (Figura 1), la cual
se divide en tres etapas; Pre-proceso, proceso y pos-proceso. Cada etapa contiene
diferentes rubros para realizar la simulación de las manchas de inundación asociadas a
los diferentes periodos de retorno de análisis.

1.- PRE-PROCESO: Entrada de Datos

- Geometría
- Malla
- Condiciones hidrodinámicas (Entrada, Salida e Iniciales)
- Rugosidad.
- Condiciones hidrológicas (Definición y asignación de hietograma y pérdidas)
- Condiciones de Sedimentos (Capa no erosionable, transporte de fondo, entrada en
suspensión e inicial de suspensión)

2.- PROCESO: Simulación hidráulica

- Ejecución de cálculo
- Parámetros de cálculo (tiempo de simulacipon, parámetros de cálculo, selección de
resultados, elección de modelo de turbulencia, selección de modelo de transporte de
sedimentos)

3 .- POS-PROCESO: Resultados
- Visualización de resultados
- Creación de gráficos
- Exportar resultados en formato raster (tirante, velocidad, erosión y sedimentación, número
de froude, caudal especifico)
PRE-PROCESO

Geometría
La geometría establece el dominio físico de análisis del modelo numérico, y también a
partir de ella se asignan algunas de las condiciones, así como el refinamiento del mallado;
por ello, es importante identificar la estructuración más conveniente con la finalidad de
facilitar tales asignaciones posteriores.
Una geometría puede ser conformada por elementos como líneas, puntos y superficies, a
los cuales se le asignan las características deseadas en conjunto o de manera individual.
Las características de la geometría dependen del grado de precisión de los análisis y
resultados que se quieran obtener.
Mallado
La malla determina el dominio de cálculo (computacional) al discretizar el dominio físico
de interés, es decir, el cálculo se realiza para cada elemento que conforma el mallado, por
ello es tan importante optimizar esta condición pues va directamente relacionada al
tiempo de cálculo, que conforme aumenta el número de elementos también el tiempo de
procesamiento.
De acuerdo a la estructuración de la geometría, se trabaja con un mallado no estructurado
e irregular, integrado por elementos triangulares. Y, el tamaño de los elementos se
determinan tomando en cuenta: las características de las zonas de interés, las inundables
y las críticas.
Para poder determinar los tamaños óptimos, se debe llevar a cabo una calibración entre
tiempo de simulación o costo computacional, y tamaño mínimo requerido para los
elementos.
Condiciones hidrodinámicas
Para ejecutar el módulo hidrodinámico se deben introducir las condiciones de contorno
hidrodinámicas (entrada y salida) y las condiciones iniciales.

 Condiciones de contorno de entrada

Está condición varía dependiendo qué se quiera simular; las condiciones de contorno de
entrada para el caso de un análisis hidrodinámico se realizan por medio de la asignación
de los hidrogramas de diseño para los diferentes periodos de retorno que se quieran
analizar.

 Condiciones de contorno de salida

Las condiciones de salida se realizan por medio de una condición de flujo (subcritico,
crítico o supercrítico), el programa se ajusta según topografía y caudal de entrada.
Además se puede proporcionar una condición de salida proporcionando un nivel del agua
conocido o una curva de variación de niveles con respecto del tiempo, para casos donde
se simule en la costa.
  Condiciones iniciales

La cantidad de agua existente en la corriente en estudio en el momento de iniciar el

cálculo, este dato se puede introducir como una elevación o como un tirante de agua,
para algunos casos prácticos se estipula que este valor puede ser “0”, para simular que la
corriente no contiene agua al inicio de la simulación.
Rugosidad
La rugosidad se establece a través de coeficiente de Manning, el programa cuenta con
valores predeterminados, los cuales se pueden definir según usos de suelo para
caracterizar las superficies que conformen la geometría o, en su caso, a los elementos de
la malla.
Procesos hidrológicos
Para realizar los procesos hidrológicos en IBER, se debe definir el hietograma de estudio,
dependiendo de los periodos de retorno que se quieran analizar, se asignan dichos
hietogramas a las superficies o mallas creadas y finalmente se asigna un modelo de
perdida por infiltración, el cual depende de los datos que se tengan.

 Definición de hietograma

La lluvia se asigna en forma de hietograma. Cada instante de tiempo del hietograma

indica el inicio de un bloque de precipitación de intensidad constante que se mantiene
hasta el próximo instante de tiempo.

 Asignación de hietograma

Se le asigna el hietograma elaborado en el rubro anterior a cada una de las geometrías

creadas, se puede tener más de un hietograma.

 Pérdidas

La versión actual de Iber permite considerar tres modelos distintos de pérdidas por
infiltración: el modelo de Horton, el modelo lineal, y el modelo de Green-Ampt. Con el
menú se accede a la introducción de los distintos parámetros de cada modelo, y a su
asignación a superficies o elementos de la malla.
Transporte de sedimentos
El módulo de sedimentos al igual que los de hidrodinámica y procesos hidrológicos,
requieren la asignación de condiciones iniciales y de frontera.
A través del menú se accede a los submenús que permiten asignar condiciones a la
geometría o a la malla para el cálculo del transporte de sedimentos. La activación de los
distintos módulos de cálculo y la declaración de los parámetros generales se debe hacer
desde el menú “Datos del problema”.

 Capa no erosionable
Es la posición capa de roca se puede fijar una cota o profundidad del lecho no
erosionable. Sirve tanto para transporte de fondo como para transporte en suspensión.

 Sedimentograma para transporte de fondo

Sirve para asignar la cantidad de transporte de fondo a la entrada de la zona de interés,

se accede con el submenú Cond Cont Transp Sed Fondo a la cual se asigna a la condición de
contorno aguas arriba para caudal sólido por arrastre de fondo. La condición se asigna a
las líneas del contorno o frontera de análisis, o lados de elementos del contorno. Se
puede escoger entre:

 Agua limpia: caudal sólido nulo en la entrada.

 Capacidad de arrastre: el caudal sólido de entrada corresponde a la
capacidad de arrastre del flujo en dicha entrada, que a su vez
depende del modelo de transporte de fondo escogido en “Datos de
problema”
 Dependiente del tiempo: permite asignar un sedimentograma a la
entrada.

 Entrada de sedimento en suspensión

Sirve para asignar la cantidad de transporte en suspensión a la entrada de la zona de

interés, se accede en el submenú Entrada Sed Suspensión asigna una concentración de
sedimento en suspensión (puede ser variable en el tiempo) a una entrada de caudal.

 Condición inicial de sedimento en suspensión

El submenú Cond Inicial Sed Suspensión permite asignar una concentración inicial de
sedimento en suspensión antes de iniciar el análisis, se puede considerar como el
material suelto en el lecho del río o de la zona de interés.

PROCESO

Ejecución de un cálculo
Para lanzar un cálculo, en primer lugar se deben fijar los parámetros de cálculo, o datos
del problema, mediante las distintas pestañas del menú “Datos del problema”. En el
siguiente sub apartado se explican las distintas opciones de dicho menú.
El cálculo se lanza con el menú Calcular. A través de dicho menú es posible lanzar
directamente el cálculo (“Calcular > Calcular”) o acceder a la ventana de cálculo (“Calcular
> Ventana de cálculo”). Esta segunda opción permite tener un mejor control del proceso
de cálculo, o acceder al archivo que muestra el estado del proceso (Botón Ver salida),
mientras el cálculo se está ejecutando.
La ventana de información del proceso también se abre, durante el proceso o cuando éste
ha finalizado, con el menú “Calcular > Ver información del proceso”. Es conveniente
consultar dicha información para detectar posibles avisos de errores.
Parámetros del cálculo
Se accede a través del menú “Datos > Datos del problema”. Los más representativos son
los parámetros de tiempo, los de cálculo, los de resultados y transporte de sedimentos.

Parámetros de tiempo
Simulación: Con el botón Nueva/Continuar se puede escoger entre lanzar una simulación
nueva o continuar una que ya se ha ejecutado hasta un cierto instante.
Incremento de tiempo máximo: Iber ajusta automáticamente el incremento de cálculo para
satisfacer la condición de Courant. Adicionalmente es posible fijar un incremento de
tiempo de cálculo máximo.
Instante inicial: Valor del instante de tiempo de inicio del cálculo.
Tiempo máximo de simulación: Valor del instante de tiempo final del cálculo.
Intervalo de resultados: Fija el incremento de tiempo entre instantes de escritura de
resultados.

Parámetros de cálculo (Pestaña General)

Con esta pestaña se fijan algunos parámetros de cálculo y del esquema numérico.
Número de procesadores: Iber puede lanzar un cálculo paralelo con el número de
procesadores que se desee. Si se indica un número de procesadores mayor al existente
Iber utilizará el máximo de procesadores posibles.
Función limitadora de flujo: Permite escoger entre un esquema numérico de primer orden,
o uno de segundo orden con diferentes limitadores de flujo: Minmod, Superbee o Van
Leer
CFL: Valor del número de Courant-Friedrichs-Levy para conseguir un esquema numérico
estable.
Límite seco-mojado: Umbral para considerar que un elemento está seco y no se realice
ningún cálculo hidrodinámico en él.
Límite almacenamiento elementos secos: Cuando en un instante de tiempo un elemento
tiene un calado negativo, Iber por defecto lo considera seco, y guarda en memoria el valor
de dicho calado negativo. Para volver a estar mojado, el elemento debe llenar antes este
“calado negativo”. Es posible fijar el valor de este “caldo negativo” o incluso limitarlo a
cero, con ello se ganará algo de precisión en el proceso de mojado y secado, pero se
reducirá el incremento de tiempo de cálculo y aumentará, por lo tanto, el tiempo de cálculo
total.
Fricción en las paredes: Permite decidir si se considera, o no, la rugosidad en las paredes
verticales del contorno, así como el valor de dicha rugosidad.

Selección de resultados
Iber sólo creara archivos de resultados para los resultados seleccionados en el menú,
esto quiere decir que se deben seleccionar las diferentes variables hidráulicas de interés
con la finalidad de visualizar y exportar dichos resultados.

Selección de un modelo de transporte de sedimentos

La pestaña Sedimentos permite activar o desactivar los módulos de cálculo de transporte

de sedimentos por arrastre de fondo y en suspensión, así como fijar los parámetros de los
mismos.
Tanto para transporte en suspensión, como para transporte de fondo, es posible activar o
desactivar el modelo de avalancha.
Es posible fijar un instante de inicio de cálculo del transporte de sedimentos distinto al del
cálculo hidrodinámico. De esta manera se puede, por ejemplo, obviar los procesos de
erosión y sedimentación hasta que las condiciones de velocidad y calado se estabilicen.

Entre los modelos de transporte en suspensión se encuentran:

 Van Rijn
 Smith McLean
 Ariathurai

En cuanto a los modelos de transporte de fondo los siguientes:

 Meyer-Peter & Müller

 Van Rijn
 Engelund - Fredsoe
 Ad hoc
 Engelund - Hasen

POS-PROCESO

Visualización de resultados
Iber dispone de multitud de opciones para visualizar y analizar los resultados, personalizar los
colores, las leyendas, mostrar etiquetas de valores, etc. Una opción para visualizar resultados es
mediante el menú “Ventana > ver resultados”. La ventana que se abre permite acceder a los
distintos grupos de resultados (o Análisis) que se pueden visualizar en cada instante de tiempo,
como áreas coloreadas, áreas coloreadas suavizadas, o vectores.

Dentro del mismo menú, con “Ventana > animar” se pueden realizar animaciones de los resultados.

Creación de gráficos
Con los distintos botones de creación de gráficos de la barra de visualización de
resultados (botones evolución de punto, gráfico de línea, gráfico punto y gráfico de borde)
es posible crear gráficos de utilidad para el análisis de resultados. Los gráficos creados de
pueden gestionar con la ventana de gráficos (“Ventana > ver gráficos”) y se pueden
exportar para abrirlos, por ejemplo, en Excel con “Archivo > Exportar > Gráfico”.
Exportar resultados en formato raster
Es la opción más útil para visualizar los resultados realizados por Iber, los cuales son
transportados a un SIG. Iber permite exportar la mayoría de resultados en formato ASCII
grid de Arc Info a través del menú “Iber_Herramientas > Resultados en Grid”. Sólo se
debe seleccionar un resultado, introducir el valor de un instante de tiempo y un tamaño de
celda. Iber creará una carpeta dentro de la carpeta del proyecto actual (carpeta Grids of
results), y dentro de la misma los archivos en formato ASCII de los resultados
seleccionados en el instante de tiempo más próximo al solicitado.