Iber v3. Manual de referencia e interfaz de usuario de las nuevas implementaciones

<!-- metadata commented in wiki content


==Iber v3. Manual de referencia e interfaz de usuario de las nuevas implementaciones==

Marcos Sanz-Ramos<sup>1</sup>, Luis Cea<sup>2</sup>, Ernest Bladé<sup>1</sup>, David López-Gómez<sup>3</sup>, Esteban Sañudo<sup>3</sup>, Georgina Corestein<sup>1</sup>, Gonzalo García-Alén<sup>2</sup>, José Luis Aragón-Hernández<sup>1,4</sup>

<sup>1</sup>Flumen Institute, Universitat Politècnica de Catalunya (UPC - Barcelona TECH) – International Centre of Numerical Methods in Engineering (CIMNE), 08034 Barcelona, Spain

<span id='_GoBack'></span><sup>2</sup>Universidade da Coruña, Water and Environmental Engineering Group, Department of Civil Engineering, School of Civil Engineering, 15071 A Coruña, Spain

<sup>3</sup>Centro de Estudios Hidrográficos, Centro de Estudios y Experimentación de Obras Públicas, 28005 Madrid, Spain

<sup>4</sup>División de Ingenierías Civil y Geomática, Facultad de Ingeniería, Universidad Nacional de México (UNAM), Ciudad de México 04510, México

''Correspondence to'': [mailto:marcos.sanz-ramos@upc.edu marcos.sanz-ramos@upc.edu]

<span id='_Toc108335426'></span>
-->
==Resumen==

Iber es un modelo numérico bidimensional de simulación de flujo turbulento en lámina libre en régimen no-permanente, y de procesos medioambientales en hidráulica fluvial. Desde el lanzamiento de la primera versión, que incluía un motor de cálculo hidrodinámico completamente acoplado con procesos de transporte de sedimento y turbulencia, Iber ha ido evolucionando hasta convertirse en una herramienta de modelización de flujo en lámina libre de procesos ambientales de elevada complejidad. En este documento se presentan los desarrollos realizados para la versión 3, donde los avances son principalmente en cuatro líneas de investigación: un nuevo módulo de drenaje urbano, un avance significativo en las capacidades del módulo de hidrología, un nuevo módulo de erosión de suelos, y un nuevo módulo para el cálculo del transporte de sedimentos considerando material no uniforme (mezclas). Asimismo, todos los trabajos han ido acompañados por una tarea transversal de mejora de la interfaz, tanto para los módulos ya existentes, como la creación de nuevas ventanas y menús para los nuevos módulos con el objetivo de mejorar todo el flujo de trabajo.

'''Palabras clave''': drenaje urbano, hidrología, erosión de suelos, granulometría no-uniforme, Iber

<span id='_Toc108335427'></span>

==Abstract==

Iber is a two-dimensional hydraulic model for the simulation of free surface flow in rivers and estuaries, and the simulation of environmental processes in fluvial hydraulics. Since the release of the first version of Iber, which included a hydrodynamic calculation engine fully coupled with sediment transport processes and turbulence, it has evolved to become a free surface flow modelling tool for highly complex environmental processes. This document presents the developments made for version 3, where the advances are applied mainly in four current research lines: a new urban drainage module, a significant advance in the capabilities of the hydrological process module, a new soil erosion module, and a new module for calculating sediment transport considering non-uniform material (mixtures). Likewise, all the work has been accompanied by a cross-cutting task of improving the interface, both for existing modules and the creation of new windows and menus for new modules aiming to improve the whole workflow.

'''Keywords''': urban drainage, hydrology, soil erosion, non-uniform granulometry, Iber

<span id='_Toc108335428'></span>

==1 Introducción==

Iber es un modelo numérico de simulación de flujo turbulento en lámina libre en régimen no-permanente, y de procesos medioambientales en hidráulica fluvial [<span id='cite-_Bib001'></span>[[#_Bib001|1]],<span id='cite-_Bib002'></span>[[#_Bib002|2]]]. El modelo fue desarrollado en primera instancia en coordinación por los grupos universitarios Flumen, ahora Instituto Flumen, de la Universitat Politècnica de Catalunya, el grupo GEAMA (Grupo de Enxeñería da Agua e do Medio Ambiente) de la Universidade da Coruña, el Centre Internacional de Mètodes Numèrics en Enginyeria (CIMNE), y el Laboratorio de Hidráulica del Centro de Estudios Hidrográficos del CEDEX.

El modelo Iber se presentó oficialmente (Iber 1.0) en junio de 2010. Entonces Iber contaba con tres módulos de cálculo principales: un módulo hidrodinámico, un módulo de turbulencia y un módulo de transporte de sedimentos (considerando sedimento de granulometría uniforme). Desde entonces, Iber ha incorporado nuevos módulos adicionales como el módulo de rotura de presas y balsas [<span id='cite-_Bib003'></span>[[#_Bib003|3]]], el módulo de calidad de aguas [<span id='cite-_Bib004'></span>[[#_Bib004|4]],<span id='cite-_Bib005'></span>[[#_Bib005|5]]], el módulo de hábitat hidráulico [<span id='cite-_Bib006'></span>[[#_Bib006|6]]], y mejorando también las capacidades de los módulos ya existentes [<span id='cite-_Bib007'></span>[[#_Bib007|7]]], así como la amigabilidad de la interfaz, y el tiempo de cálculo.

En este documento se presentan los desarrollos realizados para la versión 3, donde los avances son principalmente en cuatro líneas de contenido: un nuevo <span id='cite-_Ref106261111'></span>[[#_Ref106261111|Módulo de Drenaje Urbano]], un avance significativo en las capacidades del <span id='cite-_Ref106262222'></span>[[#_Ref106262222|Módulo de Hidrología Distribuido]], un nuevo <span id='cite-_Ref106263333'></span>[[#_Ref106263333|Módulo de Erosión de Suelos''']], y un nuevo módulo para el cálculo del <span id='cite-_Ref106264444'></span>[[#_Ref106264444|Transporte de Sedimentos considerando material no uniforme (mezclas)]]. Asimismo, todos los trabajos han ido acompañados por una tarea transversal de <span id='cite-_Ref106260000'></span>[[#_Ref106260000|mejora de la interfaz]], tanto para los módulos ya existentes, como la creación de nuevas ventanas y menús para los nuevos módulos. Por ello, antes de la descripción de cada uno de los módulos, se describen los desarrollos de interfaz generales.

Los apartados de este documento, que correspondientes a cada uno de los nuevos módulos, además del apartado de descripción general de la interfaz, se estructuran en dos partes diferenciadas: fundamentos teóricos e interfaz del usuario del propio módulo. En la presentación de las bases teóricas se describe, en primer lugar, los fenómenos objeto de los nuevos módulos, tanto su descripción física como las ecuaciones de gobierno que permiten describirlos, y, en segundo lugar, se relatan los esquemas numéricos para la resolución de las ecuaciones implementados en Iber. El apartado de interfaz es una descripción de los nuevos menús y ventanas desarrollados en Iber para cada uno de los nuevos módulos, con lo que por un lado se presentan los desarrollos de mejora de la interfaz propiamente dicha, pero por el otro se dispone ya de una guía para el uso de las nuevas herramientas.

<span id='_Toc108335430'></span>

===1.1 Breve descripción de los nuevos módulos===

Iber ha pasado de ser un modelo hidrodinámico y de procesos morfodinámicos, a un modelo numérico que, estando basado en la resolución de las ecuaciones de aguas profundas, permite simular múltiples procesos ambientales de base hidrodinámica. La incorporación de nuevos módulos, así como la mejora de los ya existentes, se ha llevado a cabo de forma totalmente acoplada, siendo el nexo de unión entre todos ellos el módulo hidrodinámico. La <span id='cite-_Ref106264389'></span>[[#_Ref106264389|Fig. 1]] representa, de forma esquemática, las interrelaciones entre los distintos módulos que incorpora Iber.

<span id='_Ref106264389'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image1.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 1. Interrelaciones entre los distintos módulos principales de cálculo de Iber v3.'''</span>


<span id='_Toc108335431'></span>

====1.1.1 Módulo de drenaje urbano====

El nuevo <span id='cite-_Ref106261111'></span>[[#_Ref106261111|Módulo de Drenaje Urbano]] se ha construido con la finalidad de ser complementario del modelo Iber ya existente, de cálculo bidimensional en superficie, para poder incorporar los procesos asociados al drenaje de agua en zona urbana a través de la red de alcantarillado.

El elemento principal de dicho módulo es el cálculo, en una dimensión, de la propagación del flujo a través de los conductos cerrados que conforman la red de alcantarillado, ya sea en lámina libre, o con flujo a presión. La consideración de flujo a presión se ha implementado a través de la metodología de la ranura de Preissmann.

Un proceso del drenaje urbano es el intercambio de agua entre la superficie y la red, ya sea solo de captación, a través de las rejas de alcantarillado, o de salida de la red a superficie en episodios de grandes lluvias. Ambos fenómenos han sido incorporados en Iber a través de dos metodologías distintas: considerando las rejas como orificios, o mediante la metodología de evaluación de la eficiencia hidráulica de las rejas desarrollada en los laboratorios del Instituto Flumen, en la UPC.

Para la simulación de las uniones de colectores, que también pueden funcionar en lámina libre o a presión, se ha optado por considerarlos como elementos puntuales de superficie, y calcularlos con un esquema bidimensional, semejante al utilizado por Iber en superficie.

En la red de alcantarillado se ha implementado una opción para poder considerar elementos especiales de la propia red (por ejemplo, depósitos de retención, uniones complejas de colectores, etc.), pudiéndose simular estos en un esquema bidimensional. Asimismo, este módulo se puede emplear para realizar la simulación únicamente en superficie, tanto en superficie como en la red, o tan solo en la red de drenaje.

<span id='_Toc108335432'></span>

====1.1.2 Nuevos desarrollos en el módulo de hidrología====

El <span id='cite-_Ref106262222'></span>[[#_Ref106262222|Módulo de Hidrología Distribuido]] de Iber permite calcular la transformación lluvia-escorrentía en toda una cuenca, haciendo posible su utilización como modelo hidrológico distribuido en cuencas de pequeño y mediano tamaño. El módulo está basado en las ecuaciones hidrodinámicas de aguas poco profundas bidimensionales, las cuales se utilizan para resolver el movimiento de la escorrentía superficial sobre el terreno, al que se le han añadido una serie de nuevos procesos, como son: el aporte de agua de lluvia, la infiltración de agua en el terreno, el flujo subsuperficial, la evapotranspiración y la percolación.

El módulo hidrológico es completamente bidimensional y distribuido en el espacio, y trabaja sobre la misma malla de cálculo que el módulo hidrodinámico. La base del módulo hidrológico es el módulo hidrodinámico, al que se le han añadido los procesos de aporte de agua por precipitación y la infiltración del agua en el terreno. De esta forma, el exceso de precipitación se mueve por el terreno y por los cauces de acuerdo con las ecuaciones de aguas someras bidimensionales. El agua infiltrada pasa a formar parte del subsuelo, en donde se puede realizar un balance de masa de agua distribuido en el espacio a nivel de elemento de la malla de cálculo. En dicho balance de agua en el terreno se consideran los procesos de infiltración, evapotranspiración, percolación y flujo en el subsuelo. Los tres primeros implican un movimiento de agua en la columna vertical, mientras que el flujo en el subsuelo se considera bidimensional y paralelo al terreno. Todos estos procesos pueden activarse o desactivarse por el usuario.

<span id='_Toc108335433'></span>

====1.1.3 Módulo de erosión en laderas por escorrentía====

El nuevo <span id='cite-_Ref106263333'></span>[[#_Ref106263333|Módulo de Erosión de Suelos]] de Iber permite determinar las zonas en las que se produciría erosión del suelo en una cuenca y cuantificar la misma, así como evaluar la evolución temporal y espacial de las concentraciones de sólidos en suspensión a lo largo de los cauces de la red fluvial.

El módulo puede trabajar con distintos niveles de complejidad de la estructura del suelo, mediante la activación o no de los diferentes procesos de erosión considerados. De esta manera, se pueden considerar una única capa de suelo erosionado (no cohesivo), o dos capas de suelo, una correspondiente a la matriz original del suelo (cohesiva) y otra correspondiente al suelo erosionado (no cohesivo). Las partículas de ambas capas pueden ponerse en movimiento por el impacto de la lluvia y por la fricción que genera el agua con el fondo. Una vez en movimiento, los procesos de transporte de sólidos considerados son por carga de fondo y en suspensión, si bien en la mayor parte de los casos el transporte en suspensión de partículas finas será predominante en las laderas de las cuencas. Asimismo, se consideran diferentes clases o tipos de sedimento, por lo que permite trabajar con granulometrías distribuidas si el usuario lo considera necesario.

<span id='_Toc54974079'></span><span id='_Toc90987590'></span><span id='_Toc108335434'></span>

====1.1.4 Módulo de transporte de sedimento no uniforme (mezclas)====

Desde el principio Iber incorporó un módulo de transporte de sedimentos por arrastre de fondo, y un módulo de transporte en suspensión, pero considerando la hipótesis, o simplificación, de material uniforme. En esta versión de Iber se ha implementado y validado un nuevo <span id='cite-_Ref106264444'></span>[[#_Ref106264444|Módulo de Transporte de Sedimentos con granulemtría no uniforme]], o granulometrías extendidas. Ello proporciona no solo mayor precisión en los resultados, sino también permite la simulación numérica de procesos que antes no se podían contemplar, como son: el acorazamiento del lecho, erosiones diferenciales por efecto de las distintas granulometrías, o la clasificación del material del fondo por zonas según la capacidad de arrastre.

Los desarrollos se basan en la teoría de la capa activa, y se ha realizado su implementación numérica mediante el método de los volúmenes finitos. De esta manera, el fondo de un cauce se puede dividir en distintos estratos, cada uno de ellos con una granulometría distinta representada por una serie de clases y su fracción. La ecuación de estimación del caudal sólido se hace para cada una de las clases, considerando el fenómeno de escondimiento para la interacción entre los distintos tamaños de grano.

<span id='_Toc108335435'></span>
<span id='_Ref106260000'></span>

==2 Interfaz de Iber 3.0==

<span id='_Toc108335436'></span>

===2.1 Generalidades===

<span id='_Toc108335437'></span>

====2.1.1 Apariencia====

La versión actual de Iber, al igual que versiones anteriores, se ha embebido dentro de la interfaz de pre- y post-proceso GiD [8–10]. Para ello se ha empleado la versión 15 de GiD, que incorpora múltiples ventajas y mejoras con respecto a versiones anteriores (más información en [http://www.gidhome.com www.gidhome.com]).

<span id='_Ref106265212'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image2.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 2. Visión general de la interfaz de Iber v3.'''</span>


La apariencia general del programa es, a priori, muy similar a la versión precedente (<span id='cite-_Ref106265212'></span>[[#_Ref106265212|Fig. 2]]). El principal cambio en la interfaz, concretamente en la parte de pre-proceso, ha sido la simplificación de los menús y el agrupamiento de los diferentes módulos de cálculo en Plug-ins como se verá en el apartado <span id='cite-_Ref106265695'></span>[[#_Ref106265695|2.1.2]].

En versiones anteriores de Iber los menús y submenús, así como las opciones generales de cada uno de los módulos de cálculo, se mostraban por defecto, aunque el/la usuario/a no los emplease durante la generación y cálculo de un modelo. Esto podía entorpecer el flujo de trabajo y, además, generar confusión a la hora de definir los parámetros.

Puesto que Iber v3 incluye nuevos módulos cálculo, se ha creído conveniente agrupar todas y cada una de las opciones, menús y parámetros de cada módulo de cálculo. Esto ha redundado en una simplificación significativa de toda la interfaz, mejorando el flujo de trabajo para que el/la usuario/a pueda generar un modelo de forma eficaz, eficiente y bajo una experiencia lo más amigable posible, sin interferencias ni distracciones.

Así, por defecto, tan solo aparecen los menús, opciones y parámetros para el cálculo de la hidrodinámica, que es el módulo básico para el resto de procesos que incorpora Iber. Como se puede ver en la <span id='cite-_Ref106265581'></span>[[#_Ref106265581|Fig. 3]]a, el menú ‘Datos’ se ha simplificado enormemente mostrando, únicamente, aquellos submenús relacionados con los parámetros generales del modelo (‘Datos del Problema’), la ‘Hidrodinámica’ (condiciones de contorno, iniciales, internas, etc.) y aquellos parámetros que guardan relación con los esfuerzos tensionales (rugosidad del fondo, turbulencia y viento). En la <span id='cite-_Ref106265581'></span>[[#_Ref106265581|Fig. 3]]b se muestra la ventana ‘Datos del Problema’ que, al igual que el menú ‘Datos’, por defecto tan solo muestra opciones generales del módulo hidrodinámico.

<span id='_Ref106265581'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image3.png|650px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 3. (a) Disposición por defecto del menú ‘Datos’. (b) Ventana de ‘Datos del Problema’ por defecto.'''</span>


Asimismo, ha habido pequeños cambios en la interfaz que también modifican ligeramente el flujo de trabajo (como se verá en el apartado <span id='cite-_Ref106265695'></span>[[#_Ref106265695|2.1.2]]). Es el caso de la herramienta ‘Brecha’, antes ubicada en el menú ‘Datos’, que ha sido trasladada al menú ‘Herramientas Iber’.

<span id='_Ref106265695'></span><span id='_Toc108335438'></span>

====2.1.2 Flujo de trabajo====

El flujo de trabajo de la versión 3 de Iber sigue siendo el mismo que en versiones anteriores, pero ahora si un/a usuario/a desea realizar cálculos con alguno de los módulos adicionales a la hidrodinámica, lo tendrá que activar. Para ello, deberá emplear la ventana de Plug-ins (<span id='cite-_Ref106266172'></span>[[#_Ref106266172|Fig. 4]]a), que es accesible desde '''Herramientas Iber >> Plug-ins''' de la barra de herramientas de Iber (<span id='cite-_Ref106266172'></span>[[#_Ref106266172|Fig. 4]]b).

<span id='_Ref106266172'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image4.png|400px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 4. (a) Ventana de selección de plug-ins. (b) Activación de los plug-ins mediante el atajo en la barra de herramientas de Iber.'''</span>


Esta ventana permite a usuarios/as seleccionar, en primera instancia, el método de paralelización cálculo mediante tarjetas gráficas (GPU) y, en segundo lugar, el módulo de cálculo deseado. La elección de un método de paralelización del cálculo mediante GPUs es opcional y, por el momento, está condicionado al tipo o tipos de módulos de cálculo que se desee emplear.

<span id='_Ref106266529'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image5.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 5. Ejemplo de activación del módulo de ‘Transporte de Sedimentos’: (a) cambios en la ventana de activación de Plug-ins; (b) incorporación de los menús y submenús en ‘Datos’; (c) actualización de la ventana ‘Datos del Problema’ con las nuevas opciones (c).'''</span>


Una vez seleccionado el módulo de cálculo, se aplicarán los cambios a la interfaz para mostrar las funcionalidades de dicho módulo. A modo de ejemplo se va a seleccionar el módulo de ‘Transporte de Sedimentos’. De este modo, automáticamente se desactivarán los módulos que no son compatibles con el que se ha activado (<span id='cite-_Ref106266529'></span>[[#_Ref106266529|Fig. 5]]a) (esto de describe en el apartado <span id='cite-_Ref106266611'></span>[[#_Ref106266611|2.2]]). Además, el menú ‘Datos’ (<span id='cite-_Ref106266529'></span>[[#_Ref106266529|Fig. 5]]b) y la ventana ‘Datos del Problema’ (<span id='cite-_Ref106266529'></span>[[#_Ref106266529|Fig. 5]]c) ahora incorporarán la opciones, parámetros y funcionalidades propias de dicho módulo. El flujo de trabajo a partir de aquí es idéntico al de versiones anteriores.

<span id='_Ref106266611'></span><span id='_Toc108335439'></span>

===2.2 Compatibilidad entre módulos de cálculo===

Iber es un modelo en continuo desarrollo en el que, a lo largo de su trayectoria, ha ido incorporando nuevos módulos de cálculo y procesos que amplían y agilizan la computación de procesos hidrodinámicos, así como otros íntimamente relacionados con el cálculo hidrodinámico en lámina libre.

Si bien todos los módulos de cálculo dependen de la hidrodinámica y, por tanto, están completamente integrados con ella, no todos son compatibles entre sí. Esto sucede, por ejemplo, con el módulo de ‘Hábitat’, el cual no es compatible con el módulo de ‘Drenaje Urbano’ ni con el de ‘Erosión de suelos’ porque el hábitat fluvial no guarda relación directa con el drenaje urbano y los procesos de erosión de suelos en cuencas. Otro caso sería el de ‘Transporte de Sedimentos’, que tampoco es compatible con los dos módulos anteriores. El motivo por el cual no es compatible con ‘Erosión de suelos’ es porque este módulo ya incorpora su propio submódulo de transporte de sedimentos; mientras que no es compatible con ‘Drenaje Urbano’ porque todavía no se ha integrado con este módulo.

Asimismo, desde la versión 2.5, Iber incorpora la posibilidad de realizar el cálculo hidrodinámico y de procesos hidrológicos mediante dos posibilidades: código secuencial parcialmente paralelizado mediante CPU; y cálculo paralelo mediante GPU (tarjetas gráficas) con arquitectura CUDA. El cálculo con CPU incorpora todos y cada uno de los últimos módulos de cálculo y funcionalidades, mientras que el cálculo con GPU incorpora algunos de los módulos de cálculo y funcionalidades del cálculo en CPU. Es por ello que cuando se activa en Plug-ins la opción de paralelización ‘Plus’ o ‘R-Iber’, este último recientemente incorporado fruto del presente desarrollo, no es posible seleccionar diversos módulos de cálculo ni tampoco utilizar otras funcionalidades que sí están disponibles en la versión secuencial.

<span id='_Toc108335440'></span>

===2.3 Nuevas herramientas y funcionalidades mejoradas===

Otros cambios relacionados con el flujo de trabajo, pero que a la vez se trata de nuevas funcionalidades y herramientas, es la consideración de flujo bajo a presión en tramos cubiertos y la eliminación de la herramienta ‘Brecha’ del menú ‘Datos’. De este modo, para ‘Cubiertas’, el/la usuario/a deberá proceder como en versiones anteriores, imponiendo la cubierta donde considere, pero ahora deberá activar su cálculo a través de su correspondiente pestaña en ‘Datos del Problema’ (<span id='cite-_Ref106265581'></span>[[#_Ref106265581|Fig. 3]]b). Aquí tendrá, además, la posibilidad de calcular cubiertas mediante la metodología tradicional de la ranura de Preissmann 2D [11,12] o el nuevo método basado en la teoría de cubiertas elásticas [13]. En lo que respecta a la herramienta ‘Brecha’ (<span id='cite-_Ref106267019'></span>[[#_Ref106267019|Fig. 6]]a), así como pasa con ‘Alcantarillas’ (<span id='cite-_Ref106267019'></span>[[#_Ref106267019|Fig. 6]]b), el/la usuario/a no tiene que activarla, tan solo tiene que definirla como se ha hecho hasta ahora en versiones anteriores de Iber y decidir si desea o no calcularla.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image6.png|450px]] </div>

<span id='_Ref106267019'></span><span style="text-align: center; font-size: 75%;">'''Fig. 6. Integración de la opción de cálculo en: (a) ‘Brechas’; (b) ‘Alcantarillas’.'''</span>

Asimismo, ha habido un ligero cambio en cómo el/la usuario/a solicita los resultados hidrodinámicos, así como aquellos que dependen de la hidrodinámica. Éstos se siguen indicando desde la pestaña ‘Resultados’ de la ventana ‘Datos del Problema’, sin embargo, en la versión 3 de Iber, se han ordenado y agrupado para una mayor agilidad y comodidad tal como muestra la <span id='cite-_Ref106267279'></span>[[#_Ref106267279|Fig. 7]]a. En la opción ‘Hidrodinámica’ se encuentran resultados propios de la resolución de las ecuaciones de flujo en lámina libre y que, por defecto, están casi todos activados. Los términos de ‘Peligrosidad’ se han separado según los criterios del RD 9/2008 [14] y del ACA 2003 [15], agilizando así el cálculo, y, además, la pestaña ‘Peligrosidad personalizada’ se ha integrado en este apartado. Asimismo, se ha creado una nueva opción para extraer resultados, llamada ‘Sondas’, que permite extraer valores hidrodinámicos (calado, velocidad, etc.) de puntos concretos del modelo a intervalos de tiempo diferentes a los de la escritura de resultados (<span id='cite-_Ref106267279'></span>[[#_Ref106267279|Fig. 7]]b). Esta opción es especialmente útil para modelos cuyo tiempo de cálculo es muy extenso y en los que no es necesario escribir todos los resultados con pasos de tiempo reducidos.

<span id='_Ref106267279'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image7.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 7. Organización de la pestaña ‘Resultados’ y nuevas opciones: (a) vista general; (b) nueva herramienta ‘sondas’.'''</span>


En cuanto al post-proceso, la principal novedad está en una mejora sustancial del tiempo de respuesta y de cálculo para con la generación de rasters de resultados ('''Herramientas Iber >> Resultados a Raster'''). Además, en la actual versión se ha posibilitado cambiar el tipo de análisis, seleccionar en qué instante de tiempo se desea extraer el resultado o resultados, y se han mejorado las opciones de exportación para que sean más rápidas y eficientes. De este modo, un/a usuario/a puede seleccionar uno o varios resultados de un análisis concreto, exportarlos en el paso de tiempo que está activo, en todos o en uno concreto que no se esté visualizando, con un tipo u otro de interpolación y un tamaño de celda determinado. Los resultados se guardarán dentro de la carpeta del modelo preservando la clasificación por medio del tipo de análisis y paso de tiempo. Por ejemplo, si se seleccionan ‘Velocidad’ y ‘Cota de agua’ (<span id='cite-_Ref106267963'></span>[[#_Ref106267963|Fig. 8]]a) y se aplica, en la carpeta del modelo se crearán los ''n'' rasters de velocidad y cota de agua asociados a los ''n'' pasos de tiempo analizados (<span id='cite-_Ref106267963'></span>[[#_Ref106267963|Fig. 8]]b).

<span id='_Ref106267963'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image8.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 8. (a) Ejemplo de la herramienta ‘Exportar rasters de resultados’. (b) conjunto de rasters de resultados exportados.'''</span>


Finalmente, para el módulo de Transporte de Sedimentos para mezclas se ha diseñado una nueva funcionalidad que permite conocer la distribución granulométrica de la muestra. Para ello, es necesario ir a '''Herramientas Iber >> Transporte de Sedimentos >> Granulometría'''. La <span id='cite-_Ref106268058'></span>[[#_Ref106268058|Fig. 9]] muestra, a modo de ejemplo, la granulometría en un instante de tiempo concreto en los elementos 23, 28 y 33. Como se observa, la gráfica está en escala logarítmica.

<span id='_Ref106268058'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image9.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 9. Visualización de la distribución granulométrica de una muestra.'''</span>


<span id='_Toc108335441'></span>
<span id='_Ref106261111'></span>

==3 Módulo de drenaje urbano==

<span id='_Toc108335442'></span>

===3.1 Fundamentos teóricos===

El flujo de agua en zona urbana presenta una complejidad considerable respecto el flujo en zona no urbana. Por un lado, el agua puede circular por la superficie (sistema de drenaje mayor) y/o por la red de colectores (sistema de drenaje menor). Por otro lado, el flujo puede ser en lámina libre o a presión. En el sistema de drenaje menor, en los colectores y en las uniones de una red de alcantarillado ocurren fenómenos hidráulicos complejos. La interacción entre el sistema de drenaje mayor (tejados, calles, áreas verdes, áreas grises, canales naturales y canales artificiales, etc.) y el sistema de drenaje menor (colectores, uniones, compuertas, etc.), conocida como drenaje urbano dual, se realiza en los puntos físicos correspondientes a la localización espacial de los elementos de captación (fuentes y sumideros).

<span id='_Toc54793149'></span>Desde el punto de vista numérico, para la resolución de las ecuaciones de flujo en lámina libre en una dimensión, las ecuaciones de flujo en presión (tanto en una dimensión como en dos) y las ecuaciones de aguas poco profundas en dos dimensiones, se emplean esquemas numéricos en volúmenes finitos, capaces de calcular discontinuidades y frentes de onda en la solución, sin necesidad de técnicas adicionales [16].

<span id='_Toc108335443'></span>

====3.1.1 Ecuaciones de la hidrodinámica====

Para el flujo en zona urbana se utilizan las ecuaciones de Saint Venant, tanto en una dimensión como en dos dimensiones. Las ecuaciones bidimensionales se utilizan sobre todo en superficie, pero también para elementos que forman parte de la red de colectores, como pueden ser depósitos de retención o uniones de colectores. Las ecuaciones unidimensionales se utilizan principalmente en los colectores.

<span id='_Toc90987600'></span><span id='_Toc108335444'></span>

====3.1.2 Ecuaciones de Saint Venant en 2D====

En el módulo de drenaje urbano, los cálculos en dos dimensiones se realizan en base a las ecuaciones de aguas someras, o ecuaciones de Saint Venant 2D, que se describen en detalle en el Manual de Referencia Básico. Se repiten aquí por completitud, incluyendo los términos que representan las componentes de precipitación e infiltración y con la fórmula de Manning para representar la fricción, siendo el sistema de ecuaciones el siguiente:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}\frac{\partial h}{\partial t}+\frac{\partial {q}_{x}}{\partial x}+\frac{\partial {q}_{y}}{\partial y}=R-f\\\frac{\partial {q}_{x}}{\partial t}+\frac{\partial }{\partial x}\left( \frac{{q}_{x}^{2}}{h}+g\frac{{h}^{2}}{2}\right) +\frac{\partial }{\partial y}\left( \frac{{q}_{x}{q}_{y}}{h}\right) =-gh\frac{\partial {z}_{b}}{\partial x}-g\frac{{n}^{2}}{{h}^{\frac{7}{3}}}\left| \mathit{\boldsymbol{q}}\right| {q}_{x}\\\frac{\partial {q}_{y}}{\partial t}+\frac{\partial }{\partial x}\left( \frac{{q}_{x}{q}_{y}}{h}\right) +\frac{\partial }{\partial y}\left( \frac{{q}_{y}^{2}}{h}+g\frac{{h}^{2}}{2}\right) =-gh\frac{\partial {z}_{b}}{\partial y}-g\frac{{n}^{2}}{{h}^{\frac{7}{3}}}\left| \mathit{\boldsymbol{q}}\right| {q}_{y}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref106269045'></span>(1)
|}

En ellas, <math display="inline">h</math> es el calado de agua, <math display="inline">{q}_{x}</math> y <math display="inline">{q}_{y}</math> son las dos componentes del caudal unitario, <math display="inline">\vert \mathit{\boldsymbol{q}}\vert</math>  es su módulo, <math display="inline">R</math> es la intensidad de precipitación, <math display="inline">f</math> es la tasa de infiltración, <math display="inline">{z}_{b}</math> es la cota del fondo, <math display="inline">g</math> es la aceleración de la gravedad y <math display="inline">n</math> es el coeficiente de Manning.

Estas ecuaciones son las que se resuelven mediante el módulo hidrodinámico de Iber desde sus inicios, y se encuentran ampliamente documentadas, tanto las propias ecuaciones como los métodos de resolución (p.ej. en Toro [17]).

<span id='_Toc108335445'></span>

====3.1.3 Ecuaciones de Saint Venant en 1D====

La deducción de las ecuaciones de Saint Venant en 1D se realiza directamente aplicando las leyes de conservación de la masa y cantidad de movimiento a un volumen de control considerando una sección arbitraria y canal no prismático [18]. El resultado es un sistema de ecuaciones que, en forma conservativa, se puede escribir como sigue:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}\frac{\partial A}{\partial t}+\frac{\partial Q}{\partial x}=0\\\frac{\partial Q}{\partial t}+\frac{\partial }{\partial x}\left( \frac{{Q}^{2}}{A}+g{I}_{1}\right) =g{I}_{2}+gA\left( {S}_{0}-{S}_{f}\right) \end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref106269418'></span>(2)
|}

En las expresiones anteriores, es <math display="inline">A</math> el área de la sección mojada, <math display="inline">Q</math> el caudal circulante, <math display="inline">{S}_{0}</math> la pendiente del terreno, <math display="inline">{S}_{f}</math> la pendiente de fricción o motriz, <math display="inline">{I}_{1}</math> la fuerza debida a la presión del agua en una sección e <math display="inline">{I}_{2}</math> la contribución de las fuerzas de presión del contorno. Las tres últimas se expresan respectivamente como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{S}_{f}=\frac{{v}^{2}{n}^{2}}{{R}_{h}^{4/3}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{1}=\int_{0}^{h}\left( h-\eta \right) b\left( x,\eta \right) d\eta</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(4)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{2}=\int_{0}^{h}\left( h-\eta \right) \frac{\partial b\left( x,\eta \right) }{\partial x}d\eta</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(5)
|}

donde <math display="inline">n</math> es el coeficiente de fricción de Manning, <math display="inline">{R}_{h}^{}</math> el radio hidráulico. Para canales prismáticos, aunque tengan una sección cualquiera, el término <math display="inline">{I}_{2}</math> es idénticamente igual a cero.

<span id='_Toc108335446'></span>

====3.1.4 Esquema numérico para el flujo en colectores====

La aplicación de un método numérico conlleva la discretización espacial del dominio en celdas o elementos de cálculo. La unión de estos elementos forma una malla. Iber, en todos sus módulos, utiliza el Método de los Volúmenes Finitos para resolver las ecuaciones, tal como se detalla en el Manual de Referencia Hidráulico, por lo que cada elemento será a su vez un Volumen Finito. Un volumen finito unidimensional es un segmento, mientras que un volumen finito bidimensional es una superficie (en Iber tienen forma de triángulo o cuadrilátero).

En el módulo de drenaje urbano conviven dominios unidimensionales y bidimensionales y, por lo tanto, también habrá mallas en 1D y mallas en 2D (<span id='cite-_Ref106270044'></span>[[#_Ref106270044|Fig. 10]]). La malla puede ser estructurada (<span id='cite-_Ref106270044'></span>[[#_Ref106270044|Fig. 10]]b) o no estructurada (<span id='cite-_Ref106270044'></span>[[#_Ref106270044|Fig. 10]]c). En las mallas estructuradas, cada elemento de la misma está inequívocamente identificado por unos índices en coordenadas cartesianas. Por su parte, en mallas no estructuradas los elementos no tienen un orden particular. Los elementos de una malla estructurada en 2D por lo general son cuadriláteros, aunque también pueden ser triángulos equiláteros. En cambio, en una malla no estructurada pueden ser triángulos, cuadriláteros o una combinación de éstos. Las mallas no estructuradas tienen la ventaja de adaptarse mejor a geometrías complejas como las que presentan la batimetría y fronteras de estudios hidráulicos. Sin embargo, éstas requieren de una estructura de datos compleja y unos requerimientos de memoria mayores.

<span id='_Ref106270044'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image10.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 10. Tipos de malla: (a) malla en una dimensión; (b) malla estructurada en dos dimensiones; (c) malla no estructurada en dos dimensiones.'''</span>


En los dominios unidimensionales, la discretización empleada para la resolución en el tiempo de las ecuaciones de Saint Venant con el método de los volúmenes finitos se puede escribir de forma genérica como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{U}}}_{i}^{n+1}={\mathit{\boldsymbol{U}}}_{i}^{n+1}-\frac{\Delta t}{\Delta x}\left( {\mathit{\boldsymbol{F}}}_{i+\frac{1}{2}}^{\ast }-\right. </math><math>\left. {\mathit{\boldsymbol{F}}}_{i-\frac{1}{2}}^{\ast }\right) +\frac{\Delta t}{\Delta x}{\mathit{\boldsymbol{H}}}_{i}^{\ast }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(6)
|}

En esta nomenclatura, los subíndices enteros hacen referencia a un volumen finito o elemento de la malla (segmento asociado a una sección transversa), mientras que los subíndices <math display="inline">i+</math><math>\frac{1}{2}</math> e <math display="inline">i-\frac{1}{2}</math> hacen referencia a los puntos de contacto entre segmentos. En Iber el cálculo de los flujos numéricos se realiza con el esquema de Roe [19] y su expresión resulta:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{F}}}_{i+\frac{1}{2}}^{\ast }=\frac{1}{2}\left( {\mathit{\boldsymbol{F}}}_{i}+\right. </math><math>\left. {\mathit{\boldsymbol{F}}}_{i+1}\right) -\frac{1}{2}\left[ \sum _{j=i}^{2}{\tilde{\alpha }}_{j}{\varphi }_{j}{\tilde{e}}_{j}+\right. </math><math>\left. \sum _{j=i}^{2}{\tilde{\gamma }}_{j}signo\left( {\tilde{\lambda }}_{j}\right) {\tilde{\mathit{\boldsymbol{e}}}}_{j}\right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(7)
|}

En la ecuación interior, y las siguientes, los términos con el superíndice <math display="inline">\hat{}</math> son valores promediados en los contornos de los elementos. Los coeficientes constantes o fuerzas de cada onda <math display="inline">{\tilde{\alpha }}_{j}</math>, los valores y vectores propios <math display="inline">{\tilde{\lambda }}_{j}</math> y <math display="inline">{\tilde{\mathit{\boldsymbol{e}}}}_{j}</math> de la matriz jacobiana definida como <math display="inline">\mathit{\boldsymbol{A}}=</math><math>\partial \mathit{\boldsymbol{F}}/\partial \mathit{\boldsymbol{U}}</math>, se pueden expresar respectivamente como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tilde{\alpha }}_{1,2}=\frac{\left( \tilde{c}\mp \tilde{u}\right) \Delta A\pm \Delta Q}{2\tilde{c}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(8)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tilde{\lambda }}_{1,2}=\tilde{c}\mp \tilde{u}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(9)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tilde{\mathit{\boldsymbol{e}}}}_{1,2}=\left[ \begin{matrix}1\\\tilde{c}\pm \tilde{u}\end{matrix}\right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(10)
|}

siendo,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\Delta A={A}_{i+1}-A</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(11)
|}

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\Delta Q={Q}_{i+1}-Q</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(12)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\tilde{u}=\frac{\frac{{Q}_{i}}{\sqrt{{A}_{i}}}+\frac{{Q}_{i+1}}{\sqrt{{A}_{i+1}}}}{\sqrt{{A}_{i}}+\sqrt{{A}_{i+1}}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(13)
|}

La corrección de entropía de Harten y Hyman <math display="inline">{\varphi }_{j}</math>, según Toro [17], es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\varphi }_{j}=\left\{ \begin{matrix}\left| {\tilde{\lambda }}_{j}\right| \, si\, \left| {\tilde{\lambda }}_{j}\right| \geq {\epsilon }_{j}\\{\epsilon }_{j}\, si\, \left| {\tilde{\lambda }}_{j}\right| <{\epsilon }_{j}\end{matrix}\right.</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(14)
|}

donde,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\epsilon }_{j,i+1/2}=max\left[ 0,\left( {\tilde{\lambda }}_{j,i+1/2}-{\tilde{\lambda }}_{j,i}\right) ,\left( {\tilde{\lambda }}_{j,i+1}-\right. \right. </math><math>\left. \left. {\tilde{\lambda }}_{j,i+1/2}\right) \right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(15)
|}

Para reproducir correctamente el salto de las fuerzas de presión, se utiliza una expresión que mantiene el significado físico de la celeridad, entendido como la variación de las fuerzas de presión en cada sección transversal con respecto a la variación del área de flujo en cada sección [20]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tilde{c}}^{2}=\left\{ \begin{matrix}g\frac{\left( {I}_{i+1}-{I}_{i}\right) -{\left( {\Delta {I}_{i}\vert }_{\tilde{A}}\right) }_{i+1/2}}{{A}_{i+1}-{A}_{i}}\, si\, {A}_{i}\not ={A}_{i+1}\\\frac{{c}_{i}^{2}+{c}_{i+1}^{2}}{2};\, o\, bien\, g\frac{\tilde{A}}{\tilde{B}}\, ;\acute{o}\, bien\, \frac{\Delta I}{\Delta A}\, si\, {A}_{i}={A}_{i+1}\end{matrix}\right.</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(16)
|}

Cuando <math display="inline">{A}_{i}={A}_{i+1}</math>, Iber emplea la segunda opción. En las anteriores expresiones:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\tilde{A}=\frac{2{A}_{i}{A}_{i+1}}{{A}_{i}+{A}_{i+1}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(17)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\tilde{B}=\frac{{B}_{i}+{B}_{i+1}}{2}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(18)
|}

La variación de las fuerzas de presión <math display="inline">{\left( {\Delta {I}_{i}\vert }_{\tilde{A}}\right) }_{i+1/2}</math> en volúmenes finitos para un área constante <math display="inline">{A}_{i+1/2}</math> es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\left( {\Delta {I}_{1}\vert }_{\tilde{A}}\right) }_{i+1/2}={\left( {{I}_{i}\vert }_{\tilde{A}}\right) }_{i+1}-</math><math>{\left( {\Delta {I}_{1}\vert }_{\tilde{A}}\right) }_{i}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(19)
|}

donde <math display="inline">{\left( {{I}_{i}\vert }_{\tilde{A}}\right) }_{i}</math> y <math display="inline">{\left( {{I}_{i}\vert }_{\tilde{A}}\right) }_{i+1}</math> son las fuerzas de presión en el volumen finito <math display="inline">i</math> e <math display="inline">i+</math><math>1</math> para un valor del área <math display="inline">\tilde{A}</math>.

Finalmente, el término <math display="inline">{\tilde{\gamma }}_{j}</math> se define como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}\, {\tilde{\gamma }}_{1}=\frac{1}{2\tilde{c}}g\Delta \left( {{I}_{i}\vert }_{\tilde{A}}\right) \\{\tilde{\gamma }}_{2}=-{\tilde{\gamma }}_{1}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(20)
|}

<span id='_Toc108335447'></span>

====3.1.5 Flujo mixto en conductos cerrados====

La ocurrencia simultánea de flujo en lámina libre y flujo en presión en un mismo conducto se conoce comúnmente como flujo mixto. Este fenómeno se puede encontrar, frecuentemente, en los colectores de aguas pluviales, en túneles y conductos de obras de toma de centrales hidroeléctricas, llenado/vaciado de conductos, colectores de almacenamiento, etc.

El proceso de transición de flujo en lámina libre a flujo en presión puede ser inducido, entre otras causas, por la variación de caudales de entrada, bloqueo o reducción de la capacidad del conducto, sumersión de la salida del conducto, operación de elementos de control (compuertas, vertederos, etc.), fallo de estaciones de bombeo, presencia de pozos de caída y falta de aire, y puede producir una variación importante de los calados, cargas de presión, velocidades, entrada o expulsión de aire.

<span id='_Ref106273836'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image11.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 11. Flujo mixto en un conducto: (a) transición gradual; (b) transición brusca.'''</span>


Los frentes de onda generados en la transición de flujo en lámina libre a flujo en presión pueden ser de dos tipos: gradual, que consiste de un frente de onda en lámina libre con condiciones de flujo en lámina libre en ambos lados del frente (<span id='cite-_Ref106273836'></span>[[#_Ref106273836|Fig. 11]]a); o brusco, que es un frente de onda a tubo lleno con condiciones de flujo en lámina libre delante y flujo en presión atrás del mismo (<span id='cite-_Ref106273836'></span>[[#_Ref106273836|Fig. 11]]b). Además, estos frentes de onda de presión se pueden clasificar como: interfase positiva (el frente de onda de presión se mueve hacia el flujo en lámina libre), interfase negativa (el frente de onda de presión se mueve hacia el flujo en presión), o interfase con cambio de dirección (el frente de onda se mueve de una interfase positiva hacia una interfase negativa o viceversa).

<span id='_Ref106275164'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image12.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 12. Método de la ranura de Preissmann: (a) flujo en lámina libre; (b) flujo en presión; (c) flujo en presión negativa.'''</span>


En el módulo de drenaje urbano en Iber, para el cálculo de los flujos a presión, tanto en una dimensión como en dos, se ha implementado el método de la ranura de Preissmann [21]. El método de la ranura de Preissmann consiste en implementar una ranura hipotética de ancho <math display="inline">{T}_{s}</math> en la clave de los conductos cerrados (<span id='cite-_Ref106275164'></span>[[#_Ref106275164|Fig. 12]]b), proporcionando la presión hidrostática extra mientras que se permite el uso de las ecuaciones de flujo en lámina libre para ambos tipos de flujo (flujo en lámina libre y flujo en presión). El ancho de la ranura se calcula bajo la consideración de que la celeridad de la onda de gravedad <math display="inline">c</math> es igual a la celeridad de la onda de presión <math display="inline">a</math>. Esencialmente, esta aproximación explota la similitud de la ecuación de onda del flujo en lámina libre y flujo en presión. Así, a partir de la definición de la celeridad de la onda en lámina libre, tenemos:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>c=\sqrt{g\frac{A}{B}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(21)
|}

Despejando el ancho de la superficie libre del agua de la expresión anterior, que al emplear el método de la ranura de Preissmann se simboliza con <math display="inline">{T}_{s}</math> (<span id='cite-_Ref106275164'></span>[[#_Ref106275164|Fig. 12]]b), utilizando el área del colector a tubo lleno denominada área máxima <math display="inline">{A}_{max}</math> y empleando la celeridad de la onda de presión <math display="inline">a</math>, se obtiene:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{T}_{s}=\frac{g{A}_{max}}{{a}^{2}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(22)
|}

Cuando en un colector ocurre flujo en presión, el calado calculado es la carga de presión. Así, con el método de la ranura de Preissmann clásico, el área en presión es <math display="inline">{A}_{p}=A-</math><math>{A}_{max}</math> y la carga debida a la presión <math display="inline">{h}_{p}</math> se obtiene como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{h}_{p}=\frac{{A}_{p}}{{T}_{s}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(23)
|}

Para mitigar los problemas relacionados con las inestabilidades numéricas debidas a los esquemas numéricos cuando el flujo cambia rápidamente de flujo en lámina libre a flujo en presión, se plantea el uso de una transición geométrica gradual entre el conducto y el ancho de la ranura de Preissmann. Barnett [22] sugiere una transición geométrica gradual de tipo exponencial, la cual ha sido modifican ligeramente, resultando:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>w={T}_{s}{e}^{k\left( 1,5{h}_{max}-h\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(24)
|}

donde <math display="inline">h</math> es el calado o carga de presión, <math display="inline">{h}_{max}</math> es el calado máximo, <math display="inline">w</math> el ancho de la ranura en la zona de la transición geométrica, <math display="inline">k</math> un parámetro que determina la forma de la transición; valores más pequeños corresponden a transiciones más elongadas. En Iber se emplea la transición geométrica gradual a partir de una carga de presión de <math display="inline">h=</math><math>0,98{h}_{max}</math> para conductos con sección transversal convergente en su parte superior y de <math display="inline">h=</math><math>{h}_{max}</math> para otro tipo de secciones transversales, ambos hasta un calado de <math display="inline">\, h=</math><math>1,5{h}_{max}</math>.

La principal limitación de este método es que asume flujo a superficie libre en todo el sistema. Por consiguiente, si la carga piezométrica cae por debajo de la coronación del conducto, se desarrolla flujo en lámina libre y, por tanto, sería incapaz de simular flujo en presión negativa. Con el fin de simular flujo en presión subatmosférica, Kerger et al. [23] ampliaron este concepto y desarrollaron la ranura de Preissmann negativa que consiste en extender la ranura por debajo de la coronación del conducto (<span id='cite-_Ref106275164'></span>[[#_Ref106275164|Fig. 12]]c). De esta forma, cuando el flujo está en presión, para cada valor de calado menor al calado máximo existen dos valores de área: uno para superficie libre y otro para flujo en presión negativa.

Cabe decir que el método de la ranura de Preissmann podría ser numéricamente inestable cuando el ancho de la ranura es pequeña, necesaria para representar grandes niveles de agua o una celeridad de onda en presión importante. En Iber se considera correctamente la calidad de la onda de presión en la reducción del incremento de tiempo por la condición de Courant.

<span id='_Toc108335448'></span>

====3.1.6 Uniones de colectores y dominios 1D y 2D====

En una red de colectores puede haber zonas adecuadas para simular en 1D, pero otras que conviene simular en 2D, como pueden ser depósitos de retención, cámaras de disipación de energía, uniones complejas en colectores, entre otras. Por ello se ha implementado una formulación para la conexión entre dominios calculados en 1D y dominios calculados en 2D.

Para uniones entre dominios 1D y 2D se utiliza el método completamente conservativo presentado por Bladé et al. [16]. En él, una unión <math display="inline">k</math> se puede discretizar en <math display="inline">i</math>'' ''elementos o volúmenes finitos en 2D y <math display="inline">j</math> elementos frontera o ficticios; cada uno de los <math display="inline">i</math> colectores ( <math display="inline">p</math> colectores de entrada y <math display="inline">q</math> colectores de salida) se puede descomponer en <math display="inline">n</math>''-2'' celdas o volúmenes finitos en 1D y 2 celdas frontera o ficticios (volúmenes finitos 1 y <math display="inline">n</math>). La conexión entre los dominios en 1D (conductos) y en 2D (unión) se realiza a través de los volúmenes finitos ubicados en la frontera de cada dominio, para lo cual existe un conjunto de elementos (volúmenes finitos en 2D que representan la unión) conectados a una celda (volumen finito en 1D que representa un tramo de colector, cuyo centro alberga los valores de las variables hidráulicas y las propiedades geométricas de la sección transversal del conducto).

<span id='_Toc108335449'></span>

====3.1.7 Uniones de colectores====

Las uniones son los elementos que conectan los tramos de ríos, las calles en una zona urbana, los canales en una red de riego o los colectores en una red de alcantarillado. En una red de alcantarillado, un pozo de registro al inicio de un colector se puede entender como una unión con un único colector de salida. Este tipo de unión, además, puede recibir un volumen de agua del escurrimiento superficial. Uniones comunes son aquellas formadas por dos, tres o cuatro conductos, pero también existen uniones con cinco o más conductos.

La unión de colectores presenta una problemática específica, ya que en una unión pueden concurrir colectores con alineaciones verticales y horizontales muy diferentes, producirse cambios de sección y pendiente, así como notables variaciones de caudal entre colectores de entrada y salida [24]. Estas singularidades geométricas del colector y de las uniones, pueden provocar que en los colectores se presenten múltiples combinaciones de los diferentes tipos de flujo.

Para el cálculo de las uniones de colectores simples (es decir, donde no se requiere la caracterización detallada del flujo en 2D), la estrategia implementada en Iber es considerar la unión como un único elemento 2D. Ello se hace internamente, sin información al usuario, y la única opción que existe es si se considera continuidad de flujo a través de la unión, o no. En el caso de considerar continuidad de flujo, que sería para el caso donde los colectores de entrada están bien alineados, los caudales de entrada y salida están relacionados por la conservación de la cantidad de movimiento. En caso contrario, que correspondería a la falta de alineación entre colectores de entrada y salida, la presencia de obstáculos, o directamente de disipadores de energía, internamente se considera que la velocidad en el colector es nula.

<span id='_Toc108335450'></span>

====3.1.8 Drenaje dual====

El drenaje urbano integra dos componentes distintos: uno superficial, o sistema de drenaje mayor, compuesto por calles, cunetas, y canales naturales o artificiales; y un sistema subsuperficial, o sistema de drenaje menor, compuesto por colectores, uniones y estructuras de control como compuertas, vertederos, etc. (<span id='cite-_Ref106275985'></span>[[#_Ref106275985|Fig. 13]]). El sistema de drenaje menor habitualmente se diseña para funcionar en lámina libre para un caudal asociado a un periodo de retorno, relacionado con los eventos más frecuentes. Cuando se presenta un evento de lluvia con periodo de retorno superior al de diseño, los colectores (sistema de drenaje menor) pueden entrar en carga, incluso el agua puede llegar a salir a la superficie urbana y ser trasportada por el sistema de drenaje mayor (calles). La conexión entre los dos sistemas es bidireccional, es decir, el flujo puede ir del sistema de drenaje mayor al sistema de drenaje menor, y viceversa.

<span id='_Ref106275985'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
[[File:Draft_Sanz-Ramos_968687760_2436_Fig13.png|100%]]</div>

<span style="text-align: center; font-size: 75%;">'''Fig. 13. Elementos y procesos hidráulicos en el drenaje urbano ante un evento de lluvia.'''</span>


Los sumideros se diseñan para recolectar el agua de lluvia que circula sobre las calles, así como de otras áreas urbanas, y conducirla hacia los colectores, o lo que es lo mismo, del sistema de drenaje mayor al sistema de drenaje menor. En ocasiones puede existir conducción del flujo del sistema de drenaje menor al sistema de drenaje mayor, cuando la carga piezométrica en el colector es mayor a la superficie libre del agua en la calle en el punto donde se ubica dicho elemento. Los elementos que funcionan como fuentes son, por lo general, las tapas de los pozos de registro, pero los sumideros también pueden desempeñar dicha función.

Según lo descrito anteriormente, hay una interacción entre el sistema de drenaje mayor y el sistema de drenaje menor, que da lugar a un intercambio de flujo en ambas direcciones a través de las fuentes y los sumideros, el cual debe ser estimado adecuadamente para modelar el drenaje urbano dual convenientemente. De esta forma, en los puntos de localización de las fuentes y sumideros el proceso de interacción de flujo da lugar a tres casos:

:1) La carga piezométrica en el colector es menor que el nivel del terreno de la superficie urbana, por lo que ésta no interviene sobre el flujo de intercambio y el flujo captado del sistema de drenaje mayor entra al sistema de drenaje menor libremente (<span id='cite-_Ref106275995'></span>[[#_Ref106275995|Fig. 14]]a).
:2) La carga piezométrica del colector se encuentra entre el nivel del terreno de la superficie urbana y la superficie libre del agua en la calle, por lo que el flujo aún entra al sistema de drenaje menor, pero bajo la influencia de la carga hidráulica del colector (<span id='cite-_Ref106275995'></span>[[#_Ref106275995|Fig. 14]]b).
:3) La carga hidráulica en el colector es mayor a la superficie libre del agua en la calle, por lo tanto, el flujo escapa del sistema de drenaje menor hacia el sistema de drenaje mayor (<span id='cite-_Ref106275995'></span>[[#_Ref106275995|Fig. 14]]c).

<span id='_Ref106275995'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image14.png|90%]] </div>

<span id='_Ref343694560'></span><span style="text-align: center; font-size: 75%;">'''Fig. 14. Esquema de flujo en los puntos de interacción (fuente/sumidero): (a) ‘sumidero libre’; (b) ‘sumidero anegado’; (c) ‘fuente’.'''</span>


La capacidad hidráulica o de intercepción (cantidad de agua interceptada por un sumidero bajo ciertas condiciones) gobierna el caudal de intercambio (caudal interceptado por los sumideros) en los casos a) y b). El caudal interceptado es función de la posición, tamaño y tipo de reja, condiciones de flujo, rugosidad de la superficie, geometría de la calzada, y el fenómeno de colmatación [25,26]. Esta definición, válida para los casos a) y b), asume que la funcionalidad de los sumideros no es afectada por las condiciones existentes en los pozos de registro y en los colectores.

Para estimar los caudales de intercambio, existen formulaciones semiempíricas tipo orificio o vertedero y formulaciones experimentales basadas en la definición de la eficiencia hidráulica, las características del flujo y la geometría de los sumideros. La primera formulación se puede utilizar en los casos a), b) y c), mientras la segunda, solamente es aplicable al caso a).

'''''Formulación orificio/vertedero'''''

Este tipo de formulación considera que el flujo de intercambio es función del calado de aproximación. Así, el caudal de intercambio para el caso a) puede ser estimado mediante una formulación tipo vertedero para calados pequeños y una formulación tipo orificio para calados grandes. Como la transición de los regímenes de flujo de vertedero a orificio y viceversa no es muy clara, se recomienda utilizar la formulación que estime el caudal mínimo. En el caso b), el flujo puede ser caracterizado con una formulación tipo orificio o tipo compuerta. Finalmente, para el caso c) una ecuación tipo orificio es la más apropiada.

Las fórmulas tradicionales para estimar el caudal a través de un orificio <math display="inline">{Q}_{o}</math> y a través de un vertedero <math display="inline">{Q}_{v}</math> son:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{Q}_{o}={C}_{d,o}{A}_{int}\sqrt{2gH}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(25)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{Q}_{v}={C}_{d,v}{L}_{int}{H}^{3/2}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(26)
|}

donde <math display="inline">{A}_{int}</math> es el área de captación, <math display="inline">{L}_{int}</math>la longitud de captación, <math display="inline">H</math> la carga hidráulica, <math display="inline">{C}_{d,o}</math> el coeficiente de descarga del orificio y <math display="inline">{C}_{d,v}</math> el coeficiente de descarga del vertedero. Valores entre 0,60 y 0,80 son recomendados para el coeficiente de descarga del orificio y valores entre 1,70 y 3,20 para el coeficiente de descarga del vertedero.

Uno de los inconvenientes del cálculo del caudal a través de un vertedero consiste en determinar la longitud de intercepción correcta a emplea. En cambio, una formulación tipo orificio emplea el área <math display="inline">{A}_{int}</math> (corresponde al área de huecos). Cuando el intercambio de flujo es del sistema de drenaje mayor al sistema de drenaje menor y el calado cubre totalmente el elemento de captación, esta longitud se caracteriza con el perímetro total del elemento. En caso contrario, esta longitud será menor al perímetro total del elemento de captación.

'''''Metodología UPC'''''

Esta metodología permite estimar el caudal interceptado por un sumidero a partir su eficiencia hidráulica <math display="inline">E</math>, definida como la proporción o porcentaje entre el caudal interceptado <math display="inline">{Q}_{int}</math> y el caudal total de aproximación al sumidero <math display="inline">Q</math>.

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>E=\frac{{Q}_{int}}{Q}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(27)
|}

Los resultados de los ensayos experimentales, realizados en el laboratorio de hidráulica del Departament d’Enginyeria Hidràulica, Marítima i Ambiental (DEHMA) de la Escola Tècnica Superior d’Enginyers de Camins, Canals i Ports (ESTCCPB) de la Universitat Politècnica de Catalunya (UPC), para analizar el comportamiento hidráulico y capacidad de captación de diferentes tipos de rejas y macro-rejas [26], determinaron una ecuación de tipo potencial que relaciona la eficiencia de captación de rejas y macro-rejas con el calado y caudal de aproximación a la reja y de dos parámetros característicos de la geometría las mismas, definida como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>E=A{\left( \frac{Q}{h}\right) }^{-B}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref106282459'></span>(28)
|}

donde,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>A=\frac{1,988{A}_{g}^{0,403}}{{p}^{0,190}{\left( {n}_{t}+1\right) }^{0,088}{\left( {n}_{l}+1\right) }^{0,012}{\left( {n}_{d}+1\right) }^{0,082}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(29)
|}

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>B=1,346\frac{{L}^{0,179}}{{W}^{0,394}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(30)
|}

donde <math display="inline">h</math> es el calado de aproximación, <math display="inline">A</math> y <math display="inline">B</math> los parámetros geométricos característicos de cada reja. Estos parámetros se obtienen a partir del área mínima que engloba los hueco <math display="inline">\, {A}_{g}</math>, el porcentaje entre el área de huecos <math display="inline">\, {A}_{H}</math> y el área mínima que engloba los huecos <math display="inline">\, p=</math><math>100\left( {A}_{H}/{A}_{g}\right)</math> , el número de barras transversales <math display="inline">{n}_{t}</math>, el número de barras longitudinales <math display="inline">{n}_{l}</math>, el número de barras diagonales <math display="inline">{n}_{d}</math>, la longitud de la reja <math display="inline">L</math> y ancho de la reja <math display="inline">W</math>.

Es importante enfatizar que esta ecuación se obtuvo experimentalmente para el caudal que circula por un ancho de calzada de 3 m, por lo que se tuvieron las precauciones necesarias al momento de implementarla en el modelo numérico. Para anchos de calzada distintos, la Ecuación <span id='cite-_Ref106282459'></span>[[#_Ref106282459|(28)]] puede generalizarse considerando una distribución uniforme de la velocidad [26], y así se hace en Iber [27].

De forma similar, para el caso de rejas transversales continuas, Gómez y Russo [28] determinaron una ley que relaciona la eficiencia hidráulica de captación de estas estructuras con el número de Froude <math display="inline">\, {F}_{r}</math>, calado de aproximación y dos coeficientes característicos de la geometría <math display="inline">\, \alpha</math>  y <math display="inline">\beta</math> , resultando en:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>E=\alpha \left[ \, {F}_{r}{\left( \frac{h}{L}\right) }^{0,812}\right] +\, \beta</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(31)
|}

En este caso, la eficiencia hidráulica de captación de rejas continuas puede ser expresada en términos de una eficiencia unitaria, definida como el cociente entre el caudal interceptado y el caudal de aproximación, ambos por unidad de ancho.

Una vez determinada la eficacia con cualquiera de las dos fórmulas anteriores, el caudal interceptado <math display="inline">\, {Q}_{int}</math> se puede determinar como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{Q}_{int}=EQ</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(32)
|}

Esta metodología tiene la ventaja de poder extrapolar la formulación a otras rejas con características similares sin necesidad de ser ensayadas [26].

<span id='_Toc108335451'></span>

===3.2 Interfaz de usuario===

Como se ha indicado, para el funcionamiento del módulo de drenaje urbano se sigue la misma estrategia que para los módulos existentes: el cálculo se realiza con una malla de cálculo que tiene una serie de propiedades y condiciones, y dónde cada volumen finito corresponde a un elemento de la malla. La malla se suele generar a partir de una geometría previa, trasfiriendo posteriormente las propiedades asignadas en ella a la malla de cálculo. Es posible, sin embargo, asignar propiedades directamente a la malla, o modificar las que ya tiene.

A continuación, se describen los diferentes menús y ventanas desde las que se pueden activar los procesos descritos en las secciones anteriores y definir los correspondientes parámetros. El módulo se activa desde el menú '''Herramientas Iber >> Plug-ins >> Drenaje urbano'''. Esto permite al/la usuario/a utilizar e implementar únicamente variables relacionadas con los módulos que interaccionan con el drenaje urbano dual.

<span id='_Toc54793165'></span><span id='_Toc90987611'></span><span id='_Toc108335452'></span>

====3.2.1 Formulaciones, opciones y parámetros generales del modelo====

Las opciones de cálculo se especifican en '''Datos >> Datos del problema >> Drenaje urbano'''. Una vez activado el módulo de drenaje urbano, automáticamente aparece la pestaña ‘Drenaje Urbano’ en la ventana de Datos del Problema (<span id='cite-_Ref106295871'></span>[[#_Ref106295871|Fig. 15]]), desde donde se deben definir el tipo de flujo (superficial, subterráneo o dual) y el tipo de método de intercambio entre el flujo superficial y el flujo subterráneo, ya sea según el método ‘UPC’ [26] o como un orificio/vertedero.

<span id='_Ref106295871'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image15.png|500px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 15. Datos generales del módulo de drenaje urbano.'''</span>


<span id='_Toc108335453'></span>

====3.2.2 Red de drenaje====

La red de drenaje puede estar compuesta por uniones, conductos (que pueden unirse mediante uniones) y zonas donde el flujo subterráneo puede fluir o almacenarse (superficie 2D). Los elementos que conectan el flujo superficial con el flujo en la red de drenaje subterránea son las rejas (o ''inlets'' en inglés), y son los que permiten el intercambio de caudales entre la parte superficial (p.ej. calles) y la parte subterránea (p.ej. la red de drenaje).

Esta red puede construirse de forma manual, empleando las herramientas disponibles en Iber para la creación de geometrías (líneas, superficies), o bien importarse desde archivos externos, pudiendo ser o no necesaria la implementación de ciertas condiciones a posteriori.

'''''Definición manual de la red'''''

Para la definición completa de la red es necesario especificar qué conductos la forman (líneas), las uniones entre conductos (puntos) y, en su caso, qué zonas existen para ser calculadas en dos dimensiones (superficies). A la definición de la red de drenaje se accede mediante el menú '''Datos >> Drenaje urbano >> Red de drenaje'''. Desde aquí es posible definir las características de las uniones, los conductos y las superficies de la red de drenaje.

Las '''Uniones''' se definen mediante sus características geométricas (<span id='cite-_Ref106296101'></span>[[#_Ref106296101|Fig. 16]]a): la cota de fondo, su diámetro (o diámetro equivalente si no son circulares), el calado inicial y su cota superior (para láminas de agua superiores a la cota superior la unión funciona a presión). Además, es posible indicar si se desea permitir la continuidad del flujo entre conducciones de entrada y salida (opción ''Continuidad del flujo en la unión'') o no. En caso de no existir continuidad, se supone que la velocidad del agua en la unión es nula, o lo que es lo mismo, que no hay continuidad entre la línea de energía de los conductos entrantes y salientes. En las uniones es posible asignar, también, un caudal de entrada.

En cuanto a las '''Conducciones''', pueden ser de tipo circular o rectangular (<span id='cite-_Ref106296101'></span>[[#_Ref106296101|Fig. 16]]b), y deben asignarse a líneas o elementos lineales. Los conductos pueden trabajar en presión, por lo que es necesario definir la celeridad de la onda (en metros por segundo).

<span id='_Ref106296101'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image16.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 16. Definición e implementación de la red de drenaje: (a) uniones; (b) conducciones; (c) zonas de la subsuperficie.'''</span>


Finalmente, está la posibilidad de que existan zonas en la red de drenaje que tengan una geometría más compleja y requieran un cálculo en 2D. Para ello se dispone de las '''Áreas de la subsuperficie''' (<span id='cite-_Ref106296101'></span>[[#_Ref106296101|Fig. 16]]c). Aquí tan solo se tiene que indicar qué superficie corresponde a este tipo de estructura, y el programa automáticamente calculará la hidrodinámica en 2D teniendo en cuenta las entradas/salidas de flujo que haya en sus contornos (más adelante se indica cómo conectar dichas superficies con el resto de la red de drenaje).

'''''Importación de la red'''''

La importación, no solo de la geometría de la red, sino también de las características de la red (uniones, conductos, áreas no conectadas con la superficie) puede realizarse desde el menú '''Datos >> Drenaje urbano >> Importar SPH'''. Para ello es necesario disponer de un archivo ESRI Shapefile (''*.shp'') que contenga la definición geométrica de la red, ya sea en 2D o 3D.

<span id='_Ref106296380'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image17.png|450px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 17. Opción ‘Crear condiciones’ para que Iber lea la base de datos del *.shp (archivo *.dbf asociado) y asigne las condiciones a la geometría.'''</span>


Si el archivo ''*.shp'' va acompañado de una base de datos (''*.dbf''), se podrá emplear para asignar las condiciones de la red (tipo de unión, dimensiones del conducto, rugosidad, etc.) también de forma automática (opción ‘Crear condiciones’ activada, <span id='cite-_Ref106296380'></span>[[#_Ref106296380|Fig. 17]]).

Este archivo debe tener una estructura concreta para que la importación de todas las condiciones sea conforme. En caso contrario, tan solo se importará la geometría. La estructura de la base de datos para la asignación automática de las condiciones de la red en Iber para elementos tipo '''Uniones''' se detalla en la <span id='cite-_Ref106296445'></span>[[#_Ref106296445|Tabla 1]]:

<span id='_Ref106296445'></span>
<span style="text-align: center; font-size: 75%;">'''Tabla 1. Estructura de la base de datos (*.dbf) que acompaña al archivo de geometría (*.shp) de las Uniones.'''</span>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 1'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 2'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 3'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 4'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 5'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 6'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 7'''
|-
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Elevación superior
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Elevación del fondo
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Calado
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Diámetro
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Continuidad flujo
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Entrada/salida flujo
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Tabla
|}


Las Columnas 5 y 6 tienen que ser un número entero, pudiendo tomar 0, en caso de que no se desee habilitar dicha opción, o 1, en caso que se desee habilitar dicha opción. Si la Columna 6 toma valor 1, se espera que en la Columna 7 se defina la tabla de intercambio de entrada/salida de flujo en la unión. Esta tabla se tiene que definir de la siguiente manera:

:* Tipo: texto
:* Campo como sigue: #N# XX A1 A2 B1 B2 C1 C2 … n1 n2

El valor del campo de la Columna 7, separado mediante espacios, queda definido en primer lugar con “#N#”, seguido del número de elementos que tiene la tabla (“XX”) y luego los valores de la tabla agrupados por parejas, es decir, tiempo-caudal. Así, por ejemplo, un hidrograma de entrada de forma triangular simétrico con un caudal punta de 100 L/s a la 1h, deberá introducirse en la Columna 7 como: #N# 6 0 0 3600 0.1 7200 0.

En este sentido, la estructura de la base de datos para la asignación automática de las condiciones de la red en Iber para elementos tipo '''Conducciones''' es la que se detalla en la <span id='cite-_Ref106296673'></span>[[#_Ref106296673|Tabla 2]]:

<span id='_Ref106296673'></span>
<span style="text-align: center; font-size: 75%;">'''Tabla 2. Estructura de la base de datos (*.dbf) que acompaña al archivo de geometría (*.shp) de las Conducciones.'''</span>

{| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 1'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 2'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 3'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 4'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 5'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 6'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Col. 7'''
|-
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Núm. conducción
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Tipo
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Diámetro
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Ancho
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Alto
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Manning
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Celeridad
|}


La Columna 2 define el tipo de conducción, por lo que debe ser un campo tipo Texto cuyo valor debe ser “Rectangular” o “Circular”. En caso que la Columna 2 sea tipo “Circular”, la Columna 3 debe tomar el valor del diámetro del conducto, mientras que las Columnas 4 y 5 pueden tomar cualquier valor, pero deben existir. En cambio, en caso que la Columna 2 sea tipo “Rectangular”, la Columna 3 podrá tomar cualquier valor, pero deberá existir, mientras que las Columnas 4 y 5 deben tomar el valor del ancho y el alto del conducto, respectivamente.

Una vez importada la geometría del modelo, es posible modificar los parámetros mediante los menús comentados para la definición manual de la misma.

<span id='_Toc108335454'></span>

====3.2.3 Sumideros====

'''''Definición manual de los sumideros'''''

La definición/modificación de los sumideros se realiza desde el menú '''Datos >> Drenaje urbano >> Rejas'''. La <span id='cite-_Ref106296970'></span>[[#_Ref106296970|Fig. 18]]a muestra la ventana para la implementación de rejas, que pueden ser tipo longitudinal, transversal o como una reja circular (similar a un pozo de registro). La <span id='cite-_Ref106296970'></span>[[#_Ref106296970|Fig. 18]]b muestra, a modo de ejemplo, cómo se representaría la reja (resaltada en amarillo) y su conexión con la red de drenaje.

Para cada reja o sumidero lo primero que hay que hacer es crear-la con el botón correspondiente y luego introducir sus características. Los parámetros que hay que introducir son:

:* Posición en superficie: Coordenadas
:* Tipo: puede ser “''To network''” (a la red) o “''Sink''” (sumidero sin que el agua se incorpore a ningún colector)
:* Posición en la red: Coordenadas del punto de conexión a la red
:* Tipo: Transversal o longitudinal
:* A: Parámetro A del método UPC
:* B: Parámetro B del método UPC
:* Longitud
:* Anchura
:* Área efectiva (no obstruida)
:* Coeficiente de desagüe de vertedero para la metodología orificio/vertedero
:* Coeficiente de desagüe de orificio para la metodología orificio/vertedero

En el caso de que exista flujo entrando a la reja, además del que intercepta de la superficie, es posible incorporarlo a través de un hidrograma.

<span id='_Ref106296970'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image18.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 18. Definición e implementación de rejas: (a) ventana para la introducción de datos; (b) representación de la reja y su conexión con la red de drenaje.'''</span>


'''''Importación de sumideros'''''

También, desde el menú '''Datos >> Drenaje urbano >> Importar archivo de rejas''', se posibilita la importación de la localización y propiedades de las rejas. Se trata de un archivo de texto que debe contener el identificador de la reja, su ubicación, el tipo (longitudinal, transversal o pozo de registro), sus características geométricas y propiedades hidráulicas. También es posible indicar un hidrograma de entrada/salida en ese elemento.

La <span id='cite-_Ref106297414'></span>[[#_Ref106297414|Fig. 19]] presenta la estructura del archivo requerido. Cada fila corresponde a una reja con los siguientes parámetros:

:* Número de reja
:* 1 si es visible (se representa en ventana), 0 en caso contrario
:* Coordenadas (x,y,z) en superficie
:* 1 si es un sumidero, 2 si hay conexión a la red
:* Coordenadas (x,y,z) de la conexión al colector
:* 1 si es longitudinal, 2 si es transversal 3 si es pozo de registro
:* A
:* B
:* Longitud
:* Anchura
:* Porcentaje de área efectiva
:* Cd Vertedero
:* Cd Orificio
:* 0 si no hay hidrograma asociado, 1 en caso que haya
:* Nombre del sumidero (sin espacios)
:* Número de valores del hidrograma (1 en caso de que no haya hidrograma)
:* Tantas filas como valores del hidrograma, con parejas de datos tiempo-caudal

<span id='_Ref106297414'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image19.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 19. Estructura del archivo de importación automática de rejas.'''</span>


El formato del archivo de importación corresponde con el formato del archivo “''Iber_Sink&Sources.dat''”, que se crea dentro de la carpeta del proyecto al ejecutar un cálculo con sumideros introducidos manualmente. De esta manera, leyendo este archivo de un proyecto ya creado, se pueden importar los sumideros a un proyecto nuevo.

<span id='_Toc108335455'></span>

====3.2.4 Condiciones de contorno y conexiones====

La red de drenaje puede presentar diferentes tipos de condiciones de contorno, así como conexiones directas con la superficie (además de las ''Rejas''). Para ello se han generado nuevas ventanas donde se pueden implementar estas condiciones.

<span id='_Ref106297881'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image20.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 20. Ventana para la introducción de datos de las condiciones de contorno de entrada de la red: (a) régimen supercrítico; (b) crítico/subcrítico.'''</span>


Por un lado, se encuentran las condiciones de contorno de la propia red. Su implementación se realiza desde el menú '''Datos >> Hidrodinámica >> Condiciones de contorno de la red'''. Se pueden indicar las condiciones de contorno de entrada, que pueden ser en régimen supercrítico (<span id='cite-_Ref106297881'></span>[[#_Ref106297881|Fig. 20]]a) o en régimen crítico/subcrítico (<span id='cite-_Ref106297881'></span>[[#_Ref106297881|Fig. 20]]b). En cada caso se muestran el tipo de condiciones que se pueden implementar, ya sea una tabla Tiempo-Caudal-Elevación/Calado o una tabla Tiempo-Caudal (hidrograma). Para las condiciones de salida, se sigue una filosofía análoga, cambiando tan solo el tipo de salida según el régimen hidráulico. Estas condiciones deben imponerse en un contorno de la red (puntos extremos de los colectores).

Por otro lado, se encuentran las condiciones de contorno para elementos de tipo superficie, que son análogos a los ya existentes en Iber para la modelización de flujo en superficie. La ventana es accesible desde el menú Datos >> Hidrodinámica >> Condiciones de contorno, y ha sido adaptada para contener opciones de conexión entre elementos 1D (conducciones de la red, <span id='cite-_Ref106298399'></span>[[#_Ref106298399|Fig. 21]]a) y elementos 2D (''Zonas de la subsuperfice'' o zonas de la superficie, <span id='cite-_Ref106298399'></span>[[#_Ref106298399|Fig. 21]]b).

La existencia de condiciones de contorno, tanto en superficie como en la red de drenaje, requiere conectar ambos tipos de flujo. Por ejemplo, es posible que el flujo que sale por una condición 2D sea el flujo que entra en una condición 1D, y viceversa. Por ese motivo existen las opciones “Conexión 1D” y “Conexión 2D” en las condiciones de contorno 2D y 1D, respectivamente.

<span id='_Ref106298399'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
[[Image:Draft_Sanz-Ramos_968687760-image21.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 21. Condiciones de contorno de la red: (a) opciones para la conexión entre condiciones salida/entrada 1D; (b) condiciones entrada/salida 2D.'''</span>

En las condiciones de contorno de la red de drenaje (1D), se permite la conexión con entidades 2D. En ese caso, se solicita el punto de la red que va a ser considerado como una entrada/salida en una entidad 2D del modelo. Del mismo modo, en una condición de contorno 2D se permite la conexión con la red de drenaje 1D, es necesario indicar el punto de la red en la que va a entrar/salid flujo. Para facilitar su selección, se permite la selección del punto de la red a través de la interfaz.

<span id='_Toc108335456'></span>
<span id='_Ref106262222'></span>

==4 Módulo de hidrología distribuido==

<span id='_Toc54965419'></span><span id='_Toc90987620'></span><span id='_Toc108335457'></span>

===4.1 Fundamentos teóricos===

El módulo de procesos hidrológicos de Iber permite calcular la transformación lluvia-escorrentía en toda una cuenca, haciendo posible su utilización como modelo hidrológico distribuido en cuencas de pequeño y mediano tamaño. Es completamente bidimensional y distribuido en el espacio, y trabaja sobre la misma malla de cálculo que el módulo hidrodinámico.

Este módulo se basa en las ecuaciones hidrodinámicas de aguas poco profundas bidimensionales, las cuales se utilizan para resolver el movimiento de la escorrentía superficial sobre el terreno. Se le han añadido una serie de nuevos procesos como son, el aporte de agua de lluvia, la infiltración de agua en el terreno, el flujo subsuperficial, la evapotranspiración y la percolación.

<span id='_Toc108335458'></span>

====4.1.1 Procesos considerados====

La base del módulo hidrológico es el módulo hidrodinámico, al que se le han añadido los procesos de aporte de agua por precipitación y la infiltración del agua en el terreno. De esta forma, el exceso de precipitación se mueve por el terreno y por los cauces de acuerdo con las ecuaciones de aguas someras bidimensionales. El agua infiltrada pasa a formar parte del subsuelo, donde se realiza un balance de masa de agua distribuido en el espacio a nivel de elemento de la malla de cálculo. En dicho balance de agua en el terreno se consideran los procesos de infiltración, evapotranspiración, percolación y flujo en el subsuelo. Los tres primeros implican un movimiento de agua en la columna vertical, mientras que el flujo en el subsuelo se considera bidimensional y paralelo al terreno. Todos estos procesos, representados de forma esquemática en la <span id='cite-_Ref106365993'></span>[[#_Ref106365993|Fig. 22]], pueden ser activados o desactivados por el usuario.

<span id='_Ref106365993'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image22.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 22. Esquema de los procesos considerados en el módulo hidrológico.'''</span>


<span id='_Toc108335459'></span>

====4.1.2 Escorrentía superficial====

El cálculo de la escorrentía superficial, y su propagación a lo largo de la cuenca hidrográfica (incluyendo laderas y cauces de agua), se realiza mediante la resolución de las ecuaciones de aguas poco profundas bidimensionales, incluyendo términos que representan las componentes de precipitación e infiltración (Ecuación <span id='cite-_Ref106269045'></span>[[#_Ref106269045|(1]])).

Este sistema de ecuaciones es el mismo utilizado en Iber para calcular la hidrodinámica en un curso fluvial. La única diferencia es que se han incluido 2 nuevos esquemas numéricos para su resolución, que suelen dar resultados numéricamente más estables en aplicaciones hidrológicas. Estos modelos son los siguientes:

:* DHD
:* DHD_basin

La descripción detallada del modelo DHD se presenta en el artículo Cea y Bladé [29]. En cuanto al modelo DHD_basin, es una variante del modelo DHD que se diferencia únicamente en el tratamiento de los frentes que separan las zonas secas de las mojadas (frente seco-mojado). En el esquema DHD, la cota de la lámina de agua en los frentes seco-mojado se calcula utilizando la ecuación (15) de Cea y Bladé [29]. En el esquema DHD_basin, la altura de la lámina en los frentes seco-mojado se calcula mediante una interpolación lineal de la cota de agua en el elemento mojado y la cota del terreno en el elemento seco.

Se recomienda la utilización de uno de estos dos esquemas (preferiblemente el DHD_basin) en aplicaciones que impliquen transformación lluvia-escorrentía en la totalidad de una cuenca hidrológica. Sin embargo, estos esquemas producen resultados menos precisos que los esquemas de Roe de orden 1 y orden 2 para el cálculo de flujo en ríos, canales o estuarios, por lo que no deben utilizarse en este tipo de aplicaciones.

<span id='_Toc108335460'></span>

====4.1.3 Precipitación====

La precipitación puede definirse a partir de series temporales de intensidad de precipitación observadas en pluviómetros de coordenadas conocidas (hietogramas), o a través de ficheros raster georreferenciados. Estos últimos se definen como campos espaciales de intensidad o profundidad de precipitación en diferentes instantes de tiempo. En ambos casos, la precipitación puede ser variable en tiempo y espacio.

En el caso de imponer la lluvia a partir de hietogramas, el usuario puede definir tantos pluviómetros como desee. Dichos pluviómetros pueden asignarse manualmente a diferentes zonas del modelo, con el fin de definir manualmente la distribución espacial de la lluvia. En el caso de que no se asignen manualmente, el programa realizará una interpolación espacial de la lluvia utilizando el método del vecino más cercano (''natural neighbour''), utilizando para ello las coordenadas espaciales de cada pluviómetro.

En el caso de definir la lluvia a partir de ficheros raster, el usuario debe suministrar una serie de ficheros georreferenciados (en el mismo sistema de referencia que el resto del modelo) en formato ASCII Grid de ArcInfo. Cada raster define la intensidad de precipitación o la profundidad de precipitación durante un intervalo de tiempo específico. Todos los ficheros raster deben tener el mismo número de columnas y filas, el mismo tamaño de pixel y el mismo origen. Los ficheros raster pueden tener cualquier nombre, pero todos ellos deben copiarse en una subcarpeta llamada ''Rain_rasters'', localizada dentro de la carpeta del proyecto (''proyecto.gid/Rain_rasters''). La carpeta ''Rain_rasters'' debe crearla el usuario.

<span id='_Ref106600646'></span><span id='_Ref107216854'></span><span id='_Toc108335461'></span>

====4.1.4 Infiltración====

Se incluyen los cuatro modelos de infiltración siguientes:

:* Lineal
:* Green-Ampt
:* Número de Curva (SCS)
:* Horton

Estas formulaciones son las más extendidas para el cálculo de la infiltración, brindando suficiente versatilidad al usuario a la hora de definir y parametrizar las pérdidas por infiltración. En todo caso, la infiltración de agua en un terreno real es un proceso complejo que puede, en ocasiones, no estar correctamente representado por ninguna de estas formulaciones. Asimismo, la determinación de los parámetros de estas formulaciones suele requerir su calibración mediante datos de campo.

El modelo lineal es el más sencillo, y permite introducir unas pérdidas iniciales seguidas de una infiltración potencial constante en el tiempo. Puede ser adecuado para eventos cortos en los cuales puede considerarse la infiltración más o menos constante en el tiempo. También puede utilizarse en casos en los que la capacidad de infiltración del suelo sea muy elevada inicialmente (superior a la precipitación), hasta que llega un momento en el que se satura el suelo y la infiltración disminuye drásticamente a un valor aproximadamente constante en el tiempo.

El modelo de Green-Ampt utiliza una formulación basada en parámetros físicos que pueden relacionarse con ciertas propiedades del terreno, como son la porosidad, la permeabilidad saturada, la succión o el contenido de humedad del suelo. Genera una infiltración potencial elevada inicialmente (dependiente de la succión y de la humedad inicial del suelo), que va decayendo de forma aproximadamente exponencial en el tiempo, hasta llegar a un valor mínimo de infiltración igual a la permeabilidad saturada del terreno.

El modelo del Número de Curva del SCS es una formulación ampliamente conocida y utilizada en múltiples modelos hidrológicos. A unas pérdidas iniciales le sigue una capacidad de infiltración variable, que disminuye de forma aproximadamente exponencial en el tiempo. Su principal ventaja es que depende de un único parámetro, cuyo valor puede ser relacionado con la pendiente del terreno, el tipo de suelo y los usos del suelo. Existe, además, numerosa bibliografía para estimar su valor.

El modelo de Horton es similar al de Green-Ampt en cuanto a la capacidad de infiltración generada, pero sus parámetros no están directamente relacionados con las características del terreno.

'''''Modelo de infiltración lineal'''''

En el modelo lineal, la tasa de infiltración potencial ( <math display="inline">{f}_{pot}</math>) se asume constante en el tiempo. El usuario introduce su valor en mm/h como un parámetro del modelo, que puede ser variable en espacio. La infiltración real ( <math display="inline">f</math>) en cada instante de tiempo es igual a la infiltración potencial siempre y cuando la intensidad de lluvia sea superior a <math display="inline">{f}_{pot}</math>, o si el calado sobre la superficie es suficiente para satisfacer su valor.

Además, el modelo lineal considera unas pérdidas iniciales ( <math display="inline">{I}_{a}</math>) definidas por el usuario en mm. Las pérdidas iniciales se satisfacen con la lluvia inicial o con el calado de agua sobre la superficie del terreno.

'''''Modelo de infiltración de Green-Ampt'''''

En el modelo de Green-Ampt [30] se asume que existe un frente de saturación que separa una capa de suelo saturado, situada inmediatamente bajo la superficie del terreno, de una capa inferior de suelo no saturado. A medida que el agua se infiltra en el terreno, el frente saturado desciende y el espesor de la capa saturada ( <math display="inline">b</math>) aumenta.

En la implementación del modelo de Green-Ampt en Iber, la tasa de infiltración potencial ( <math display="inline">{f}_{pot}</math>) y la infiltración real ( <math display="inline">f</math>) se calculan como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{f}_{pot}={k}_{s}\left( 1+\frac{\left( h+\psi \right) }{b}\right) \quad \quad f=</math><math>\mathrm{min}\,\left( {f}_{pot},\, \frac{h}{\Delta t}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(33)
|}

donde <math display="inline">{k}_{s}</math> es la permeabilidad saturada del suelo en la dirección vertical, <math display="inline">h</math> es el calado de agua sobre el terreno, <math display="inline">\psi</math>  es la succión en la capa de suelo no saturada y <math display="inline">b</math> es el espesor de la capa de suelo saturada. La succión ( <math display="inline">\psi</math> ) y la permeabilidad saturada ( <math display="inline">{k}_{s}</math>) son parámetros del modelo que debe introducir el usuario, mientras que el espesor de la capa de suelo saturada se actualiza a medida que se infiltra más agua en el terreno mediante las siguientes expresiones:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>b={b}_{0}+\frac{F}{\phi }\quad \quad \, F=\int_{0}^{t}f\, dt</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(34)
|}

donde <math display="inline">{b}_{0}</math> es el espesor inicial de la capa de suelo saturada, que depende del contenido de humedad inicial del suelo al comienzo de la simulación, y <math display="inline">\phi</math>  es la porosidad del suelo.

Dado que Iber es un modelo 2D en el que no se resuelve el flujo en dirección vertical, el espesor de la capa de suelo saturada ( <math display="inline">b</math>) se relaciona con el contenido de humedad promediado verticalmente en la capa de suelo ( <math display="inline">\theta</math> ), mediante la siguiente relación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{b}{{d}_{s}}=\frac{\theta }{\phi }\quad \quad {b}_{0}={\theta }_{0}\frac{{d}_{s}}{\phi }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(35)
|}

donde <math display="inline">{d}_{s}</math> es el espesor total de la capa de suelo, <math display="inline">\theta</math>  es el contenido de humedad promediado verticalmente en la capa de suelo. Es decir, el valor de <math display="inline">b</math> se obtiene agrupando todo el contenido de agua de la capa de suelo en una capa horizontal situada justo bajo el terreno (<span id='cite-_Ref106367812'></span>[[#_Ref106367812|Fig. 23]]). Los valores del espesor total de la capa de suelo ( <math display="inline">{d}_{s}</math>), su porosidad ( <math display="inline">\phi</math> ) y la humedad inicial ( <math display="inline">{\theta }_{0}</math>), son parámetros del modelo que debe introducir el usuario.

<span id='_Ref106367812'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image23.png|400px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 23. Esquema de la distribución de humedad en la dirección vertical en una capa de suelo.'''</span>


Cuando el espesor de la capa saturada alcanza el espesor total de la capa de suelo ( <math display="inline">b=</math><math>{d}_{s}</math>), ya no es posible infiltrar más agua. Esta situación se corresponde a la de un suelo totalmente saturado ( <math display="inline">\theta =</math><math>\phi</math> ).

El modelo de Green-Ampt permite, asimismo y de forma opcional, definir unas pérdidas iniciales ( <math display="inline">{I}_{a}</math>) definidas por el usuario en mm. Las pérdidas iniciales se satisfacen con la lluvia inicial o con el calado de agua sobre la superficie del terreno.

En la implementación del modelo de Green-Ampt se consideran, por tanto, un máximo de seis parámetros: <math display="inline">{k}_{s}</math>, <math display="inline">\psi</math> , <math display="inline">\phi</math> , <math display="inline">{d}_{s}</math>, <math display="inline">{\theta }_{0}</math>, <math display="inline">{I}_{a}</math>.

La principal ventaja del modelo de Green-Ampt es que se puede recuperar capacidad de infiltración realizando un balance de masa de agua en el suelo, incorporando los procesos de percolación hacia capas más profundas, de evapotranspiración y de flujo subsuperficial. En efecto, teniendo en cuenta la relación entre la humedad del suelo y el espesor de la capa saturada, la infiltración potencial puede escribirse como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{f}_{pot}={k}_{s}\left( 1+\frac{\left( h+\psi \right) }{\theta }\, \frac{\phi }{{d}_{s}}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(36)
|}

De esta forma, si la percolación y la evapotranspiración superan a la infiltración, la humedad del suelo disminuirá y, por tanto, la capacidad de infiltración aumentará. En el caso de que no se consideren los procesos de percolación o evapotranspiración, el contenido de humedad siempre aumentará y la capacidad de infiltración disminuirá de forma continua.

El modelo puede simplificarse reduciendo el número de parámetros. Si el usuario define las pérdidas iniciales y la humedad inicial igual a cero ( <math display="inline">{I}_{a}=</math><math>0</math> y <math display="inline">{\theta }_{0}=0</math>), y un espesor de suelo suficientemente elevado ( <math display="inline">{d}_{s}\rightarrow \infty</math> ), los parámetros del modelo se reducen a los tres utilizados en la implementación clásica del modelo de Green-Ampt: <math display="inline">{k}_{s}</math>, <math display="inline">\psi</math> , <math display="inline">\phi</math> , y la infiltración se calculará como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{f}_{pot}={k}_{s}\left( 1+\frac{\left( h+\psi \right) \, \phi }{F}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(37)
|}

'''''Modelo de infiltración del SCS'''''

En el método del número de curva del SCS [31–33], se define la capacidad máxima de retención del suelo ( <math display="inline">S</math>) como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>S=\frac{25400}{CN}-254</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(38)
|}

donde <math display="inline">CN</math> es el Número de Curva, a definir por el usuario en función del tipo de suelo, usos del suelo y pendientes del terreno. La abstracción inicial ( <math display="inline">{I}_{a}</math>) se define como un porcentaje de la capacidad máxima de retención, como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{a}=\alpha \, S</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(39)
|}

En la implementación clásica del método del número de curva, la abstracción inicial es un 20% de la capacidad máxima de retención ( <math display="inline">{I}_{a}=</math><math>0,2\, S</math>). Sin embargo, diversos estudios han encontrado que en muchos casos las pérdidas iniciales son inferiores a dicho valor, por lo que en la implementación en Iber se introduce el valor de <math display="inline">\alpha</math>  como un parámetro a definir por el usuario.

A partir de los parámetros anteriores, la precipitación neta acumulada en cada instante se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{P}_{n}=\frac{{\left( P-\alpha \, S\right) }^{2}}{P+\left( 1-\alpha \right) \, S}\quad P>\alpha \, S</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(40)
|}

donde <math display="inline">P</math> es la profundidad de precipitación bruta acumulada, y <math display="inline">{P}_{n}</math> es la profundidad de precipitación neta acumulada. La ecuación anterior es válida una vez que se han satisfecho las pérdidas iniciales ( <math display="inline">P>{I}_{a}</math>). Para valores de precipitación bruta inferiores a <math display="inline">{I}_{a}</math>, la precipitación neta es igual a cero.

La infiltración acumulada se calcula en cada paso de tiempo como la diferencia entre la precipitación neta y la bruta:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>F=P-{P}_{n}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(41)
|}

A partir de la infiltración acumulada en dos pasos de tiempo consecutivos, la tasa de infiltración se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>f=\frac{{F}_{(t+\Delta t)}-{F}_{\left( t\right) }}{\Delta t}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(42)
|}

Una diferencia con la implementación clásica del método del número de curva en la implementación de las pérdidas iniciales en Iber, es que pueden satisfacerse con la precipitación inicial o con la escorrentía superficial que haya sobre el terreno antes de que comience la precipitación. Es decir, puede haber infiltración antes de que comience a llover. Una vez que se han satisfecho las pérdidas iniciales, las pérdidas por infiltración se satisfacen únicamente con la precipitación (a diferencia de los otros métodos de infiltración, que pueden continuar infiltrando agua incluso aunque no llueva).

'''''Modelo de infiltración del SCS continuo'''''

Esta es una variante del modelo SCS que permite recuperar capacidad de infiltración durante períodos de tiempo en los que no llueve. La formulación implementada en Iber está basada en la formulación propuesta en Kannan et al. [34] para el cálculo de infiltración a nivel diario. La principal diferencia es que en Iber se implementa con un paso de tiempo igual al paso de tiempo de cálculo, mientras que en Kannan et al. se plantea únicamente para valores diarios de precipitación e infiltración.

El método se basa en recalcular, en cada paso de tiempo, la capacidad de retención del suelo ( <math display="inline">S</math>) a partir de la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{S}^{n}={S}^{n-1}+\left( {e}^{n-1}+{p}^{n-1}-{f}^{n-1}\right) \cdot \Delta t</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(43)
|}

donde <math display="inline">S</math> es la capacidad de retención del suelo en un instante determinado, <math display="inline">f</math> es la infiltración real, <math display="inline">e</math> es la evapotranspiración y <math display="inline">p</math> es la percolación. Los superíndices <math display="inline">n</math> y <math display="inline">n-</math><math>1</math> hacen referencia a los instantes de tiempo en los que se calcula cada término. Por lo tanto, la capacidad de retención del suelo disminuye a medida que infiltra el agua, y aumenta a medida que evapotranspira o que percola. Si el usuario no activa la percolación ni la evapotranspiración, el modelo no recuperará capacidad de infiltración.

A partir del valor de la capacidad de retención del suelo, se calcula la infiltración en cada paso de tiempo como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>f=\frac{h-{h}^{\ast }}{\Delta t}\, \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(44)
|}

donde <math display="inline">h</math> es el calado sobre el terreno antes de producirse la infiltración, y <math display="inline">{h}^{\ast }</math> es el calado restante una vez satisfecha la capacidad de infiltración en este paso de tiempo, el cual se evalúa como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{h}^{\ast }=\left\{ \begin{matrix}\frac{{\left( h-\alpha \, {S}_{\Delta t}\right) }^{2}}{h+\left( 1-\alpha \right) \, {S}_{\Delta t}}&&h>\alpha \, {S}_{\Delta t}\\&&\\0&&h<\alpha \, {S}_{\Delta t}\end{matrix}\right.</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(45)
|}

con:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{S}_{\Delta t}=S\, \frac{\Delta t}{86400}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(46)
|}

siendo <math display="inline">\Delta t</math> el paso de tiempo de cálculo expresado en segundos. El valor de la capacidad de retención <math display="inline">S</math> se actualiza a cada paso de tiempo y se limita entre 0 y <math display="inline">{S}_{max}</math>, siendo <math display="inline">{S}_{max}</math> la máxima capacidad de retención del terreno, calculada a partir del número de curva para condiciones antecedentes de humedad secas ( <math display="inline">C{N}_{dry}</math>), como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{S}_{max}=\frac{25400}{C{N}_{dry}}-254\\C{N}_{dry}=\frac{4,2\, CN}{10-0,058\, CN}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(47)
|}

'''''Modelo de infiltración de Horton'''''

En el modelo de Horton la tasa de infiltración potencial se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{f}_{pot}={f}_{r}+\left( {f}_{0}-{f}_{r}\right) \, {e}^{-kt}\\\, \quad \quad f=\mathrm{min}\,\left( {f}_{pot},\, \frac{h}{\Delta t}\right) \end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(48)
|}

donde <math display="inline">t</math> es el tiempo de comienzo de la precipitación, <math display="inline">{f}_{o}</math> es la tasa de infiltración inicial (máxima), <math display="inline">{f}_{r}</math> es la tasa de infiltración una vez que el suelo está completamente saturado (mínima) y <math display="inline">k</math> es un coeficiente de decaimiento que define la variación temporal de la tasa de infiltración.

La ecuación anterior no se puede aplicar directamente en un modelo como Iber en el que la precipitación puede ser intermitente, ya que la capacidad de infiltración disminuiría de forma constante con el tiempo incluso durante períodos en los cuales no hubiese precipitación ni infiltración. Por ello es necesario implementar una ecuación en la cual la capacidad de infiltración disminuya a medida que el agua se infiltra en el terreno.

En Iber se implementa de la siguiente forma. En cada paso de tiempo se calcula la infiltración potencial como la integral de <math display="inline">{f}_{pot}</math> entre <math display="inline">t</math> y <math display="inline">t+</math><math>\Delta t</math>, y se actualiza el valor de <math display="inline">{f}_{0}</math> según la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\Delta {F}_{pot}={f}_{r}\, \Delta t+\frac{\left( {f}_{0}-{f}_{r}\right) }{k}\left( 1-\right. </math><math>\left. {e}^{-k\, \Delta t}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(49)
|}

Si <math display="inline">\Delta {F}_{pot}<h</math>:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}\Delta F=\, \Delta {F}_{pot}\\{f}_{o}={f}_{r}+\left( {f}_{0}-{f}_{r}\right) \, {e}^{-k\, \Delta t}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(50)
|}

Si <math display="inline">\Delta {F}_{pot}>h</math>:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}\Delta F=h\\{f}_{o}={f}_{r}+\left( {f}_{0}-{f}_{r}\right) \, {e}^{-k\, \frac{\Delta F}{f}}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(51)
|}

Este proceso se repite en cada paso de tiempo de cálculo. En el caso de que siempre haya agua suficiente en la superficie para satisfacer la capacidad de infiltración, esta implementación se corresponde con una discretización temporal de la ecuación clásica de Horton, con un decaimiento exponencial continuo en el tiempo de la capacidad de infiltración. Sin embargo, en el caso de que no haya agua disponible para infiltrarse porque no llueve y el terreno está seco ( <math display="inline">h=</math><math>0</math>), entonces <math display="inline">\Delta F=0</math> y <math display="inline">{f}_{0}=</math><math>cte</math>. Es decir, en esta situación no se infiltra agua y la tasa de infiltración potencial deja de disminuir porque no está aumentando la humedad del suelo.

'''''Resumen de los modelos de infiltración implementados'''''

En la <span id='cite-_Ref106369539'></span>[[#_Ref106369539|Tabla 3]] se resumen las principales características de los modelos de infiltración implementados en Iber.

<span id='_Ref106369539'></span>
<span style="text-align: center; font-size: 75%;">'''Tabla 3. Formulaciones implementadas en Iber expresadas como capacidad potencial de infiltración ('''</span> <math display="inline">{\mathit{\boldsymbol{f}}}_{}^{\mathit{\boldsymbol{pot}}}</math><span style="text-align: center; font-size: 75%;">''').'''</span>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: middle;"|'''Modelo de infiltración'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: middle;"|'''Parámetros'''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: middle;"|'''Notas'''
|-
|  style="vertical-align: middle;"|Lineal
|  style="text-align: center;vertical-align: middle;"|2
|  style="vertical-align: middle;"|Incluye abstracción inicial. Capacidad de infiltración constante en el tiempo.
|-
|  style="vertical-align: middle;"|Green-Ampt
|  style="text-align: center;vertical-align: middle;"|Entre

3 y 6
|  style="vertical-align: middle;"|Incluye abstracción inicial. Si se activan la percolación y/o evapotranspiración, se considera la recuperación de capacidad de infiltración en el tiempo. En caso contrario la infiltración potencial disminuye en el tiempo.
|-
|  style="vertical-align: middle;"|SCS
|  style="text-align: center;vertical-align: middle;"|2
|  style="vertical-align: middle;"|Incluye abstracción inicial. Se incluye un parámetro para calibrar las pérdidas iniciales. La infiltración potencial disminuye en el tiempo.
|-
|  style="vertical-align: middle;"|SCS continuo
|  style="text-align: center;vertical-align: middle;"|3
|  style="vertical-align: middle;"|Incluye abstracción inicial. Se incluye un parámetro para calibrar las pérdidas iniciales. Si se activan la percolación y/o evapotranspiración, se considera la recuperación de capacidad de infiltración en el tiempo. En caso contrario la infiltración potencial disminuye en el tiempo.
|-
|  style="border-bottom: 2pt solid black;vertical-align: middle;"|Horton
|  style="border-bottom: 2pt solid black;text-align: center;vertical-align: middle;"|3
|  style="border-bottom: 2pt solid black;vertical-align: middle;"|No incluye pérdidas iniciales. La infiltración potencial disminuye en el tiempo.
|}


<span id='_Ref106600683'></span><span id='_Ref106601435'></span><span id='_Toc108335462'></span>

====4.1.5 Evapotranspiración====

La evapotranspiración potencial se calcula mediante la formulación propuesta por Blaney y Criddle, posteriormente modificada por Doorenbos y Pruitt [35]. Se ha optado por esta formulación por su sencillez y porque únicamente necesita un parámetro. La evapotranspiración potencial se calcula mediante la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>ETP=-2+b\, \left[ \xi \, \left( 0,46\, {T}_{a}+8,13\right) \right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(52)
|}

donde <math display="inline">ETP</math> es la evapotranspiración potencial en mm/día, <math display="inline">{T}_{a}</math> es la temperatura del aire en grados centígrados, <math display="inline">\xi</math>  es el porcentaje de horas diurnas totales en el día considerado sobre el total de horas diurnas en un año (365 x 12) y <math display="inline">b</math> es un parámetro del modelo que debe calibrarse. Como término de referencia para el valor de <math display="inline">\xi</math> , en el equinocio de primavera u otoño (12 horas de día), se obtiene un valor de <math display="inline">\xi =</math><math>0,27</math>. En la versión actual de Iber valor de <math display="inline">\xi</math>  está prefijado a 0,27 (12 horas de día de sol), y no puede ser modificado por el usuario.

En cuanto a los valores habituales de <math display="inline">b</math>, suelen estar comprendidos entre 0,7 y 2. El usuario puede modificar este valor a su conveniencia, con el fin de calibrar el modelo.

'''''Evapotranspiración real con el modelo de infiltración de Green-Ampt'''''

En el caso de que se seleccione el modelo de Green-Ampt la evapotranspiración real, en mm/día, se calcula a partir de la potencial como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>e=ETP\frac{\theta }{\phi }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(53)
|}

siendo <math display="inline">\phi</math>  la porosidad del suelo y <math display="inline">\theta</math>  la humedad del suelo, que debe estar comprendida entre 0 y <math display="inline">\phi</math> . Es decir, si el suelo está completamente saturado, la evapotranspiración real es igual a la potencial, y si el suelo está totalmente seco, la evapotranspiración real es igual a cero. Entre estos dos estados extremos, se establece una relación lineal.

Además, se incluye una formulación para considerar el estrés hídrico en el cálculo de la evapotranspiración real. Cuando la humedad del subsuelo baja por debajo de un determinado valor umbral, las plantas tienen una gran dificultad para extraer el agua disponible y se dice que están bajo estrés hídrico. En estas condiciones se produce una disminución de la tasa de evapotranspiración. Si se activa la opción de ‘Estrés Hídrico’, se utiliza el método conocido como Coeficiente Único del Cultivo (''single crop coefficient'') para calcular la evapotranspiración como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>e={K}_{s}\, {K}_{c}\, ETP</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(54)
|}

dónde <math display="inline">{K}_{s}</math> es el coeficiente de estrés hídrico, un factor adimensional comprendido entre 0 y 1 que es función del contenido en agua del suelo, y <math display="inline">{K}_{c}</math> es un coeficiente que depende del tipo de cultivo y del estado del cultivo (época del año), comprendido entre 0,3 y 1,2.

El valor de <math display="inline">{K}_{c}</math> lo debe introducir el usuario en función del tipo de cultivo, y puede ser variable en el tiempo. El coeficiente <math display="inline">{K}_{s}</math> se evalúa en función del punto de marchitez (''wilting point'') y de la capacidad de campo (''field capacity''), como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{K}_{s}=\left\{ \begin{matrix}0&si&\theta \leq {\theta }_{wp}\\\frac{\theta -{\theta }_{wp}}{\left( {\theta }_{fc}-{\theta }_{wp}\right) \left( 1-p\right) }&si&\theta <{\theta }_{t}\\1&si&\theta \geq {\theta }_{t}\end{matrix}\right. \, \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(55)
|}

siendo <math display="inline">{\theta }_{t}={\theta }_{fc}\left( 1-p\right) +{\theta }_{wp}\, p</math>, la humedad del suelo por debajo de la cual comienza a haber estrés hídrico, <math display="inline">{\theta }_{wp}</math> el punto de marchitez, <math display="inline">{\theta }_{fc}</math> la capacidad de campo, y <math display="inline">p</math> el factor de agotamiento (''depletion factor''), que debe estar comprendido entre 0 y 1.

'''''Evapotranspiración real con el modelo de infiltración de SCS continuo'''''

En el caso de que se active el modelo de infiltración del SCS continuo, la evapotranspiración real se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>e=ETP\, \cdot \, exp\left( -B\frac{S}{{S}_{max}}\right) \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(56)
|}

donde <math display="inline">B</math> es un coeficiente de agotamiento a definir por el usuario, <math display="inline">S</math> es la capacidad de retención del terreno en el instante considerado, y <math display="inline">{S}_{max}</math> es la máxima capacidad de retención del terreno, calculada a partir del número de curva ( <math display="inline">CN</math>) para condiciones antecedentes de humedad secas ( <math display="inline">C{N}_{dry}</math>), como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{S}_{max}=\frac{25400}{C{N}_{dry}}-254\\C{N}_{dry}=\frac{4,2\, CN}{10-0,058\, CN}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(57)
|}

donde <math display="inline">CN</math> es el número de curva introducido por el usuario como parámetro del modelo SCS continuo. El coeficiente de agotamiento <math display="inline">B</math> suele estar comprendido en el rango <math display="inline">\left[ 0,2\right]</math> .

'''''Evapotranspiración real con los modelos de infiltración Lineal, SCS y Horton'''''

Cuando se activa un modelo de infiltración diferente a Green-Ampt o a SCS continuo, únicamente se puede calcular la evapotranspiración si al mismo tiempo se activa el flujo subsuperficial y se selecciona el modelo agregado. En caso de que no se active el flujo subsuperficial no tiene sentido calcular la evapotranspiración con los modelos de infiltración Lineal, SCS u Horton, ya que ninguno de ellos hace un seguimiento de la humedad del suelo ni puede recuperar capacidad de infiltración.

En caso de que se active el modelo agregado de flujo subsuperficial en combinación con los modelos de infiltración Lineal, SCS o Horton, la evapotranspiración real se calcula a partir de la potencial como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>e=ETP\frac{{V}_{s}}{{V}_{s,max}\, }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(58)
|}

siendo <math display="inline">{V}_{s}</math> el volumen de agua en el subsuelo, y <math display="inline">{V}_{s,max}</math> el volumen de agua máximo en el subsuelo. Para realizar el cálculo de la evapotranspiración, el cálculo del volumen de agua en el subsuelo se realiza a nivel de depósito subterráneo, tal como se detalla en la descripción del modelo agregado de flujo subsuperficial.

<span id='_Ref106600689'></span><span id='_Ref106601442'></span><span id='_Toc108335463'></span>

====4.1.6 Percolación====

Al igual que la evapotranspiración, la percolación se calcula de forma diferente en función del modelo de infiltración seleccionado.

'''''Percolación con el modelo de infiltración de Green-Ampt'''''

Si se selecciona el modelo de infiltración de Green-Ampt, la percolación se calcula mediante la ecuación propuesta por Famiglietti y Wood [36]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>p={k}_{p}{\left( \frac{\theta }{\phi }\right) }^{\left( 3+\frac{2}{\lambda }\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(59)
|}

donde <math display="inline">p</math> es la tasa de percolación, <math display="inline">{k}_{p}</math> es la permeabilidad de las capas inferiores del suelo en las cuales se produce la percolación, <math display="inline">\phi</math>  la porosidad del suleo, <math display="inline">\theta</math>  la humedad del suelo (que debe estar comprendida entre 0 y <math display="inline">\phi</math> ) y <math display="inline">\lambda</math>  es un parámetro del modelo (''soil size distribution index''). Los valores habituales de <math display="inline">\lambda</math>  están comprendidos entre 0,05 y 0,5.

'''''Percolación con el modelo de infiltración de SCS continuo'''''

Si se selecciona el modelo de infiltración SCS continuo, la percolación se calcula mediante la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>p={k}_{p}{\left( 1-\frac{S}{{S}_{max}}\right) }^{\left( 3+\frac{2}{\lambda }\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(60)
|}

donde <math display="inline">S</math> es la capacidad de retención del terreno en el instante considerado, y <math display="inline">{S}_{max}</math> es la máxima capacidad de retención del terreno, calculada a partir del número de curva para condiciones antecedentes de humedad secas.

'''''Percolación con los modelos de infiltración Lineal, SCS y Horton'''''

Cuando se activa un modelo de infiltración diferente a Green-Ampt o a SCS continuo, únicamente se puede calcular la percolación si al mismo tiempo se activa el flujo subsuperficial y se selecciona el modelo agregado. En caso de que no se active el flujo subsuperficial no tiene sentido calcular la percolación con los modelos de infiltración Lineal, SCS o Horton, ya que ninguno de ellos hace un seguimiento de la humedad del suelo ni puede recuperar capacidad de infiltración.

En el caso de que se active el modelo agregado de flujo subsuperficial en combinación con los modelos de infiltración Lineal, SCS o Horton, la percolación se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>p={k}_{p}{\left( \frac{{V}_{s}}{{V}_{s,max}}\right) }^{\left( 3+\frac{2}{\lambda }\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(61)
|}

siendo <math display="inline">{V}_{s}</math> el volumen de agua en el subsuelo, y <math display="inline">{V}_{s,max}</math> el volumen de agua máximo en el subsuelo. Para realizar el cálculo de la percolación, el cálculo del volumen de agua en el subsuelo se realiza a nivel de depósito subterráneo, tal como se detalla en la descripción del modelo agregado de flujo subsuperficial.

<span id='_Toc108335464'></span>

====4.1.7 Flujo superficial====

En el módulo hidrológico de Iber se incluyen dos modelos de flujo subsuperficial: un modelo agregado y un modelo distribuido. El modelo distribuido debe utilizarse conjuntamente con el modelo de infiltración de Green-Ampt, ya que hereda de dicho modelo varios de sus parámetros. El modelo agregado puede utilizarse con cualquiera de las formulaciones de infiltración definidas en la sección <span id='cite-_Ref106600646'></span>[[#_Ref106600646|4.1.4]].

'''''Modelo agregado'''''

Este modelo es más sencillo que el modelo distribuido. En muchos casos, el modelo agregado puede ser más atractivo de cara al usuario que quiere considerar el flujo subterráneo en el cálculo de forma sencilla, sin necesidad de resolver un modelo complejo que depende de nuevos parámetros variables espacialmente y, en general, difíciles de caracterizar.

El modelo agregado realiza un balance global de agua en el subsuelo de la cuenca considerando los aportes por infiltración desde el terreno y las pérdidas por percolación, evapotranspiración y exfiltración. La componente de exfiltración representa la salida del agua subterránea hacia la red de cauces fluviales, mientras que las componentes de percolación y evapotranspiración representan agua que se pierde del sistema. De esta forma, el agua infiltrada en la cuenca pasa a formar parte de un depósito subterráneo desde el cual puede percolar a acuíferos inferiores (en cuyo caso sale del sistema y se pierde su seguimiento), perderse por evapotranspiración, o salir a los cauces fluviales de la cuenca (exfiltración), en cuyo caso pasa a aumentar el caudal de los cauces superficiales.

Se puede considerar toda la cuenca como un único depósito, o puede dividirse la superficie de la cuenca en varios depósitos con propiedades diferenciadas. Independientemente del número de depósitos en los que se divida la cuenca, el balance de masa de agua en cada depósito se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{d{V}_{s}}{dt}={Q}_{inf}-{Q}_{per}-{Q}_{evap}-{Q}_{out}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(62)
|}

donde <math display="inline">{V}_{s}</math> es el volumen de agua en el depósito (subsuelo), <math display="inline">{Q}_{out}</math> es el caudal que sale desde el subsuelo hacia los cauces fluviales (exfiltración), <math display="inline">{Q}_{inf}</math> es el caudal infiltrado desde el terreno al subsuelo, <math display="inline">{Q}_{per}</math> es el caudal que percola desde el subsuelo hacia capas más profundas del terreno, <math display="inline">{Q}_{evap}</math> es el caudal de agua que se pierde por evapotranspiración.

El caudal de infiltración ( <math display="inline">{Q}_{inf}</math>) se calcula a partir de cualquiera de los modelos de infiltración presentados en la sección <span id='cite-_Ref106600646'></span>[[#_Ref106600646|4.1.4]]. Los caudales de percolación ( <math display="inline">{Q}_{per}</math>) y de evapotranspiración ( <math display="inline">{Q}_{evap}</math>) se calculan a partir de sus respectivos valores por unidad de superficie (formulaciones descritas en la secciones <span id='cite-_Ref106600683'></span>[[#_Ref106600683|4.1.5]] y <span id='cite-_Ref106600689'></span>[[#_Ref106600689|4.1.6]]), simplemente multiplicándolos por la superficie de la cuenca asociada al depósito correspondiente ( <math display="inline">{S}_{tot}</math>).

El caudal exfiltrado ( <math display="inline">{Q}_{out}</math>) se calcula según la siguiente expresión:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{Q}_{out}=A\, \left( {V}_{s}-B\, {S}_{tot}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(63)
|}

donde <math display="inline">{S}_{tot}</math> es la superficie de la cuenca asociada al depósito correspondiente y <math display="inline">A</math> y <math display="inline">B</math> son parámetros asociados a dicho depósito. El caudal exfiltrado se aporta a la escorrentía superficial, en tramos de ríos definidos previamente por el usuario.

Para aplicar el modelo agregado el usuario debe definir la extensión espacial de cada depósito (a partir de la cual el modelo calculará <math display="inline">{S}_{tot}</math>), los parámetros <math display="inline">A</math> y <math display="inline">B</math> para cada depósito, y la zona de exfiltración de cada depósito (en la cual se repartirá el caudal exfiltrado <math display="inline">{Q}_{out}</math>).

Además, el usuario debe definir, para cada depósito, el volumen máximo de agua en el subsuelo a efectos de flujo subsuperficial y el contenido inicial de agua. Estos valores se definen como una profundidad máxima de subsuelo ( <math display="inline">{d}_{max}</math>) y una saturación inicial ( <math display="inline">{S}_{r,0}</math>). A partir de estos valores se calcula el volumen inicial y máximo de agua en cada depósito como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{V}_{s,max}={S}_{tot}\, {d}_{max}\\{V}_{s,0}={S}_{r,0}{V}_{s,max}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(64)
|}

donde <math display="inline">{V}_{s,max}</math> es el volumen máximo de agua en el depósito y <math display="inline">{V}_{s,0}</math> el volumen inicial de agua en el depósito.

'''''Modelo distribuido'''''

El modelo de flujo subterráneo distribuido se basa en la ecuación de Boussinesq para flujo 2D en medio poroso saturado no confinado. Para evaluar el flujo de agua subsuperficial y la distribución espacio-temporal del agua en el suelo se resuelve la siguiente ecuación diferencial de conservación de la masa:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\phi \frac{\partial b}{\partial t}=\frac{\partial }{\partial x}\left( T\frac{\partial H}{\partial x}\right) +</math><math>\frac{\partial }{\partial y}\left( T\frac{\partial H}{\partial y}\right) +f-p-e-</math><math>g</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(65)
|}

siendo <math display="inline">\phi</math>  la porosidad del suelo, <math display="inline">b</math> el espesor saturado de suelo, <math display="inline">T</math> la transmisividad horizontal, <math display="inline">H</math> la carga hidráulica, <math display="inline">f</math> la infiltración de agua desde el terreno al subsuelo, <math display="inline">p</math> la percolación, <math display="inline">e</math> la evapotranspiración y <math display="inline">g</math> la exfiltración de agua del subsuelo a la superficie.

Dado que en Iber no se resuelve el flujo en dirección vertical, el espesor de la capa de suelo saturada ( <math display="inline">b</math>) se relaciona con el contenido de humedad promediado verticalmente en la capa de suelo ( <math display="inline">\theta</math> ), mediante la siguiente relación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>b=\theta \frac{{d}_{s}}{\phi }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(66)
|}

donde <math display="inline">{d}_{s}</math> es el espesor total de la capa de suelo. Es decir, el valor de <math display="inline">b</math> se obtiene agrupando todo el contenido de agua de la capa de suelo en una capa horizontal situada justo bajo el terreno (<span id='cite-_Ref106367812'></span>[[#_Ref106367812|Fig. 23]]). Esta simplificación es únicamente válida para espesores de suelo pequeños.

La transmisividad horizontal se considera variable en función del contenido de humedad del suelo, calculándose mediante la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>T={k}_{h}\, {d}_{s}\, {e}^{-\left( \frac{\phi -\theta }{m}\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(67)
|}

donde <math display="inline">{k}_{h}</math> es la permeabilidad saturada horizontal del suelo y <math display="inline">m</math> es un parámetro que en la versión actual se toma constante e igual a 0,1. Los valores de porosidad del suelo ( <math display="inline">\phi</math> ) y de espesor total de la capa de suelo ( <math display="inline">{d}_{s}</math>) se heredan del modelo de Green-Ampt, mientras que la humedad del suelo ( <math display="inline">\theta</math> ) es una variable que se va actualizando a medida que se resuelve la ecuación de conservación de la masa de agua en el subsuelo, y que también se comparte con el modelo de Green-Ampt.

En un suelo completamente saturado ( <math display="inline">\theta =\phi</math> ), la transmisividad es máxima e igual a la permeabilidad saturada por el espesor de suelo <math display="inline">T=</math><math>{k}_{h}\, {d}_{s}</math>.

La carga hidráulica se calcula como la cota inferior de la capa de suelo más el espesor saturado:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>H={z}_{b}-{d}_{s}+b</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(68)
|}

donde <math display="inline">{z}_{b}</math> es la cota del terreno.

Cuando se utiliza el modelo de flujo subterráneo distribuido, la tasa de infiltración ( <math display="inline">f</math>) debe calcularse siempre con la formulación de Green-Ampt descrita anteriormente, de la cual hereda los siguientes parámetros: porosidad del suelo ( <math display="inline">\phi</math> ), espesor total de la capa de suelo ( <math display="inline">{d}_{s}</math>) y humedad inicial del suelo ( <math display="inline">{\theta }_{0}</math>). La permeabilidad saturada del suelo utilizada para el cálculo de la transmisividad horizontal ( <math display="inline">{k}_{h}</math>) no se hereda del modelo de Green-Ampt ya que la infiltración es un fenómeno fundamentalmente en la dirección vertical, mientras que el flujo subsuperficial representado por la ecuación de Boussinesq se produce fundamentalmente en dirección horizontal (paralela al terreno). Debido a la anisotropía de muchos suelos, se permite utilizar un valor diferente de permeabilidad saturada en el modelo de Green-Ampt y en la ecuación de Boussinesq. Evidentemente, el usuario puede utilizar el mismo valor si lo considera oportuno.

Las tasas de evapotranspiración ( <math display="inline">e</math>) y percolación ( <math display="inline">p</math>) se calculan según las formulaciones descritas en las secciones <span id='cite-_Ref106601435'></span>[[#_Ref106601435|4.1.5]] y <span id='cite-_Ref106601442'></span>[[#_Ref106601442|4.1.6]] para el modelo de Green-Ampt. La exfiltración ( <math display="inline">e</math>) se calcula imponiendo que el espesor de la capa de suelo saturada debe ser inferior o igual al espesor total de la capa de suelo ( <math display="inline">b\leq {d}_{s}</math>), lo cual equivale a decir que la humedad del suelo no puede ser superior a su porosidad ( <math display="inline">\theta \leq \phi</math> ). De esta forma, si en algún momento del cálculo se obtiene un valor de <math display="inline">b</math> superior a <math display="inline">{d}_{s}</math>, la diferencia entre ambos se impone como un caudal de exfiltración que sale del subsuelo hacia los cauces fluviales, al mismo tiempo que el valor de <math display="inline">b</math> se hace igual a <math display="inline">{d}_{s}</math>.

La ecuación de conservación de la masa de agua en el subsuelo se resuelve mediante el método de volúmenes finitos en la misma malla no estructurada en la que se resuelven las ecuaciones de aguas someras.

<span id='_Ref106601634'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image24.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 24. Esquema de los procesos considerados en el modelo distribuido de flujo subsuperficial.'''</span>


En cuanto a la implementación y acoplamiento del flujo subsuperficial y del flujo de escorrentía superficial, el proceso de cálculo (<span id='cite-_Ref106601634'></span>[[#_Ref106601634|Fig. 24]]) en cada paso de tiempo es el siguiente:

:1. Resolver las ecuaciones de aguas someras 2D (escorrentía superficial) sin considerar aportes de lluvia ni de infiltración. Actualización de los calados y velocidades del agua en superficie.

:2. Resolver la ecuación de flujo subsuperficial 2D. Actualización de la humedad media en la capa de suelo ( <math display="inline">\theta</math> ) y del espesor saturado equivalente ( <math display="inline">b</math>).

:3.  Si el espesor saturado es superior al espesor de la capa de suelo ( <math display="inline">b\geq \, {d}_{s}</math>) se está produciendo exfiltración (agua subsuperficial saliendo a la superficie), en cuyo caso se actualizan las variables de flujo de la siguiente manera:
:* <math display="inline">b\, =\, {d}_{s}</math>
:* <math display="inline">h\, =\, h\, +\, \left( b-{d}_{s}\right) \, \phi</math>
:* El caudal de exfiltración en este caso sería igual a <math display="inline">g=</math><math>\frac{\left( b-{d}_{s}\right) \, \phi }{\Delta t}</math>

:4. Añadir el aporte por lluvia
:* <math display="inline">h\, =\, h\, +\, R\, \Delta t</math>

:5. Añadir las pérdidas por infiltración según el modelo de Green-Ampt
:* <math display="inline">{f}_{pot}={k}_{s}\left( 1+\frac{\left( h+\psi \right) }{b}\right)</math>
:* <math display="inline">{F}_{max}=\mathrm{min}\,\left[ h,\, \left( {d}_{s}-b\right) \phi \right]</math>
:* <math display="inline">F=\mathrm{min}\,\left[ {f}_{pot}\, \Delta t,\, {F}_{max}\right] \,</math>
:* <math display="inline">h=h-F</math>
:* <math display="inline">b=b+F/\phi</math>

:6. Añadir las pérdidas por evapotranspiración y por percolación

<span id='_Toc108335465'></span>

====4.1.8 Resultados de las variables hidrológicas====

Además de las variables hidrodinámicas de calado y velocidad del agua superficial (escorrentía superficial), al activar los procesos hidrológicos se escriben los siguientes resultados:

<span id='OLE_LINK12'></span><span id='OLE_LINK13'></span>
:* '''Hidrogramas''': el usuario puede definir diferentes contornos de salida de agua, así como secciones de aforo internas. Los hidrogramas calculados a través de dichos contornos de salida y de las secciones internas de aforo se escribirán como series temporales en el fichero ''Q_hydrograph.rep''.
:* '''Series temporales de variables hidrológicas''': las series temporales de las siguientes variables se escriben en el fichero ''Hydrology_timeseries.rep'':
::* Profundidad de escorrentía superficial (mm) instantánea y promediada en toda la cuenca
::* Profundidad de precipitación (mm) acumulada y promediada en toda la cuenca
::* Intensidad de precipitación (mm/h) instantánea y promediada en toda la cuenca
::* Profundidad de infiltración (mm) acumulada y promediada en toda la cuenca
::* Tasa de infiltración (mm/h) instantánea y promediada en toda la cuenca
::* Volumen de agua en el subsuelo (m<sup>3</sup>) instantáneo
::* Espesor saturado del subsuelo (mm) instantáneo y promediado en toda la cuenca
::* Saturación del subsuelo (-) instantánea y promediada en toda la cuenca
::* Profundidad de evapotranspiración (mm) acumulada y promediada en toda la cuenca
::* Tasa de evapotranspiración (mm/h) instantánea y promediada en toda la cuenca
::* Profundidad de percolación (mm) acumulada y promediada en toda la cuenca
::* Tasa de percolación (mm/h) instantánea y promediada en toda la cuenca
::* Caudal exfiltrado (m<sup>3</sup>/s) instantáneo
:* '''Campos espaciales de variables hidrológicas''': en la interfaz de post-proceso, se muestra un nuevo grupo de variables (''Hydrology'') en el que se incluyen los campos espaciales de las siguientes variables:
::* Intensidad de precipitación (mm/h)
::* Profundidad de precipitación (mm) acumulada desde el inicio del cálculo
::* Tasa de infiltración (mm/h)
::* Profundidad de infiltración (mm) acumulada desde el inicio del cálculo
::* Espesor saturado del suelo (m)
::* Altura piezométrica en el subsuelo (m)
::* Caudal unitário subsuperficial (m<sup>2</sup>/s)

<span id='_Toc108335466'></span>

===4.2 Interfaz de usuario===

En esta sección se describen los diferentes menús y ventanas desde las que se pueden activar los procesos hidrológicos descritos en las secciones anteriores y definir los correspondientes parámetros.

<span id='_Toc108335467'></span>

====4.2.1 Precipitación====

La precipitación puede definirse de forma variable en espacio y tiempo, ya sea mediante asignación manual de hietogramas a diferentes zonas del modelo o de forma automática mediante la definición de estaciones pluviométricas o de rasters de precipitación.

La activación y el modo de la lectura pluviométrica se realiza a través del menú '''Datos >> Datos del Problema >> Hidrología >> Precipitación ''' (<span id='cite-_Ref107216544'></span>[[#_Ref107216544|Fig. 25]]). Dentro de dicho menú, se dispone de 4 opciones:

:* Sin precipitación
:* Hietograma
:* Rasters
:* Interpolación raster

Aunque se definan y se asignen los hietogramas o los rasters de lluvia, éstos nos serán utilizados en el cálculo si se mantiene la opción “Sin precipitación”.

<span id='_Ref107216544'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image25.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 25. Ventana de selección del modo de lectura de datos pluviométricos.'''</span>


'''''Lluvia definida por hietogramas'''''

La definición de hietogramas se realiza a través del menú '''Datos >> Procesos hidrológicos >> Lluvia >> Definición de hietogramas''' (<span id='cite-_Ref107216543'></span>[[#_Ref107216543|Fig. 26]]a). En la ventana correspondiente, el usuario puede introducir las series temporales de datos pluviométricos a través del instante temporal (s) y la precipitación (mm/h). Además, es necesario establecer la ubicación de la estación pluviométrica en la que se registran dichos datos (<span id='cite-_Ref107216543'></span>[[#_Ref107216543|Fig. 26]]a). Si se selecciona la ubicación de la estación a través de la interfaz de GiD, no es posible establecer un punto fuera de la superficie del modelo, pero escribiendo las coordenadas que se desean, sí se puede fijar la ubicación fuera de la geometría del modelo.

<span id='_Ref107216543'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image26.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 26. Definición y asignación de hietogramas: (a) ventana de definición de hietogramas; (b) ventana de selección del modo de lectura de datos pluviométricos.'''</span>


Una vez definidos los hietogramas, el usuario dispone de dos opciones:

:* Se puede proceder a la asignación manual de cada hietograma en diferentes regiones del modelo, a través del menú '''Datos >> Procesos hidrológicos >> Lluvia >> Asignación de hietogramas''' (<span id='cite-_Ref107216543'></span>[[#_Ref107216543|Fig. 26]]b). De esta manera, se asigna la precipitación definida en el hietograma de forma homogénea en las regiones del modelo seleccionadas. Si ya se ha generado la malla, también es posible asignar los hietogramas directamente a los elementos de la malla.
:* En las regiones en las cuales no se haya asignado manualmente un hietograma, el programa atribuirá automáticamente la pluviometría a cada zona del modelo a través de una interpolación del tipo ''Nearest Neighbour'', considerando para ello la ubicación de cada pluviómetro.

'''''Lluvia definida por rasters'''''

Si seleccionamos esta opción, la lluvia se definirá a partir de archivos raster en formato ASCII Grid de ArcInfo (<span id='cite-_Ref107216545'></span>[[#_Ref107216545|Fig. 27]]). En cada archivo se definirá un campo de precipitación distribuido espacialmente. Todos los rasters deben figurar en una carpeta llamada ''Rain_rasters'', ubicada dentro de la carpeta del proyecto, para garantizar su correcta lectura. Se aclara, por tanto, que esta carpeta debe ser creada por el usuario y que no se genera por defecto con el modelo. Además, todos los ficheros rasters deben tener el mismo número de columnas y filas, el mismo tamaño de pixel y el mismo origen.

Dentro de la definición de la lluvia por rasters tenemos dos alternativas:

:* '''Rasters'''. Se toma el mismo valor de lluvia entre cada dos rasters, constante en el tiempo. Esta es la forma análoga a la definición temporal de la lluvia en hietogramas, i.e. la variación temporal de la precipitación consiste en bloques de intensidad constante. Si se selecciona esta opción es posible elegir el tipo de datos de lluvia incluido en los rasters, que pueden ser valores de intensidad de precipitación (mm/h) o de profundidad de precipitación (mm) entre cada dos rasters.
:* '''Rasters interpolados'''. Se realiza una interpolación lineal en el tiempo de la precipitación a partir de los datos de dos archivos raster consecutivos. En este caso, los raster presentan siempre valores de intensidad de lluvia (mm/h).

En cualquiera de los casos anteriores, el usuario debe definir una tabla en la que se asigna cada archivo raster a un instante de tiempo. Los archivos raster pueden tener cualquier nombre, y en la tabla se debe incluir el nombre completo del archivo (incluyendo la extensión).

<span id='_Ref107216545'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image27.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 27. Ventana de definición de series temporales de archivos raster de precipitación.'''</span>


<span id='_Toc108335468'></span>

====4.2.2 Modelo y parámetros de infiltración====

Iber permite la asignación de parámetros de infiltración de manera manual o automática a partir de un archivo raster. En ambos casos se permite elegir entre diferentes modelos de infiltración (Lineal, SCS, Horton, Green-Ampt, SCS continuo) y los parámetros de cada modelo pueden ser variables espacialmente.

Como siempre en Iber, la asignación automática de propiedades se realiza directamente sobre la malla de cálculo, que por lo tanto deberá estar creada con anterioridad. Asimismo, cada vez que se cree la malla de cálculo, deben realizarse todas las asignaciones automáticas de nuevo.

La elección del modelo de pérdidas por infiltración se realiza desde el menú '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración''' (<span id='cite-_Ref106602994'></span>[[#_Ref106602994|Fig. 28]])'''. '''Los parámetros relativos a cada modelo de infiltración se pueden asignar de forma manual o automática, según se detalla a continuación.

<span id='_Ref106602994'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image28.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 28. Ventana de selección del modelo de infiltración.'''</span>


'''''Asignación manual de parámetros de infiltración'''''

Para realizar una asignación manual de los parámetros de infiltración se debe seleccionar la opción Activar en '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración. '''Una vez hecho esto aparecerán los modelos de infiltración disponibles.

Una vez seleccionado un modelo, la asignación manual de parámetros se realiza desde el menú '''Datos >> Procesos Hidrológicos >> Pérdidas por infiltración >> Asignación manual'''. En esta ventana el usuario puede introducir los parámetros correspondientes a cada modelo de infiltración (<span id='cite-_Ref106603675'></span>[[#_Ref106603675|Fig. 29]]a) y asignarlos manualmente a diferentes superficies de la geometría o elementos de la malla. Desde la ventana que se muestra en la <span id='cite-_Ref106603675'></span>[[#_Ref106603675|Fig. 29]]a pueden asignarse los parámetros de cualquier modelo de infiltración, pero Iber calculará las pérdidas únicamente con el modelo seleccionado en '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración'''.

<span id='_Ref106603675'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image29.png|550px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 29. (a) Ventana de selección manual del modelo de infiltración. (b) Ventana de definición de automática de zonas de infiltración.'''</span>


'''''Asignación automática a partir de archivos raster con los valores de cada parámetro'''''

Iber permite la asignación automática de parámetros de infiltración a partir de archivos raster georreferenciados en formato ASCII Grid de ArcInfo. Se debe preparar un archivo raster por cada parámetro, que defina la distribución espacial del valor de dicho parámetro en la zona de estudio.

Para realizar la asignación automática a partir de dichos archivos raster se debe seleccionar la opción'' Activar ''en '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración'''. Una vez hecho esto, aparecerán los modelos de infiltración disponibles. Una vez seleccionado un modelo de infiltración, la asignación automática de parámetros se realiza desde el menú '''Datos >> Procesos Hidrológicos >> Pérdidas por infiltración >> Asignación automática por parámetros'''.

Se presenta, en primer lugar, la opción de escoger el modelo de infiltración para el cual queremos asignar parámetros. La asignación automática puede realizarse para los modelos Lineal, SCS y Green-Ampt. Una vez que se elige el modelo de infiltración, se abre un directorio para escoger el archivo raster correspondiente al parámetro que se introducirá. A continuación, se describen los casos posibles:

:a) '''Si se elige el modelo SCS''', únicamente se pedirá cargar un archivo raster, correspondiente al número de curva. En dicho archivo se especificarán directamente en cada celda los valores de número de curva correspondientes.
:b) '''Si se elige el modelo Lineal''', se pedirá cargar dos archivos raster en el siguiente orden:
::*   Archivo correspondiente a las pérdidas iniciales (mm)
::*   Archivo de tasas de infiltración potencial constantes en tiempo (mm/h)
:c) '''Si se elige el modelo Green-Ampt''', se pedirá cargar seis archivos raster, en el siguiente orden:
::*   Succión en región no saturada (mm)
::*   Porosidad total
::*   Saturación inicial
::*   Permeabilidad saturada del suelo en la dirección vertical (mm/h)
::*   Pérdidas iniciales (mm)
::*   Espesor total de la capa de suelo (m)

Al igual que en la asignación manual de parámetros, pueden asignarse de forma automática los parámetros de cualquier modelo de infiltración, pero Iber calculará las pérdidas únicamente con el modelo seleccionado en '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración.'''

'''''Asignación automática de parámetros de infiltración por zonas'''''

Adicionalmente a la asignación automática de parámetros de infiltración descrita en la sección anterior, Iber también permite asignar los parámetros de infiltración de forma distribuida espacialmente por zonas homogéneas. De esta forma, en lugar de definir cada parámetro por separado, se definen manualmente zonas con un conjunto de parámetros homogéneos (en cada zona) y posteriormente se asignan estas zonas al modelo mediante un archivo raster en el cual se incluye la distribución espacial de dichas zonas (y no de cada parámetro independiente). Esta forma de asignar parámetros se llama en Iber ''Asignación automática por zonas''.

Para asignar los parámetros de esta manera, en primer lugar, es necesario seleccionar la opción'' Encendido [por zonas] ''en '''Datos >> Datos del Problema >> Hidrología >> Pérdidas por infiltración'''. Una vez hecho esto aparecerán los modelos de infiltración disponibles, debiendo seleccionar uno de ellos.

Una vez seleccionado el modelo de infiltración, se definen las zonas de infiltración desde el menú '''Datos >> Procesos Hidrológicos >> Pérdidas por infiltración >> Asignación automática por zonas >> Definición de zonas de infiltración'''. Desde aquí se accede a una ventana en la que se pueden crear tantas zonas como se quiera. Para cada zona hay que seleccionar el modelo de infiltración y los parámetros correspondientes (<span id='cite-_Ref106603675'></span>[[#_Ref106603675|Fig. 29]]b). Una vez definidas las zonas de infiltración se deben asignar al modelo mediante el menú '''Datos >> Procesos Hidrológicos >> Pérdidas por infiltración >> Asignación automática por zonas >> Asignación de zonas de infiltración'''. Al seleccionar esta opción, se abrirá una ventana desde la cual será necesario seleccionar un archivo raster en el que se incluya la distribución espacial de las zonas de infiltración. De esta forma, cada zona del modelo tendrá asignada una zona de infiltración con sus parámetros correspondientes.

<span id='_Toc108335469'></span>

====4.2.3 Evapotranspiración====

Solo puede activarse la evapotranspiración si se activa alguna de las siguientes opciones en '''Datos >> Datos del Problema >> Hidrología''':

:* Modelo de infiltración de Green-Ampt
:* Modelo de infiltración de SCS continuo
:* Modelo de flujo subterráneo agregado

Desde esta ventana también se puede determinar el coeficiente de corrección para el cálculo de la evapotranspiración potencial (comprendido entre 0,7 y 2) y, en el caso de haber elegido el modelo de infiltración de Green-Ampt, activarse también la formulación de estrés hídrico (''Soil Water Stress'').

Para el cálculo de la evapotranspiración potencial es necesario definir la temperatura del aire. La definición de series temporales de temperatura y su asignación a diferentes zonas del modelo puede realizarse desde los menús '''Datos >> Atmosféricas >> Variables Atmosféricas''' y '''Datos >> Atmosféricas >> Asignación Var Atmosféricas.'''

'''''Estrés hídrico'''''

La formulación de estrés hídrico (''Soil Water Stress'') solo está disponible en el caso de que se seleccione Green-Ampt como modelo de infiltración y se active la evapotranspiración. Éstas se seleccionan desde '''Datos >> Datos del Problema >> Hidrología.'''

Para utilizar esta formulación, una vez activada en '''Datos >> Datos del Problema >> Hidrología,''' es necesario definir uno o varios tipos de cultivo desde el menú '''Datos >> Procesos Hidrológicos >> Estrés hídrico del suelo >> Definición de cultivo'''. Desde aquí se accede la ventana de definición de cultivos, en la que se pueden crear tantos cultivos como se quiera. Para cada cultivo creado hay que definir los siguientes parámetros: capacidad de campo, punto de marchitez, factor de agotamiento y <math display="inline">{K}_{c}</math>. Dado que <math display="inline">{K}_{c}</math> es un coeficiente que depende del tipo de cultivo y de su estado, y por lo tanto puede ser variable en función de la época del año, el programa permite definir una serie temporal de dicho parámetro, con el fin de poder incluir su variación estacional.

Una vez definidos los cultivos que sean necesarios, se deben asignar al modelo mediante el menú '''Datos >> Procesos Hidrológicos >> Estrés hídrico del suelo >> Asignación de cultivo. '''Esta opción da acceso a la ventana de asignación de cultivos, desde la cual podemos seleccionar cualquiera de los cultivos creados y asignarlo manualmente a diferentes zonas del modelo.

<span id='_Toc108335470'></span>

====4.2.4 Percolación====

Solo puede activarse la percolación si se activa alguna de las siguientes opciones en '''Datos >> Datos del Problema >> Hidrología''':

:* Modelo de infiltración de Green-Ampt
:* Modelo de infiltración de SCS continuo
:* Modelo de flujo subterráneo agregado

El cálculo de la percolación se activa o desactiva a través del menú '''Datos >> Datos del Problema >> Hidrología''', donde también puede determinarse el índice de distribución de poros (comprendido entre 0,05 y 0,5) y la permeabilidad para el cálculo de la percolación.

<span id='_Toc108335471'></span>

====4.2.5 Flujo subsuperficial====

Iber dispone de dos modelos para evaluar el flujo de agua a través del subsuelo: un modelo agregado y un modelo distribuido. El modelo distribuido solo puede utilizarse si se activa el modelo de Green-Ampt para calcular la infiltración, ya que hereda sus parámetros de dicho modelo. El modelo agregado puede utilizarse con cualquiera de las formulaciones de infiltración definidas en la sección <span id='cite-_Ref107216854'></span>[[#_Ref107216854|4.1.4]].

En todo caso, es necesario tener activado algún modelo de infiltración para poder calcular el flujo subsuperficial. La elección del modelo de flujo subsuperficial se realiza en el menú '''Datos >> Datos del Problema >> Hidrología'''.

Si se selecciona el modelo de infiltración de Green-Ampt, en la ventana '''Datos >> Datos del Problema >> Hidrología '''aparecerá una nueva opción llamada '''Flujo Subsuperficial''' que permite elegir entre el modelo agregado o el distribuido (o desactivar el flujo subsuperficial). Si se selecciona otro modelo de infiltración (SCS, Lineal o Horton) aparecerá una opción llamada ''Groundwater flow'' que únicamente permitirá seleccionar el modelo agregado (o desactivar el flujo subsuperficial).

Una vez seleccionado un modelo de flujo subsuperficial, la asignación de los parámetros correspondientes a cada modelo se realiza a través de los menús '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo agregado''' y '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo distribuido'''.

'''''Modelo agregado de flujo subsuperficial'''''

Para aplicar el modelo agregado el usuario debe definir la extensión espacial de cada depósito subterráneo, los parámetros de cada depósito (constantes <math display="inline">A</math> y <math display="inline">B</math>, profundidad máxima y saturación inicial), y la zona de exfiltración de cada depósito.

La creación de depósitos y la definición de sus parámetros se realiza a través del menú '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo agregado >> Definición de depósitos''', desde donde se pueden crear tantos depósitos como queramos con parámetros diferentes. Las superficies de infiltración y exfiltración asignadas a dichos depósitos pueden asignarse a través de '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo agregado >> Asignación de depósitos''',''' '''haciendo referencia a los depósitos creados anteriormente. Es necesario aclarar que se debe asignar tanto las áreas de infiltración, como de exfiltración correspondientes a cada depósito.

'''''Modelo distribuido de flujo subsuperficial'''''

El modelo distribuido de flujo subsuperficial hereda la mayoría de sus parámetros del modelo de Green-Ampt. Únicamente es necesario definir adicionalmente la permeabilidad saturada horizontal del suelo, lo cual puede hacerse desde '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo distribuido >> Parámetros '''(<span id='cite-_Ref106615520'></span>[[#_Ref106615520|Fig. 30]]a)'''. '''De esta forma se puede asignar de forma manual una permeabilidad saturada a diferentes zonas del modelo (i.e. distribuida espacialmente).

Asimismo, es necesario definir los contornos por donde puede entrar o salir el agua de la zona de estudio a través del subsuelo. Si no se define ninguna condición de contorno para el modelo distribuido el agua no podrá salir del dominio de estudio a través del subsuelo (únicamente a través de la superficie). En dichos contornos se puede imponer como condición la descarga unitaria o la carga hidráulica. Dichas condiciones de contorno se definen a través del menú '''Datos >> Procesos hidrológicos >> Flujo subsuperficial >> Modelo distribuido >> Condiciones de contorno, '''siendo necesario asignarlas a los contornos del modelo correspondientes (<span id='cite-_Ref106615520'></span>[[#_Ref106615520|Fig. 30]]b).

<span id='_Ref106615520'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image30.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 30. (a) Ventana de asignación de los parámetros del modelo distribuido. (b) Ventana de asignación de condiciones de contorno del modelo distribuido.'''</span>


Las condiciones iniciales de humedad del subsuelo se heredan de los parámetros definidos en el modelo de Green-Ampt. Concretamente se utiliza la saturación inicial y la porosidad definidas en Green-Ampt para determinar la humedad inicial del subsuelo.

<span id='_Toc108335472'></span>

====4.2.6 Secciones de aforo y series temporales====

Iber permite definir secciones de aforo y sondas temporales, a través de las cuales obtener salidas de resultados de forma localizada. En las secciones de aforo, Iber calcula el caudal que pasa a través de ellas, generando una salida de hidrogramas en cada sección de aforo definida antes del cálculo.

Las secciones de aforo se definen a través del menú '''Datos >> Hidrodinámica >> Condiciones Internas. '''En esta ventana, en la opción '''Tipo de condición''' deberemos elegir '''Estación de aforo '''(<span id='cite-_Ref106616124'></span>[[#_Ref106616124|Fig. 31]]a). Es necesario asignarle un número de estación diferente a cada estación de aforo.

Los resultados se escriben en forma de series temporales en el fichero ''Q_hydrograph.rep'' (<span id='cite-_Ref106616124'></span>[[#_Ref106616124|Fig. 31]]b). Este archivo que también incluye los hidrogramas correspondientes a otras condiciones internas que se hayan definido en el modelo, como puede ser el caso de puentes, vertederos y/o compuertas.

La definición de sondas temporales se realiza desde el menú '''Datos >> Datos del Problema >> Resultados >> Sondas'''. Desde esta ventana el usuario puede introducir las coordenadas de las sondas, así como el intervalo temporal de salida de resultados (que no tiene por qué coincidir con el intervalo de salida temporal de resultados en toda la malla de cálculo). En las sondas definidas se almacenan las variables hidráulicas interpoladas linealmente a partir de su valor en las celdas de la malla de cálculo. Las series temporales resultantes se escriben en los ficheros ''res_gauge_NN.rep'', donde NN es el número de la sonda.

<span id='_Ref106616124'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image31.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 31. (a) Ventana de la función de asignación de secciones de aforo. (b) Ejemplo de formato de archivo Q_hydrograph.rep.'''</span>


<span id='_Toc108335473'></span>
<span id='_Ref106263333'></span>

==5 Módulo distribuido de erosión de suelos en laderas por escorrentía==

<span id='_Toc108335474'></span>

===5.1 Fundamentos teóricos===

El módulo de erosión de suelos es un modelo distribuido y bidimensional, que resuelve una serie de ecuaciones con base física dependientes de variables utilizadas habitualmente en hidrología. La ecuación principal es la de conservación de la masa de sedimento en el suelo, la cual se resuelve en cada elemento de la malla de cálculo y para cada clase de sedimento que se defina (en el caso de granulometrías distribuidas). Además, se resuelve la ecuación de conservación de sedimento en suspensión, también para cada clase de sedimento. Ambas ecuaciones están relacionadas por un término que modela la sedimentación y resuspensión de partículas de sedimento.

<span id='_Toc108335475'></span>

====5.1.1 Procesos considerados====

En la <span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]] se representa de forma esquemática la conceptualización de la estructura del suelo que se realiza en el módulo de erosión de suelos.  La capa superior representa el agua de escorrentía (o corriente de un río) que transporta sedimento en suspensión, la cual se extiende entre la cota del terreno ( <math display="inline">{z}_{b}</math>) y la superficie de la lámina de agua ( <math display="inline">{z}_{s}</math>), con un calado <math display="inline">h=</math><math>{z}_{s}-{z}_{b}</math>. La concentración de sedimentos en suspensión en dicha capa (promediada en profundidad) es igual a <math display="inline">C</math>, y en ella existe un transporte de sedimentos en suspensión igual <math display="inline">{q}_{s}=</math><math>C\, h\, U</math>.

Por debajo de la capa de escorrentía existe una primera capa de suelo formada por sedimento no-cohesivo que ya ha sido erosionada previamente de la matriz original de suelo. Las partículas de sedimento de dicha capa pueden ponerse en suspensión y pasar a la capa superior, ya sea por el impacto de las gotas de lluvia o por la fricción de fondo ejercida por la escorrentía superficial. Estos dos procesos se representan en la <span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]] por los términos <math display="inline">{D}_{rdrd}</math> y <math display="inline">{D}_{fdre}</math>, respectivamente. Al mismo tiempo, una parte de las partículas de suelo que se encuentran en suspensión pueden depositarse en la capa de suelo no-cohesivo debido a su peso. Este proceso se representa en la <span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]] por el término <math display="inline">{D}_{dep}</math>.

Debajo de la capa de suelo no-cohesivo puede existir una segunda capa de suelo correspondiente a la matriz original de suelo que aún no ha sido erosionado. Suele corresponder a un suelo con fuerzas de cohesión importantes, por ello nos referiremos a dicha capa como suelo cohesivo, por contraposición a la capa de suelo no-cohesivo. El impacto de las gotas de lluvia y la fuerza de fricción de fondo generada por la escorrentía superficial pueden poner en suspensión el sedimento de la capa cohesiva. Estos procesos se representan en la <span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]] por los términos <math display="inline">{D}_{rdd}</math> y <math display="inline">{D}_{fde}</math> respectivamente, y son análogos a los términos <math display="inline">{D}_{rdrd}</math> y <math display="inline">{D}_{fdre}</math>.

<span id='_Ref106809760'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image32.png|450px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 32. Esquema de la conceptualización de la estructura del suelo considerada en el módulo de erosión de suelos.'''</span>


Nótese en la <span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]] que el proceso de deposición ( <math display="inline">{D}_{dep}</math>) siempre aporta sedimento a la capa no-cohesiva, ya que no existen fuerzas de cohesión entre las partículas en suspensión que se depositan.

En cada una de las dos capas de suelo pueden existir distintas clases de sedimento, en general diferenciadas entre ellas por su diámetro. Esto permite considerar granulometrías distribuidas o diferentes tipos de suelo. La proporción de cada clase de sedimento en la capa de suelo erosionada se denotará con el factor <math display="inline">{f}_{k}</math>, que representa el tanto por uno de la clase <math display="inline">k</math> en la capa no-cohesiva. Al representar la proporción de cada clase con respecto al sedimento total, la suma de estos factores debe ser igual a 1:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\sum _{k=1}^{N}{f}_{k}=1\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(69)
|}

siendo <math display="inline">N</math> el número de clases de sedimento consideradas. La distribución de clases de sedimento en la capa erosionada (factores <math display="inline">{f}_{k}</math>) se actualiza durante el cálculo en función del sedimento erosionado y depositado en dicha capa.

<span id='_Toc54964599'></span><span id='_Toc90987664'></span><span id='_Toc108335476'></span>

====5.1.2 Definición de variables utilizadas en varias formulaciones====

Las siguientes variables se utilizan en varias de las formulaciones implementadas en el modelo:

:* Densidad relativa sumergida de las partículas del sedimento:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>R=\frac{{\rho }_{s}-\rho }{\rho }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(70)
|}

:* Diámetro adimensional de las partículas de la clase <math display="inline">k</math>:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{\ast ,k}={D}_{k}{\left( \frac{Rg}{{\nu }^{2}}\right) }^{0,33}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(71)
|}

:* Velocidad de fricción de fondo:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{u}_{\ast }=\sqrt{\frac{{\tau }_{b}}{\rho }}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(72)
|}

:* Tensión de fondo adimensional:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{k}^{\ast }=\frac{{u}_{\ast }^{2}}{Rg{D}_{k}}=\frac{\tau }{\rho Rg{D}_{k}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(73)
|}

:* Velocidad de sedimentación de las partículas según van Rijn [37]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{D}_{k}<{10}^{-4}m&\rightarrow &{w}_{s,k}=\frac{Rg{D}_{k}^{2}}{18\, \nu }\\{10}^{-4}m<{D}_{k}<{10}^{-3}m&\rightarrow &{w}_{s,k}=\frac{10\, \nu }{{D}_{k}}\left( \sqrt{1+0,01{D}_{\ast ,k}^{3}}-1\right) \\{D}_{k}>{10}^{-3}m&\rightarrow &{w}_{s,k}=1,1\, \sqrt{Rg{D}_{k}}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(74)
|}

siendo <math display="inline">\nu</math>  la viscosidad cinemática del agua ( <math display="inline">{10}^{-6}</math> m<sup>2</sup>/s).

<span id='_Toc108335477'></span>

====5.1.3 Ecuaciones de la hidrodinámica====

Para definir los campos de profundidad y velocidad se utilizan las ecuaciones de aguas poco profundas bidimensionales, incluyendo los términos que representan las componentes de precipitación e infiltración (Ecuación <span id='cite-_Ref106269045'></span>[[#_Ref106269045|(1]])). La simulación de la fricción de fondo se realiza mediante la inclusión de la ecuación de Manning.

<span id='_Ref107212053'></span><span id='_Toc108335478'></span>

====5.1.4 Transporte por carga en suspensión====

'''''Conservación de la masa de sedimento en suspensión'''''

El cálculo de la concentración de sedimento en suspensión se realiza resolviendo la siguiente ecuación de transporte en suspensión promediada en profundidad, para cada clase ( <math display="inline">k</math>) de sedimento considerada en el modelo:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{\partial h{C}_{k}}{\partial t}+\left( \frac{\partial {q}_{x}\, {C}_{k}}{\partial x}+\right. </math><math>\left. \frac{\partial {q}_{y}\, {C}_{k}}{\partial y}\right) ={D}_{rdd,k}+{D}_{rdrd,k}+</math><math>{D}_{fde,k}+{D}_{fdre,k}-{D}_{dep,k}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(75)
|}

donde <math display="inline">{C}_{k}</math> [kg/m<sup>3</sup>] es la concentración en suspensión de la clase de sedimento <math display="inline">k</math>. Es decir, se resuelve una ecuación de transporte en suspensión para cada clase de sedimento. Los términos a la derecha de la igualdad representan cada uno de los procesos de puesta en suspensión y deposición considerados en el modelo (<span id='cite-_Ref106809760'></span>[[#_Ref106809760|Fig. 32]]). Sus unidades son [kg/m<sup>2</sup>/s]:

:* <math display="inline">{D}_{rdd}</math> representa la puesta en suspensión de partículas de suelo procedentes de la capa cohesiva debido al impacto de las gotas de lluvia (''rainfall-driven detachment'').
:* <math display="inline">{D}_{rdrd}</math> representa la puesta en suspensión de partículas de suelo procedentes de la capa no-cohesiva debido al impacto de las gotas de lluvia (''rainfall-driven re-detachment'').
:* <math display="inline">{D}_{fde}</math> representa la puesta en suspensión de partículas de suelo procedentes de la capa cohesiva debido a la fricción de fondo ejercida por la escorrentía superficial (''flow-driven entrainment'').
:* <math display="inline">{D}_{fdre}</math> representa la puesta en suspensión de partículas de suelo procedentes de la capa no-cohesiva debido a la fricción de fondo ejercida por la escorrentía superficial (''flow-driven re-entrainment'').
:* <math display="inline">{D}_{dep}</math> representa la deposición de partículas en suspensión en la capa no-cohesiva, debido a su velocidad de sedimentación.

A continuación, se detallan las formulaciones empleadas para calcular los términos anteriores.

'''''Producción de sedimento en suspensión por impacto de lluvia'''''

Los términos de puesta en suspensión de sedimento debido al impacto de las gotas de lluvia ( <math display="inline">{D}_{rdd}</math> y <math display="inline">{D}_{rdrd}</math>) expresados en [kg/m<sup>2</sup>/s] se calculan como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{D}_{rdd,k}={\alpha }_{d,k}\, R\, {\sigma }_{h}\left( 1-\, \epsilon \right) &\rightarrow &{D}_{rdrd,k}={\alpha }_{rd,k}\, R\, {\sigma }_{h}\, \epsilon \, {f}_{k}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(76)
|}

donde <math display="inline">R</math> es la intensidad de lluvia [m/s], <math display="inline">{\alpha }_{d,k}</math> y <math display="inline">{\alpha }_{rd,k}</math> [kg/m<sup>3</sup>] son los coeficientes de erosionabilidad para la capa cohesiva y no-cohesiva respectivamente (''detachment'' y ''redetachment''), <math display="inline">{\sigma }_{h}</math> es un factor que considera el efecto protector de la lámina de escorrentía sobre el suelo, y <math display="inline">\epsilon</math>  es un factor adimensional de escudo que representa el efecto protector que la capa ya erosionada (no-cohesiva) ejerce sobre la matriz original de suelo (capa cohesiva). El término <math display="inline">{D}_{rdrd,k}</math> incluye la proporción de la clase de sedimento <math display="inline">k</math> en la capa erosionada. De esta forma, si a lo largo del cálculo en algún momento se erosiona toda una clase de sedimento, su proporción será igual a cero ( <math display="inline">{f}_{k}=</math><math>0</math>) y, por lo tanto, dicha clase no contribuirá a generar más sedimento en suspensión ( <math display="inline">{D}_{rdrd,k}=</math><math>0</math>).

En las ecuaciones anteriores, también se puede expresar la lluvia en [mm/h] y los coeficientes de erosionabilidad en [(kg/s/m<sup>2</sup>)/(mm/h)]. De esta forma se puede interpretar el coeficiente de erosionabilidad como los kg de sedimento erosionado por unidad de tiempo (segundo), por unidad de superficie (m<sup>2</sup>) y por unidad de intensidad de precipitación (mm/h).

El factor <math display="inline">{\sigma }_{h}</math> controla el efecto protector de la lámina de escorrentía sobre la erosión del suelo. Si el calado de la lámina de escorrentía es elevado, las gotas de lluvia no serán capaces de poner el sedimento en suspensión ya que su energía será disipada por la lámina de agua. El término <math display="inline">{\sigma }_{h}</math> se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\sigma }_{h}={\left( \frac{{z}_{m}}{h}\right) }^{0,8}\quad \, \, {\sigma }_{h}\in \left( 0,1\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(77)
|}

siendo <math display="inline">{z}_{m}</math> un calado umbral por debajo del cual no existe efecto protector, que debe ser definido como parámetro de entrada del modelo. En efecto, si el calado de la escorrentía es inferior a dicho calado umbral ( <math display="inline">h<{z}_{m}</math>) el valor de <math display="inline">{\sigma }_{h}=</math><math>1</math>. A medida que el calado aumenta por encima de <math display="inline">{z}_{m}</math>, el valor de <math display="inline">{\sigma }_{h}</math> comienza a disminuir, tendiendo a cero para calados elevados.

Por su parte, el coeficiente protector <math display="inline">\epsilon</math> , cuyo valor está comprendido entre 0 y 1, se puede interpretar como un factor de reparto de la energía de la lluvia entre la capa cohesiva y la no cohesiva. Se calcula considerando una relación lineal en base a la masa de sedimento por unidad de superficie en la capa erosionada ( <math display="inline">{M}_{s}</math>) y a un valor crítico de dicha masa de sedimento ( <math display="inline">{M}_{s,cr}</math>), expresados ambos en [kg/m<sup>2</sup>]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\, \epsilon =min\left[ \frac{{M}_{s}}{{M}_{s,cr}},1\right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(78)
|}

Cuando <math display="inline">{M}_{s}</math> alcanza o supera el valor crítico ( <math display="inline">{M}_{s}\geq {M}_{s,cr}</math>), el efecto protector es máximo ( <math display="inline">\epsilon =</math><math>1</math>), siendo nulo el material erosionado en la capa cohesiva ( <math display="inline">{D}_{rdd}=</math><math>0</math>). Cuando no existe material en la capa no-cohesiva ( <math display="inline">{M}_{s}=</math><math>0)</math> el efecto protector es nulo ( <math display="inline">\epsilon =0</math>).

'''''Producción de sedimento en suspensión por escorrentía superficial (fricción de fondo)'''''

El término de puesta en suspensión desde la capa cohesiva por efecto de la fricción de fondo ( <math display="inline">{D}_{fde}</math>) se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{fde,k}=\left\{ \begin{matrix}{K}_{k}\, \left( \tau -{\tau }_{crc}\right) \left( 1-\epsilon \right) &si&\, {\tau >\tau }_{cr}\\0&si&{\tau \leq \tau }_{cr}\end{matrix}\right.</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(79)
|}

donde <math display="inline">\tau</math>  en [N/m<sup>2</sup>] es la tensión de fondo, <math display="inline">{\tau }_{crc}</math> en [N/m<sup>2</sup>] es la tensión crítica de la matriz original de suelo (cohesiva), por debajo del cual no se produce puesta en suspensión desde la capa cohesiva por fricción de fondo y <math display="inline">{K}_{k}</math> en [kg/s/N] es un coeficiente a definir por el usuario (''flow-driven detachability coefficient''). Nótese que el valor de <math display="inline">{\tau }_{crc}</math> es el mismo para todas las clases de sedimento, ya que mientras que la tensión de fondo no supere dicho valor no se producirá la erosión de las partículas de dicha capa. Una vez que se supera la tensión crítica de la capa cohesiva ( <math display="inline">\tau >{\tau }_{crc}</math>), la contribución de cada clase de sedimento puede ser diferente, por lo que el coeficiente <math display="inline">{K}_{k}</math> puede ser diferente para cada clase.

Respecto al flujo de partículas desde la capa no-cohesiva a la columna de agua, su cálculo está ligado a los términos <math display="inline">{D}_{fdre}</math> y <math display="inline">{D}_{dep}</math>, siendo el flujo neto de sedimento entre la capa de suelo erosionada y la columna de agua igual a la diferencia entre dichos términos. En Iber se contemplan dos formulaciones diferentes, siendo el usuario el que debe elegir cuál utilizar: Hairsine y Rose o van Rijn.

:* Hairsine y Rose

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{fdre,k}=\left\{ \begin{matrix}\frac{{\rho }_{s}{r}_{f}{F}_{k}}{\left( {\rho }_{s}-\rho \right) g}\left( \frac{{\Omega -\Omega }_{cr}}{h}\right) \epsilon \, {f}_{k}&si&\, {\Omega >\Omega }_{cr}\\0&si&{\Omega \leq \Omega }_{cr}\end{matrix}\right.</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(80)
|}

donde <math display="inline">\Omega</math>  en [W/m<sup>2</sup>] es la potencia disipada por fricción de fondo y por unidad de superficie, <math display="inline">{\Omega }_{cr}</math> en [W/m<sup>2</sup>] es un valor crítico de <math display="inline">\Omega</math> , por debajo del cual no se produce puesta en suspensión desde la capa no-cohesiva por fricción de fondo, <math display="inline">{F}_{k}</math> es un coeficiente adimensional que expresa el tanto por uno de potencia disipada (en exceso sobre <math display="inline">{\Omega }_{cr}</math>) que contribuye a la puesta en suspensión de sedimento desde la capa erosionada, <math display="inline">{\rho }_{s}</math> es la densidad de las partículas de sedimento (en torno a 2400 kg/m<sup>3</sup>, generalmente), <math display="inline">\rho</math>  el la densidad del agua (1000 kg/m<sup>3</sup>), <math display="inline">h</math> [m] es el calado, <math display="inline">g</math> es la aceleración de la gravedad (9,81 m/s<sup>2</sup>), y <math display="inline">{r}_{f}</math> es el factor de Rouse que determina si el perfil de concentración en profundidad es más o menos homogéneo. Para flujos de escorrentía con calados muy bajos, se puede asumir un factor de Rouse igual a 1 ( <math display="inline">{r}_{f}=</math><math>1</math>).

En caso de utilizar el modelo de Hairsine y Rose, el término de deposición se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{dep,k}={r}_{f}{w}_{e,k}{C}_{k}\quad \quad {w}_{e,k}=\mathrm{max}\,\left( {w}_{s,k}-{u}_{\ast }\, ,\, 0\right) \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(81)
|}

siendo <math display="inline">{w}_{s,k}</math> en [m/s] la velocidad de sedimentación de las partículas en suspensión y <math display="inline">{u}_{\ast }</math> la velocidad de fricción de fondo. La velocidad efectiva de sedimentación <math display="inline">{w}_{e,k}</math> tiene en cuenta el efecto de la turbulencia generada por fricción de fondo, que tiende a mantener en flotación a las partículas de sedimento. La velocidad de sedimentación de las partículas de sedimento se obtiene de su densidad y diámetro mediante la formulación de van Rijn [37].

El factor de Rouse se toma igual a 1 ( <math display="inline">{r}_{f}=1</math>).

:* van Rijn

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{fdre,k}={\rho }_{s}\, {C}_{k}^{\ast }\, {w}_{s,k}\, \epsilon \, {f}_{k}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(82)
|}

donde <math display="inline">{\rho }_{s}</math> es la densidad de las partículas de sedimento, <math display="inline">{w}_{s,k}</math> en [m/s] es la velocidad de sedimentación de las partículas en suspensión calculada a partir de su diámetro y densidad mediante la formulación de van Rijn, y <math display="inline">{C}_{k}^{\ast }</math> es una concentración de equilibrio, la cual se puede calcular mediante alguna de las formulaciones siguientes:

::* van Rijn [38]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{C}_{k}^{\ast }=\frac{0,005\, {T}_{k}^{1,5}}{{D}_{\ast ,k}^{0,3}}\quad \quad {T}_{k}=</math><math>\frac{{\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\, }{{\tau }_{c,k}^{\ast }}\quad \, \quad {\tau }_{c,k}^{\ast }=</math><math>\frac{{w}_{s,k}^{2}}{Rg{D}_{k}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(83)
|}

::* Smith [39]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{C}_{k}^{\ast }=\frac{0,00156\, {T}_{k}}{1+0,0024{T}_{k}}\quad \quad {T}_{k}=\frac{{\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\, }{{\tau }_{c,k}^{\ast }}\, \, \quad \, \, {\tau }_{c,k}^{\ast }=</math><math>\frac{{w}_{s,k}^{2}}{Rg{D}_{k}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(84)
|}

::* García y Parker [40]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{C}_{k}^{\ast }=\frac{1,3\cdot {10}^{-7}\, {Z}_{k}^{5}}{1+4,33\cdot {10}^{-7}\, {Z}_{k}^{5}}\quad \quad {Z}_{k}=</math><math>\frac{{u}_{\ast }}{{w}_{s,k}}{D}_{\ast ,k}^{0,9}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(85)
|}

::* Zyserman y Fredsoe [41]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{C}_{k}^{\ast }=\frac{0,331\, {\left( {\tau }_{k}^{\ast }-0,045\right) }^{1,75}}{1+0,72\, {\left( {\tau }_{k}^{\ast }-0,045\right) }^{1,75}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(86)
|}

En caso de utilizar el modelo de van Rijn, el término de deposición se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{D}_{dep,k}={r}_{f}{w}_{s,k}{C}_{k}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(87)
|}

siendo <math display="inline">{w}_{s,k}</math> en [m/s] la velocidad de sedimentación de las partículas en suspensión calculada mediante la formulación de van Rijn [37]. El factor de Rouse se toma igual a 1 ( <math display="inline">{r}_{f}=</math><math>1</math>).

<span id='_Ref107215785'></span><span id='_Toc108335479'></span>

====5.1.5 Transporte por carga de fondo====

El transporte por carga de fondo representa el movimiento de partículas de sedimento que se arrastran por la superficie del terreno como consecuencia de la fuerza de fricción ejercida por el agua sobre ellas. No son partículas que se transporten en suspensión en el seno de la columna de agua. No debe confundirse el transporte por carga de fondo con la puesta en suspensión de partículas por fricción de fondo, representado por los términos <math display="inline">{D}_{fde}</math> y <math display="inline">{D}_{fdre}</math> definidos anteriormente. En efecto, la fricción de fondo puede poner en suspensión ciertas partículas (las más finas y menos pesadas) y posteriormente estas partículas transportarse ''flotando'' en el seno de la columna de agua. En este caso estaríamos hablando de transporte en suspensión, y se calcularía con las formulaciones detalladas en el apartado <span id='cite-_Ref107212053'></span>[[#_Ref107212053|5.1.4]]. Pero cuando la fricción de fondo no tiene suficiente energía para poner en suspensión las partículas más pesadas, éstas se arrastrarán por la superficie del terreno. En este caso el transporte sería por carga de fondo.

El transporte por carga de fondo no suele considerarse en los modelos de erosión por lluvia-escorrentía, ya que generalmente se asume que este tipo de fenómenos producen básicamente el transporte de sedimento muy fino y ligero, que por lo tanto se transporta en suspensión. La mayor parte del suelo erosionado por lluvia-escorrentía son limos o arcillas. En el caso de que exista cierto transporte de fondo de las partículas más gruesas (arenas), la masa de suelo erosionada de estas fracciones de suelo suele ser mucho menor en magnitud que la erosionada por transporte en suspensión.

De todas formas, en el modelo de erosión implementado en Iber, sí que se ha considerado la componente de transporte de fondo, la cual puede ser activada o desactivada por el usuario. De esta manera se consigue una mayor versatilidad en el rango de aplicaciones del modelo, así como su posible aplicación a estudios de investigación en los que sea necesario cuantificar el transporte de fondo.

'''''Condiciones de transporte de fondo de equilibrio y no-equilibrio'''''

Para calcular el transporte de partículas de sedimento por carga de fondo se contemplan dos opciones: una formulación de equilibrio y una formulación de no-equilibrio. El usuario debe elegir qué formulación utilizar, siendo sensiblemente más sencilla la formulación de equilibrio.

:* Formulación de equilibrio

En el caso de utilizar la formulación de equilibrio, se asume que las formulaciones de capacidad de transporte de fondo expuestas en el apartado siguiente son válidas para condiciones no-estacionarias (de no-equilibrio). En este caso, se asume que en todo momento el caudal sólido de fondo real es igual a la capacidad de transporte de fondo, simplificándose la resolución del modelo.

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{b,k}={q}_{pot,k}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(88)
|}

donde <math display="inline">{q}_{b,k}</math> es el transporte de fondo de las partículas de la clase <math display="inline">k</math>, y <math display="inline">{q}_{pot,k}</math> es la capacidad de transporte de fondo de la clase de sedimento <math display="inline">k</math>.

:* Formulación de no-equilibrio

En el caso de utilizar la formulación de no-equilibrio, el caudal sólido de fondo se calcula para cada clase de sedimento ( <math display="inline">{q}_{b,k}</math>) mediante la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>-\frac{\partial }{\partial x}\left( {q}_{bx,k}\right) -\frac{\partial }{\partial y}\left( {q}_{by,k}\right) +</math><math>\frac{{q}_{b,k}-{q}_{pot,k}}{{L}_{s}}=0</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(89)
|}

siendo,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{L}_{s}=3{D}_{50}{D}_{\ast }^{0,6}{T}^{0,9}\quad \quad {D}_{\ast }={D}_{50}{\left( \frac{gR}{{\nu }^{2}}\right) }^{\frac{1}{3}}\quad \, \, T=</math><math>\frac{{\tau }^{\ast }-{\tau }_{c}^{\ast }}{{\tau }_{c}^{\ast }}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(90)
|}

donde <math display="inline">{q}_{bx},{q}_{by}</math> son las proyecciones en los ejes X e Y del caudal sólido por transporte de fondo, y <math display="inline">{L}_{s}</math> es la longitud de adaptación a las condiciones de equilibrio. Un valor pequeño de <math display="inline">{L}_{s}</math> implica que el transporte de fondo se adapta muy rápidamente a los cambios en las condiciones de flujo, en cuyo caso <math display="inline">{q}_{b,k}=</math><math>{q}_{pot,k}</math> y la capacidad de transporte de fondo depende únicamente de las condiciones de flujo locales. Al contrario, un valor grande de <math display="inline">{L}_{s}</math> implica que la adaptación de la carga de fondo a las condiciones de flujo es relativamente lenta, evitándose transiciones bruscas en la capacidad de transporte de fondo de un punto a otro del dominio. Como se ha comentado, la utilización de la formulación de no-equilibrio complica de forma considerable el cálculo del transporte por carga de fondo.

<span id='_Ref38265327'></span><span id='_Toc54964607'></span><span id='_Toc90987672'></span>
'''''Capacidad de transporte de fondo. Formulaciones'''''

Tanto si se utiliza la formulación de equilibrio como la de no-equilibrio, es necesario calcular la capacidad de transporte de fondo en condiciones de equilibrio ( <math display="inline">{q}_{pot,k}</math>) para cada clase de sedimento. Para ello, se implementan las siguientes formulaciones, siendo el usuario el que debe elegir cuál utilizar:

:* Genérica

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{pot,k}^{\ast }=A\, {\left( {\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\right) }^{B}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(91)
|}

donde <math display="inline">A</math> y <math display="inline">B</math> son parámetros impuestos por el usuario y <math display="inline">{\tau }_{c,k}^{\ast }</math> es la tensión crítica adimensional, que puede ser calculada mediante el ábaco de Shields o la formulación de Parker et al. [42].

:* Meyer-Peter y Müller [43]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{pot,k}^{\ast }=8\, {\left( {\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\right) }^{1,5}\quad \quad {\tau }_{c,k}^{\ast }=</math><math>0,047</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(92)
|}

donde <math display="inline">{\tau }_{c,k}^{\ast }=0,047</math> es la tensión crítica adimensional.

:* Wong y Parker [44]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{pot,k}^{\ast }=3,97\, {\left( {\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\right) }^{1,5}\quad \quad {\tau }_{c,k}^{\ast }=</math><math>0,0495</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(93)
|}

donde <math display="inline">{\tau }_{c,k}^{\ast }=0,0495</math> es la tensión crítica adimensional.

:* Einstein y Brown [45]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{\tau }_{k}^{\ast }<0,18&{q}_{pot,k}^{\ast }=2,15\mathrm{exp}\,\left( \frac{-0,391}{{\tau }_{k}^{\ast }}\right) \\0,18<{\tau }_{k}^{\ast }<0,52&{q}_{pot,k}^{\ast }=40\, {\left( {\tau }_{k}^{\ast }\right) }^{3}\\0,52<{\tau }_{k}^{\ast }&{q}_{pot,k}^{\ast }=15\, {\left( {\tau }_{k}^{\ast }\right) }^{1,5}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(94)
|}

En la formulación de Einstein-Brown la tensión crítica adimensional es igual a cero ( <math display="inline">{\tau }_{c,k}^{\ast }=</math><math>0)</math>.

:* van Rijn [37]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{T}_{k}<0,3&{q}_{pot,k}^{\ast }=0,053\frac{{\left( {T}_{k}\right) }^{2,1}}{{D}_{\ast ,k}^{0,3}}&\\&&siendo\, {T}_{k}=\frac{{\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\, }{{\tau }_{c,k}^{\ast }}\\{T}_{k}>0,3&{q}_{pot,k}^{\ast }=0,100\frac{{\left( {T}_{k}\right) }^{1,5}}{{D}_{\ast ,k}^{0,3}}&\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(95)
|}

donde <math display="inline">{\tau }_{c,k}^{\ast }</math> es la tensión crítica adimensional, a definir por el usuario según alguna de las formulaciones expuestas al final de esta sección.

:* Engelund y Hansen [45]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{pot,k}^{\ast }=\frac{0,1\, {\left( {\tau }_{k}^{\ast }\right) }^{2,5}}{\frac{2\, {u}_{\ast }^{2}\, {D}_{k}}{h\, {V}^{2}}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(96)
|}

En la formulación de Engelund-Hansen la tensión crítica adimensional es igual a cero ( <math display="inline">{\tau }_{c,k}^{\ast }=</math><math>0)</math>.

:* Ashida y Michiue [45]

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{pot,k}^{\ast }=17\left( {\tau }_{k}^{\ast }-{\tau }_{c,k}^{\ast }\right) \left( \sqrt{{\tau }_{k}^{\ast }}-\right. </math><math>\left. \sqrt{{\tau }_{c,k}^{\ast }}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(97)
|}

donde <math display="inline">{\tau }_{c,k}^{\ast }</math> es la tensión crítica adimensional, a definir por el usuario según alguna de las formulaciones expuestas al final de esta sección.

'''''Tensión crítica para transporte de fondo'''''

Algunas de las formulaciones anteriores (Meyer-Peter y Müller, Wong y Parker, Einstein y Brown, Engelund y Hansen) tienen un valor de la tensión crítica de fondo bien definido en la propia formulación, mientras que en otras (Genérica, van Rijn, Ashida y Michiue) este valor debe ser definido por el usuario (<span id='cite-_Ref106296446'></span>[[#_Ref106296446|Tabla 4]]). En este último caso, el usuario puede elegir entre imponer un valor constante de la tensión crítica adimensional, calcularla mediante el ábaco de Shields o mediante la formulación de Parker et al. [42].

A estos efectos, el ábaco de Shields se implementa según la parametrización propuesta por Hoffmans y Verheij [46]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}0<{D}_{\ast }<4&{\tau }_{c}^{\ast }={D}_{\ast }^{-1}\\4<{D}_{\ast }<10&{\tau }_{c}^{\ast }={D}_{\ast }^{-1}\\10<{D}_{\ast }<20&{\tau }_{c}^{\ast }={0,04\, D}_{\ast }^{-0,10}\end{matrix}</math>

<math>\begin{matrix}20<{D}_{\ast }<150&{\tau }_{c}^{\ast }={0,013\, D}_{\ast }^{0,29}\\{D}_{\ast }>150&{\tau }_{c}^{\ast }=0,055\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(98)
|}

Alternativamente a la anterior parametrización del ábaco de Shields, puede utilizarse la siguiente formulación propuesta por Parker et al. para el cálculo de la tensión crítica, la cual proporciona valores de umbral del movimiento sensiblemente inferiores al ábaco de Shields original:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{c}^{\ast }=0,5\, \left( 0,22\, {R}_{p}^{-0,6}+0,06\cdot {10}^{\left( -7,7\, {R}_{p}^{-0,6}\right) }\right) \quad \, \, {R}_{p}=</math><math>\frac{\sqrt{Rg{D}^{3}}}{\nu }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(99)
|}

<span id='_Ref106296446'></span>
<span style="text-align: center; font-size: 75%;">'''Tabla 4. Tensión crítica adimensional para diferentes formulaciones de capacidad de transporte de fondo.'''</span>

{| style="width: 57%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|'''Formulación '''
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|'''Tensión crítica adimensional'''
|-
|  style="vertical-align: top;"|Genérica
|  style="vertical-align: top;"|Constante, Shields o Parker et al. [42]
|-
|  style="vertical-align: top;"|Meter-Peter y Müller
|  style="vertical-align: top;"|0,047
|-
|  style="vertical-align: top;"|Wong y Parker
|  style="vertical-align: top;"|0,0495
|-
|  style="vertical-align: top;"|Einstein y Brown
|  style="vertical-align: top;"|0,0
|-
|  style="vertical-align: top;"|van Rijn
|  style="vertical-align: top;"|Constante, Shields o Parker et al. [42]
|-
|  style="vertical-align: top;"|Engelund y Hansen
|  style="vertical-align: top;"|0,0
|-
|  style="border-bottom: 2pt solid black;vertical-align: top;"|Ashida y Michiue
|  style="border-bottom: 2pt solid black;vertical-align: top;"|Constante, Shields o Parker et al. [42]
|}


'''''Factor de escondimiento (hiding factor)'''''

En sedimentos con granulometría distribuida (multiclase), la tensión crítica adimensional ( <math display="inline">{\tau }_{c}^{\ast }</math>) se calcula para el diámetro medio del sedimento. A partir de dicho valor, se calcula la tensión crítica para cada clase de partículas como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{c,k}^{\ast }={\tau }_{c}^{\ast }{\left( \frac{{D}_{50}}{{D}_{k}}\right) }^{\gamma }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(100)
|}

siendo <math display="inline">\gamma</math>  un factor de ocultamiento (''hiding factor'') con un valor comprendido entre 0 y 1. Dicho factor represente la interacción entre partículas de diferente diámetro, de forma que las partículas gruesas protegen a las finas (aumentando su umbral de movimiento) y las finas exponen a las gruesas (facilitando su movilidad).

Teniendo en cuenta la definición de la tensión adimensional, la ecuación anterior se puede escribir como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{c,k}={\tau }_{c}{\left( \frac{{D}_{k}}{{D}_{50}}\right) }^{1-\gamma }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(101)
|}

De forma que si <math display="inline">\gamma =0</math> la interacción entre partículas es mínima (nula) y si <math display="inline">\gamma =</math><math>1</math> la interacción entre partículas es máxima:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}Si\, \gamma =0&{\tau }_{c,k}^{\ast }={\tau }_{c}^{\ast }&{\tau }_{c,k}={\tau }_{c}\frac{{D}_{k}}{{D}_{50}}\\Si\, \gamma =1&{\tau }_{c,k}^{\ast }={\tau }_{c}^{\ast }\frac{{D}_{50}}{{D}_{k}}&{\tau }_{c,k}={\tau }_{c}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(102)
|}

'''''Corrección del transporte de fondo por formas de fondo'''''

La tensión de fondo total en el lecho de un río está generada por la rugosidad de grano del sedimento (la cual es proporcional al diámetro del sedimento) y por las formas de fondo (rizos, dunas o antidunas). Únicamente la tensión por grano contribuye al movimiento de sedimentos por carga de fondo. Si bien este fenómeno es principalmente importante en ríos, se ha incluido como opción en el módulo de erosión, de forma que el usuario puede activarlo o desactivarlo a su voluntad.

Si se activa esta opción, previamente al cálculo del caudal sólido de fondo se estima la tensión de fondo debida a grano utilizando para ello la partición de tensiones de Einstein, en la cual se calcula la tensión de grano a partir de la tensión total como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{s}=\tau {\left( \frac{{n}_{s}}{n}\right) }^{1,5}\quad \, \, {n}_{s}=\frac{{\left( {K}_{s\, \left( m\right) }\right) }^{\frac{1}{6}}}{25}\quad \quad \, {K}_{s}=</math><math>\alpha \, {D}_{50}\quad \alpha \approx 2-3</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(103)
|}

siendo <math display="inline">n</math> el coeficiente de Manning total (utilizado en el cálculo de la hidrodinámica e introducido por el usuario), <math display="inline">{n}_{s}</math> el coeficiente de Manning equivalente debido a grano, <math display="inline">{D}_{50}</math> el diámetro medio del sedimento expresado en metros, <math display="inline">{K}_{s}</math> la altura de rugosidad de grano expresada en metros (calculada a partir del diámetro del sedimento), <math display="inline">\tau</math>  la tensión total de fondo, <math display="inline">{\tau }_{s}</math> la tensión de fondo debida a grano y <math display="inline">\alpha</math>  un factor que relaciona la altura de rugosidad con el diámetro medio del sedimento (''grain roughness factor''). El factor <math display="inline">\alpha</math>  lo debe introducir el usuario, si bien generalmente suelen utilizarse un valor en torno a 3.

Si se activa la corrección por formas de fondo, en el cálculo del transporte de fondo potencial ( <math display="inline">{q}_{pot,k}^{\ast }</math>) mediante las ecuaciones presentadas anteriormente, se utiliza la tensión de grano ( <math display="inline">{\tau }_{s}</math>) en lugar de la tensión total ( <math display="inline">\tau</math> ).

'''''Corrección del transporte de fondo por pendiente del terreno'''''

Las ecuaciones de transporte de fondo expuestas anteriormente son válidas únicamente para fondo horizontal. En el caso de que el fondo no sea horizontal es necesario introducir una corrección de pendiente que tenga en cuenta el efecto de la gravedad sobre la tensión crítica de movimiento y sobre el caudal sólido de fondo. El efecto de la gravedad puede ser positivo o negativo dependiendo de si esta actúa a favor o en contra de la velocidad del agua.

La corrección del caudal sólido de fondo por pendiente del terreno consiste, por un lado, en aumentar o minorar la tensión sobre el fondo del lecho (se mayora o minora dependiendo de si la fuerza de gravedad debida a la pendiente del fondo actúa a favor o en contra de la fuerza ejercida por el agua) y, por otro lado, en minorar la tensión crítica de movimiento (tensión máxima que puede soportar el sedimento sin que haya movimiento). De esta forma, se obtiene una tensión de fondo efectiva y una tensión crítica efectiva. En función de si la gravedad actúa a favor o en contra del movimiento de las partículas, la tensión efectiva sobre el fondo se mayora o se minora, mientras que la tensión crítica efectiva con fondo inclinado siempre es inferior a la tensión crítica con fondo horizontal.

La tensión crítica efectiva que es capaz de soportar el sedimento se evalúa como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{c,e}={\tau }_{c}\mathrm{cos}\,\alpha</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(104)
|}

siendo <math display="inline">\alpha</math>  el ángulo que forma el terreno con la horizontal, y <math display="inline">{\tau }_{c}</math> la tensión crítica de arrastre con fondo horizontal. Por otro lado, la tensión efectiva sobre el fondo se calcula como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{x,e}={\tau }_{x}-{\tau }_{c}\frac{\mathrm{sin}\,{\alpha }_{x}}{\mathrm{tan}\,\varphi }\quad \quad {\tau }_{y,e}=</math><math>{\tau }_{y}-{\tau }_{c}\frac{\mathrm{sin}\,{\alpha }_{y}}{\mathrm{tan}\,\varphi }\quad \, \, {\tau }_{e}=</math><math>\sqrt{{\tau }_{x,e}^{2}+{\tau }_{y,e}^{2}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(105)
|}

siendo <math display="inline">{\tau }_{x},{\tau }_{y}</math> las tensiones de arrastre fondo en las direcciones X e Y (tensión total o de grano en función de si se activa o no la corrección por formas de fondo), <math display="inline">{\tau }_{e}</math> la tensión de arrastre efectiva sobre el fondo, <math display="inline">{\alpha }_{x}</math>, <math display="inline">{\alpha }_{y}</math> las proyecciones sobre los ejes X e Y del ángulo que forma el fondo con la horizontal y <math display="inline">\varphi</math>  el ángulo de rozamiento interno del sedimento. Una vez calculadas la tensión crítica efectiva ( <math display="inline">{\tau }_{c,e}</math>) y la tensión de arrastre efectiva ( <math display="inline">{\tau }_{e}</math>), se utilizan estos valores en lugar de <math display="inline">{\tau }_{c}</math> y <math display="inline">\tau</math>  en las ecuaciones de transporte correspondientes.

<span id='_Toc54964612'></span><span id='_Toc90987677'></span><span id='_Toc108335480'></span>

====5.1.6 Evolución temporal de la masa erosionada====

Para calcular la evolución temporal de la masa de sedimento en la capa no-cohesiva (erosionada) se resuelve la siguiente ecuación de balance de masa en la capa erosionada, para cada clase de sedimento considerada:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{\partial {M}_{s,k}}{\partial t}=-{D}_{rdrd,k}-{D}_{fdre,k}+{D}_{dep,k}-{\rho }_{s}\left( \frac{\partial {q}_{bx,k}\, }{\partial x}-\right. </math><math>\left. \frac{\partial {q}_{by,k}\, }{\partial y}\right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(106)
|}

siendo <math display="inline">{M}_{s,k}</math> la masa de sedimento de la clase <math display="inline">k</math> por unidad de superficie, expresada en [kg/m<sup>2</sup>]. Tanto <math display="inline">{M}_{s}</math>, como los términos <math display="inline">{D}_{rdrd},\, {D}_{fdre}</math> y <math display="inline">{D}_{dep}</math>, son variables en espacio y en tiempo.

En cada paso de tiempo se actualiza el valor de <math display="inline">{M}_{s,k}</math> para cada clase de sedimento, y a partir de estos valores la masa total de sedimento en la capa erosionada ( <math display="inline">{M}_{s}</math>) y la nueva granulometría en dicha capa (factores <math display="inline">{f}_{k}</math>), como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{M}_{s}=\sum _{k}^{}{M}_{s,k}\quad \quad {f}_{k}=\frac{{M}_{s,k}}{{M}_{s}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(107)
|}

<span id='_Toc108335481'></span>

====5.1.7 Evolución temporal de la topografía====

La evolución temporal de la topografía debido a los procesos de erosión descritos anteriormente vendrá dada por la aplicación del principio de conservación de la masa para todas las clases de sedimento consideradas y en las dos capas de sedimento:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\left( 1-\phi \right) \frac{\partial {z}_{b}}{\partial t}=\sum _{k}^{}\left( -\left( \frac{{D}_{rdd,k}+{D}_{rdrd,k}+{D}_{fde,k}+{D}_{fdre,k}-{D}_{dep,k}}{{\rho }_{s}\, }\right) -\left( \frac{\partial {q}_{bx,k}\, }{\partial x}-\frac{\partial {q}_{by,k}\, }{\partial y}\right) \right)</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(108)
|}

siendo <math display="inline">\phi</math>  la porosidad del sedimento.

En cada paso de tiempo, la nueva topografía obtenida a partir de esta ecuación se utiliza en las ecuaciones hidrodinámicas para asegurar un correcto acoplamiento entre los procesos de erosión y el flujo de agua.

<span id='_Toc108335482'></span>

====5.1.8 Simplificaciones del modelo====

Las ecuaciones y formulaciones implementadas en el módulo de erosión de suelos, descritas en las secciones anteriores, son relativamente complejas e incluyen un elevado número de parámetros. Esto puede complicar la aplicación del modelo en casos reales en los cuales sea difícil definir todos los parámetros y opciones del modelo. Aun disponiendo datos observados para la calibración del modelo, puede ser complicado ajustar todos sus parámetros debido al alta no-linealidad de las formulaciones y a la posible equifinalidad entre parámetros [47]. Por último, puede ser interesante en algunas aplicaciones limitar la complejidad del modelo, por ejemplo, en su aplicación al lavado de sedimentos en entornos urbanos [48].

A continuación, se proponen algunas simplificaciones del modelo que pueden ser interesantes para su aplicación a casos reales.

'''''Una única capa de sedimento y lavado de finos'''''

La pérdida de suelo en cuencas rurales suele producirse principalmente por lavado de partículas finas del suelo. En este caso puede asumirse que el principal proceso de transporte de sólidos es por carga en suspensión, y despreciarse el transporte por carga de fondo ( <math display="inline">{q}_{b}=</math><math>0</math>), lo cual simplifica en gran medida el número de ecuaciones y parámetros del modelo.

Asimismo, puede simplificarse la estructura del suelo a una única capa (la capa erosionada), con el fin de evitar parametrizar dos capas de suelo con diferentes propiedades erosivas. En este caso los términos de producción de sedimento en suspensión <math display="inline">{D}_{rdd}</math> y <math display="inline">{D}_{fde}</math> son iguales a cero y no es necesario definir los parámetros correspondientes.

Si además:

:* asumimos disponibilidad ilimitada de sedimento en el suelo, definiendo un valor inicial de masa de sedimento ( <math display="inline">{M}_{s}</math>) elevado y un valor de masa crítica ( <math display="inline">{M}_{s,cr}</math>) reducido, lo cual implica <math display="inline">\epsilon =</math><math>1</math> siempre;
:* despreciamos el factor protector de la lámina de escorrentía imponiendo un valor de <math display="inline">{z}_{m}</math> muy elevado, lo cual implica <math display="inline">{\sigma }_{h}=</math><math>1</math>;
:* utilizamos la formulación de Harisine y Rose [49,50] para el lavado por escorrentía;

obtenemos un modelo con solamente 3 parámetros: <math display="inline">{\alpha }_{rd}</math>, <math display="inline">{\Omega }_{cr}</math> y <math display="inline">F</math>.

Podemos dar un paso más y despreciar la potencia crítica para la puesta en suspensión por escorrentía ( <math display="inline">{\Omega }_{cr}</math>), argumentando que su valor suele ser varios órdenes de magnitud inferiores a <math display="inline">\Omega</math> , con lo cual nos quedaría un modelo de únicamente 2 parámetros: <math display="inline">{\alpha }_{rd}</math> y <math display="inline">F</math>.

Este tipo de modelo fue utilizado en Cea et al. [47], proporcionando resultados satisfactorios para el cálculo del lavado de suelo en una parcela agrícola. Debe tenerse en cuenta que, en ciertas aplicaciones, es más interesante considerar la posible variabilidad espacial de estos coeficientes, que introducir parámetros adicionales que compliquen en exceso el modelo.

'''''Una única capa de sedimento y lavado de finos por impacto de lluvia'''''

Además de las simplificaciones realizadas anteriormente, en cuencas en las que no existe erosión en cárcavas (''rill erosion ''o'' gully erosion'') puede asumirse que el sedimento se pone en suspensión principalmente por acción de la lluvia, con lo cual podemos despreciar en primera aproximación el término <math display="inline">{D}_{fdre}</math>.

De esta forma nos quedaría un modelo con <math display="inline">{\alpha }_{rd}</math> como único parámetro (que podría ser variable espacialmente).

'''''Lavado de sólidos en zonas urbanas'''''

Este tipo de modelos de erosión de suelos se puede aplicar, también, para la modelización del lavado de sólidos en zonas urbanas. Para este tipo de aplicación se asume una única capa de sedimento no-cohesivo que representaría los sólidos sobre el pavimento. El pavimento sería la capa inferior, que en este caso es no erosionable, por lo que no es necesario incluirla en el cálculo.

Bajo estas condiciones, los términos de producción desde la capa inferior son igual a cero ( <math display="inline">{D}_{rdd}=</math><math>{D}_{fde}=0</math>). Si se utiliza la formulación de Harisine-Rose para el lavado por escorrentía obtenemos un modelo con 5 parámetros: <math display="inline">{\alpha }_{rd},\, {z}_{m},\, {M}_{s,cr},\, {\Omega }_{cr},\, F</math>, similar al que se utilizó en Naves et al. [51]. En dicha publicación se pueden, además, encontrar rangos de variación realistas de dichos parámetros en un entorno urbano controlado en laboratorio.

<span id='_Toc108335483'></span>

===5.2 Interfaz de usuario===

A continuación, se describen los diferentes menús y ventanas desde las que se pueden activar los procesos descritos en las secciones anteriores y definir los correspondientes parámetros.

<span id='_Toc108335484'></span>

====5.2.1 Formulaciones, opciones y parámetros generales del modelo====

Para definir las opciones y parámetros generales del caso a modelar, debe accederse a la ventana de ‘Datos del problema’, desde el menú de '''Datos >> Datos del problema >> Erosión de suelos'''. Desde dicha ventana se puede activar el módulo de erosión de suelos y seleccionar los diferentes modelos y procesos a considerar en el cálculo, así como los resultados que se desean ver en la interfaz de post-proceso. Las principales opciones a las que se accede desde esta ventana son las siguientes (<span id='cite-_Ref107218712'></span>[[#_Ref107218712|Fig. 33]]):

:* Activar/desactivar el módulo de erosión de suelos y seleccionar la asignación manual o automática de las zonas de erosión.
:* Estructura vertical del suelo o modelo de mezcla.
:* Propiedades físicas del suelo y parámetros de transporte (para cada clase de sedimento si en el cálculo se considera granulometría distribuida).
:* Activar/desactivar el transporte en suspensión, así como definir las formulaciones y parámetros correspondientes.
:* Activar/desactivar el transporte de fondo, así como definir las formulaciones y parámetros correspondientes.
:* Activar/desactivar diferentes correcciones a aplicar en las formulaciones utilizadas en el modelo (opcionales).
:* Seleccionar las variables que se guardarán en los ficheros de resultados para poder verlas posteriormente en el post-proceso.

<span id='_Ref107218712'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image33.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 33. Datos generales del módulo de erosión de suelos.'''</span>


<span id='_Toc108335485'></span>

====5.2.2 Asignación de zonas erosionables====

Una vez definidas las formulaciones y datos generales que se utilizarán en el cálculo, deben definirse y asignarse las zonas del dominio susceptibles de ser erosionadas (zonas erosionables). Esto se hace desde el menú '''Datos >> Erosión de suelos >> Zonas erosionables'''.

A través de la ventana que se muestra en la <span id='cite-_Ref107219070'></span>[[#_Ref107219070|Fig. 34]]a pueden crearse tantas zonas erosionables como se desee. Es recomendable darle un nombre a cada zona erosionable diferente al nombre que se genera por defecto. Para cada zona erosionable debe definirse la cantidad de material erosionable (kg/m<sup>2</sup>) disponible inicialmente para cada tipo de sedimento en la capa superior de suelo. Esta cantidad inicial de sedimento variará a lo largo del cálculo en función de los procesos de erosión y sedimentación que tengan lugar. Asimismo, se debe establecer si este material erosionable está situado sobre la cota definida por el MDT o por debajo de dicha cota. Por último, es necesario definir la concentración inicial de sólidos en suspensión en dicha zona (kg/m<sup>3</sup>). Lo más común en cálculos hidrológicos es que las zonas estén secas al comienzo del cálculo, en cuyo caso la concentración inicial de sólidos en suspensión será igual a cero.

Si un tipo de sedimento no está presente en dicha zona no es necesario incluirlo en la tabla.

'''''Asignación manual'''''

La asignación manual de zonas erosionables (ya deben estar creadas) se realiza mediante el menú '''Datos >> Erosión de suelos >> Asignación manual de zonas de erosión'''. Este menú da acceso a la ventana que se muestra en la <span id='cite-_Ref107219070'></span>[[#_Ref107219070|Fig. 34]]b, en la cual se puede elegir una zona erosionable creada anteriormente y asignarla a diferentes regiones del modelo. Al igual que otras propiedades de Iber, la asignación puede realizarse en la geometría o en la malla. En general es preferible realizar la asignación sobre la geometría del modelo y posteriormente generar la malla de cálculo.

<span id='_Ref107219070'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image34.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 34. (a) Definición de zonas erosionables. (b) Ventana de signación manual de zonas erosionables.'''</span>


'''''Asignación automática'''''

Para asignar las zonas erosionables de forma automática es necesario disponer de un archivo raster en formato ASCII Grid de ArcInfo (''*.asc'') que contenga la distribución espacial de las zonas en base a su identificador (ID). El identificador de cada zona puede conocerse cambiando el nombre de dicha zona (como se ha comentado anteriormente es recomendable asignarle un nombre a cada zona erosionable diferente al que se genera por defecto), apareciendo automáticamente a modo de prefijo del nombre asignado como puede verse en la <span id='cite-_Ref107219430'></span>[[#_Ref107219430|Fig. 35]]a. En el ejemplo mostrado en la <span id='cite-_Ref107219430'></span>[[#_Ref107219430|Fig. 35]]b se le ha dado el nombre ''Gravera'', y al darle a Aceptar el programa la guarda como ''ID2 – Gravera'' (ya que existía previamente otra zona con el identificador ID1). Si se crea una nueva zona tendrá el prefijo ID3.

<span id='_Ref107219430'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image35.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 35. (a) Ventana de cambio de nombre de una zona erosionable. (b) Ventana de asignación según el prefijo ID de la zona erosionable predefinida.'''</span>


Sabiendo cual es el ID de cada zona y disponiendo del archivo ráster con su distribución espacial, accedemos a la asignación automática mediante el menú '''Datos >> Erosión de suelos >> Asignación automática de zonas de erosión''', desde dónde se deberá seleccionar el archivo ráster. Como siempre en Iber, la asignación automática de propiedades se realiza directamente sobre la malla de cálculo, por lo que esta deberá estar creada con anterioridad. Asimismo, cada vez que se cree la malla de cálculo, deben realizarse todas las asignaciones automáticas de nuevo.

<span id='_Toc108335486'></span>

====5.2.3 Condiciones de contorno====

Una vez definidos los datos generales y asignado las zonas erosionables, deben asignarse las condiciones de contorno, tanto para el transporte en suspensión como para el transporte de fondo. Dichas condiciones deben asignarse únicamente en los contornos de entrada de caudal. En caso de no asignarse se considerará que no entra sedimento por dichos contornos.

Las condiciones de contorno se asignan desde el menú '''Datos >> Erosión de suelos >> Condiciones de contorno'''. Este menú da acceso a la ventana donde se puede elegir entre asignar las condiciones de contorno para el transporte en suspensión o de fondo.

'''''Condiciones de contorno para el transporte en suspensión'''''

Para asignar las condiciones de contorno de transporte en suspensión debe definirse la concentración de sólidos en suspensión (para cada clase de sedimento considerada) en los contornos de entrada de caudal. Para ello se define la concentración de sólidos totales en suspensión (kg/m<sup>3</sup>) como una serie temporal (<span id='cite-_Ref107220428'></span>[[#_Ref107220428|Fig. 36]]a). En el caso de que se trabaje con varias clases de sedimentos, dicha concentración corresponde a la suma de todas las clases, por lo que también será necesario definir la fracción de cada clase de sedimento en dicho contorno, en la columna ''Fracción de material''. La suma de los valores introducidos para las clases de sedimento en la columna ''Fracción de material'' debe ser igual a 1.

<span id='_Ref107220428'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image36.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 36. Ventana de asignación de condiciones de contorno para transporte de sedimentos por: (a) carga en suspensión; (b) carga de fondo.'''</span>


'''''Condiciones de contorno para el transporte de fondo'''''

Para el transporte de fondo se debe de asignar, en cada contorno de entrada de caudal, la evolución temporal del caudal sólido total que entra por dicho contorno (kg/s), y la fracción de dicho caudal que corresponde a cada tipo de sedimento (<span id='cite-_Ref107220428'></span>[[#_Ref107220428|Fig. 36]]b). Al igual que en el caso de transporte en suspensión, la suma de las fracciones de sedimento (columna ''Fracción de material'') debe ser igual a 1.

<span id='_Toc108335487'></span>
<span id='_Ref106264444'></span>

==6 Módulo de transporte de sedimento de fondo para mezclas==

<span id='_Toc108335488'></span>

===6.1 Fundamentos teóricos===

En este apartado se presenta el módulo de transporte de fondo para material no uniforme, es decir, cuando el sedimento transportado por un río no puede caracterizarse por un único diámetro. Es un hecho que un río que es alimentado con sedimento con distinto tamaño de grano, tiene la capacidad de clasificarlo. Esto es especialmente significativo en los ríos con lecho formado por una mezcla de gravas y arenas, pero también ocurre para otros tamaños, como pueden ser bolos, o gravas con distintos tamaños característicos. Esta clasificación acaba resultando en que, para muchos ríos, el fondo se encuentra formado por estratos con distintas granulometrías, con su interfaz sensiblemente horizontal, y una capa de sedimento grueso (coraza) en su parte superior.

Los aspectos relacionados con la variedad de diámetros de sedimento de un río tienen repercusiones importantes en distintas facetas de la ingeniería de ríos, como pueden ser las restauraciones fluviales, la preservación de hábitos, la propia evolución geomorfológica del río o en la interacción entre el material transportado por el fondo y las estructuras construidas en el mismo.

<span id='_Toc108335489'></span>

====6.1.1 Concepto de la capa activa====

En el transporte de sedimentos por arrastre de fondo para material distribuido, es decir, para mezclas de sedimento, es un fenómeno que se produce en un espesor determinado del fondo del río [52]. Es común explicar este fenómeno a través del concepto de la ''capa activa''. Aunque existen distintas interpretaciones de la capa activa, la más utilizada es que se puede considerar como capa activa el estrato superior del fondo de un río que se mueve y que es donde la clasificación del material realmente tiene lugar (<span id='cite-_Ref107221871'></span>[[#_Ref107221871|Fig. 37]]a). De esta manera, se suele llamar capa activa a la capa superior del fondo del río. Desde el punto de vista de la modelización numérica, su espesor se especifica mediante un parámetro, N, multiplicado por D<sub>50</sub> o D<sub>90</sub> del tamaño más grande utilizado en la simulación. La capa se suele considerar constante durante una simulación, aunque la composición de la capa activa en cada sección transversal es una función del tiempo e incorpora los efectos acumulativos de la erosión selectiva y la deposición [52,53].

<span id='_Ref107221871'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image37.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 37. Teoría de la capa activa: (a) esquema de la conceptualización de la estructura del fondo de un río bajo la idea de transporte por capa activa (fuente: Spasojevic y Holly [54]); (b) división de la estructura del fondo de un río en capa activa y estratos.'''</span>


<span id='_Toc108335490'></span>

====6.1.2 Discretización del sustrato====

Para poder calcular el transporte de fondo utilizando el concepto de capa activa el terreno erosionable, el sustrato se divide en una serie de estratos (<span id='cite-_Ref107221871'></span>[[#_Ref107221871|Fig. 37]]b). Cada capa, en cada instante y en cada punto, tendrá su granulometría. La definición de la granulometría se define en base a una serie de ''clases'' o diámetros característicos del sedimento. Es decir, en cada punto de una capa o estrato, habrá un porcentaje (en peso) de cada diámetro característico.

<span id='_Toc108335491'></span>

====6.1.3 Ecuaciones====

El transporte sólido total (caudal solido) en cada punto <math display="inline">{q}_{b}\,</math> es la suma del caudal sólido para cada una de las clases o diámetros <math display="inline">{q}_{ks}</math> . La proporción o fracción de caudal de cada clase se puede expresar como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\, {\mathit{\boldsymbol{q}}}_{b}=\sum _{k=1}^{n}{q}_{{k}_{s}}\, \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(109)
|}

donde <math display="inline">{q}_{{k}_{s}}</math> es la capacidad de transporte calculada con alguna fórmula de transporte sólido.

La ecuación de Exner, que representa la continuidad del caudal sólido en el fondo de un río, para el caso de mezclas, se escribe como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\left( 1-p\right) \frac{\partial {z}_{b}}{\partial t}+\sum _{{k}_{s}=1}^{n{^\circ}classes}\frac{\partial {q}_{{k}_{s},x}}{\partial x}+\frac{\partial {q}_{{k}_{s},y}}{\partial y}=0\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(110)
|}

La ecuación que representa el balance de masa total en un punto de un río, si se realiza por cada una de las clases o diámetros, es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\left( 1-p\right) \frac{\partial \left( {\beta }_{{k}_{s}}{h}_{a}\right) }{\partial t}+</math><math>\frac{\partial {q}_{{k}_{s},x}}{\partial x}+\frac{\partial {q}_{{k}_{s},y}}{\partial y}=</math><math>{\boldsymbol{F}}_{{k}_{s}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(111)
|}

En las expresiones anteriores <math display="inline">{q}_{{k}_{s}}</math>, es el caudal sólido unitario de la clase <math display="inline">{k}_{s}</math>'', '' <math display="inline">{\beta }_{{k}_{s}}</math> la fracción de dicha clase <math display="inline">{k}_{s}</math> en la capa active y <math display="inline">{F}_{{k}_{s}}</math> la fuente de sedimentos de la clase <math display="inline">{k}_{s}</math> hacia la capa activa procedente del substrato inferior y <math display="inline">p</math> representa la porosidad.

El intercambio de material entre la capa activa (móvil) y la capa inferior (estacionaria) se puede determinar según:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{F}_{{k}_{s}}=-\left( 1-p\right) \frac{\partial \left[ {\beta }^{'}\left( {z}_{b}-{h}_{a}\right) \right] }{\partial t}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(112)
|}

siendo <math display="inline">{\beta }^{'}</math>:

:* <math display="inline">{\beta }_{s{k}_{s}}</math>, si el flujo es ascendente, del substrato hacia la capa activa
:* <math display="inline">{\beta }_{{k}_{s}}</math>, si el flujo es descendente, desde la capa activa al sustrato

<math display="inline">{\beta }_{s{k}_{s}}</math>, es la fracción de la clase <math display="inline">{k}_{s}</math> en el estrato activo (estrato en el que se encuentra la capa activa en cada instante y punto). Debe cumplirse que:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\sum _{{k}_{s}=1}^{n{^\circ}\, classes}{\beta }_{{k}_{s}}=1\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(113)
|}

El caudal sólido se puede calcular con cualquiera de las fórmulas de transporte que incluye Iber (Meyer-Peter y Müller, Engelund y Fredsoe, Engelund y Hansen, van Rijn, Recking, y una formulación ad-hoc donde el usuario puede modificar los coeficientes). En la versión 3 se ha añadido la de Ashida y Michiue [55], por ser una de las comúnmente utilizadas para mezclas. Con ella, el caudal sólido adimensional para una clase <math display="inline">i</math> se define como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{i}^{\ast }=17\left( {\tau }_{{s}_{i}}^{\ast }-{\tau }_{{sc}_{i}}^{\ast }\right) \left( \sqrt{{\tau }_{{s}_{i}}^{\ast }}-\right. </math><math>\left. \sqrt{{\tau }_{{sc}_{i}}^{\ast }}\right) \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(114)
|}

A partir de allí, el caudal sólido de la clase es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{q}_{i}={q}_{i}^{\ast }\left( {F}_{i}\sqrt{\gamma g{D}_{i}}{D}_{i}\right) \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(115)
|}

Siendo, por tanto, el caudal sólido total:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>q=\sum _{}^{}{q}_{i}\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(116)
|}

<span id='_Toc108335492'></span>

====6.1.4 Escondimiento (hiding)====

En el cálculo del caudal sólido para mezclas de sedimentos es importante considerar el hecho de que una partícula de un diámetro característico no se mueve de la misma manera cuando forma parte de una mezcla, como lo haría si el resto de partículas tuvieran el mismo tamaño. Este efecto se conoce como escondimiento, ampliamente estudiado por Ashida y Michiue [55], y más adelante por otros muchos autores.

Por consistencia, en Iber se ha implementado la ecuación de escondimiento de los primeros, pero también la interpretación de la misma según otros autores, incluso la posibilidad de modificar los parámetros a criterio del usuario. Según Ashida y Michiue, la relación entre la tensión tangencial crítica adimensional para una clase <math display="inline">i</math>, <math display="inline">{\tau }_{{sc}_{i}}^{\ast }</math>, se relaciona con la tensión tangencial crítica adimensional que se tendría con un diámetro característico <math display="inline">{D}_{g}^{}</math>, <math display="inline">{\tau }_{{sc}_{g}}^{\ast }.\, {D}_{g}^{}\,</math>  corresponde a la media geométrica de la distribución granulométrica del estrato según Egiazaroff [56]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{{\tau }_{{sc}_{i}}^{\ast }}{{\tau }_{{sc}_{g}}^{\ast }}=\frac{{D}_{i}^{}}{{D}_{g}^{}}\, {\left[ \frac{\mathrm{ln}\,(19)}{\mathrm{ln}\,\left( 19\frac{{D}_{i}^{}}{{D}_{g}^{}}\right) }\right] }^{2}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(117)
|}

Según Ashida y Michiue [55], la expresión sería:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}Si\, \frac{{D}_{i}^{}}{{D}_{g}^{}}<0,4&\frac{{\tau }_{{sc}_{i}}^{\ast }}{{\tau }_{{sc}_{g}}^{\ast }}={0,843\left( \frac{{D}_{i}^{}}{{D}_{g}^{}}\, \right) }^{-1}\\Si\, \frac{{D}_{i}^{}}{{D}_{g}^{}}\geq 0,4&\frac{{\tau }_{{sc}_{i}}^{\ast }}{{\tau }_{{sc}_{g}}^{\ast }}=\frac{{D}_{i}^{}}{{D}_{g}^{}}\, {\left[ \frac{\mathrm{ln}\,(19)}{\mathrm{ln}\,\left( 19\frac{{D}_{i}^{}}{{D}_{g}^{}}\right) }\right] }^{2}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(118)
|}

Según Parker, Klingeman y McLean [57]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{{\tau }_{{sc}_{i}}^{\ast }}{{\tau }_{{sc}_{g}}^{\ast }}={\left( \frac{{D}_{i}^{}}{{D}_{g}^{}}\, \right) }^{-0,982}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(119)
|}

Según Wilcock y Crowe [<span id='cite-_Bib058'></span>[[#_Bib058|58]]]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{{c}_{i}}^{\ast }={\tau }_{{c}_{g}}^{\ast }{\left( \frac{{D}_{i}}{{D}_{g}}~\right) }^{-~\left( \frac{0,67}{1-{e}^{1,5-\frac{{D}_{i}}{{D}_{g}}}}\right) }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(120)
|}

Asimismo, se ha implementado la posibilidad de considerar una fórmula ad-hoc, con una expresión de forma similar a la anterior, pero con coeficientes ajustables por el usuario ( <math display="inline">\alpha</math>  y <math display="inline">\beta</math> ):

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{{\tau }_{{sc}_{i}}^{\ast }}{{\tau }_{{sc}_{g}}^{\ast }}=\alpha {\left( \frac{{D}_{i}^{}}{{D}_{g}^{}}\, \right) }^{\beta }</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(121)
|}

<span id='_Toc54965864'></span><span id='_Toc90987698'></span><span id='_Toc108335493'></span>

====6.1.5 Implementación de las ecuaciones en el esquema numérico====

La ecuación de Exner presentada anteriormente no es nada más que la ecuación de conservación de la masa para el sedimento arrastrado por el fondo del río. Mediante el método de los volúmenes finitos, presentado en el manual de referencia hidráulico de Iber, resulta que la ecuación para actualizar cota del fondo del río <math display="inline">{z}_{i}^{n}</math> en cada instante se puede escribir como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{z}_{i}^{n+1}={z}_{i}^{n}-\frac{\Delta t}{\left( 1-p\right) {V}_{i}}\sum _{{k}_{s}=1}^{classes}\sum _{w=1}^{~lados}\left( {\overline{q}}_{{k}_{s},i,w}{\overline{n}}_{i,w}\right) {l}_{i,w}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(122)
|}

donde <math display="inline">{\overline{q}}_{{k}_{s},i,w}</math>, es el caudal sólido unitario para la clase <math display="inline">{k}_{s}</math>, en el lado (o pared) ''w'' del elemento (o volumen finito) <math display="inline">i</math>. Para la capa activa, la actualización del porcentaje de cada clase, que se puede considerar como una ecuación de Exner local. Se expresará entonces como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\left( {h}_{a}{\beta }_{{k}_{s}i}^{n+1}\right) =\left( {h}_{a}{\beta }_{{k}_{s}i}^{n}\right) -</math><math>\frac{\Delta t}{\left( 1-p\right) {V}_{i}}\sum _{w=1}^{lados}\left( {\overline{q}}_{{k}_{s},i,w}{\overline{n}}_{i,w}\right) -{F}_{{k}_{s}i}^{n}\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(123)
|}

La fuente de fondo, o intercambio entre el estrato activo y la capa activa, responde a:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{F}_{{k}_{s}i}^{n+1}=-\left( 1-p\right) \left[ \frac{\left( {\beta }_{i}^{n+1}\left( {z}_{b}^{n+1}-{h}_{a}^{n+1}\right) \right) -\left( {\beta }_{i}^{n}\left( {z}_{b}^{n}-{h}_{a}^{n}\right) \right) }{\Delta t}\right] \,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(124)
|}

<math display="inline">\beta</math> , desde el punto de vista numérico, toma la forma:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}{\beta }_{{sk}_{s\, }\, }\, \, si\quad \, \, \frac{\left( {z}_{b}^{n+1}-{h}_{a}^{n+1}\right) -\left( {z}_{b}^{n}-{h}_{a}^{n}\right) }{\Delta t}>0\, \, \\{\beta }_{{k}_{s\, }\, }\, \, si\quad \, \, \frac{\left( {z}_{b}^{n+1}-{h}_{a}^{n+1}\right) -\left( {z}_{b}^{n}-{h}_{a}^{n}\right) }{\Delta t}<0\, \end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(125)
|}

Para la actualización de la fracción de cada clase en el estrato activo, las expresiones resultantes son:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\begin{matrix}Si\, {\, \left( {z}_{b}-{h}_{a}\right) }^{n+1}<{\left( {z}_{b}-{h}_{a}\right) }^{n}&\rightarrow &{\beta }_{{sk}_{s}}^{n+1}={\beta }_{{sk}_{s}}^{n}\\Si\, \, {\left( {z}_{b}-{h}_{a}\right) }^{n+1}>{\left( {z}_{b}-{h}_{a}\right) }^{n}&\rightarrow &{\beta }_{{sk}_{s}}^{n+1}={\beta }_{{sk}_{s}}^{n}{\left( {z}_{b}-{h}_{a}-{z}_{{s}_{1}}\right) }^{n}\end{matrix}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(126)
|}

siendo,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\beta }_{{sk}_{s}}^{n+1}=\frac{{F}_{{k}_{s}}^{n+1}\Delta t+{\alpha }_{{sk}_{s}}^{n}}{{\left( {z}_{b}-{h}_{a}-{z}_{{s}_{1}}\right) }^{n+1}}\,</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(127)
|}
<span id='_Toc108335494'></span>

===6.2 Interfaz de usuario===

A continuación, se describen los diferentes menús y ventanas desde las que se pueden activar los procesos descritos en las secciones anteriores, así como definir los correspondientes parámetros. Este nuevo módulo se ha creado en base a tablas dinámicas y, por tanto, requiere de un proceso secuencial para su correcto funcionamiento. Para evitar problemas para con la implementación de las condiciones, se han implementado mensajes de ayuda para que el usuario sepa cómo debe proceder.

<span id='_Toc54965866'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image38.png|550px]] </div>

<span id='_Ref107224388'></span><span style="text-align: center; font-size: 75%;">'''Fig. 38. Datos generales del módulo de transporte de sedimentos (mezclas).'''</span>


<span id='_Toc90987700'></span><span id='_Toc108335495'></span>

====6.2.1 Formulaciones, opciones y parámetros generales del modelo====

Para definir las opciones y parámetros generales del caso a modelar, debe accederse a la ventana de ‘Datos del problema’, desde el menú de '''Datos >> Datos del problema >> Sedimentos'''. Desde dicha ventana se puede activar el módulo de transporte de sedimentos, ya sea por suspensión como por fondo. En el caso de mezclas de sedimentos, esta opción requiere la activación del transporte de fondo y seleccionar en ''Tipo de sedimento'' la opción '''Mezclas'''. A continuación, se requiere la introducción de diferentes parámetros, que también son comunes para granulometría uniforme. Las principales opciones a las que se accede desde esta ventana son las siguientes (<span id='cite-_Ref107224388'></span>[[#_Ref107224388|Fig. 38]]):

:* Activar/desactivar el módulo de transporte de sedimentos.
:* Selección del tipo de sedimento.
:* Propiedades del sedimento, como son la porosidad, el ángulo de fricción y la densidad relativa del sedimento.
:* Asimismo, se ofrece la posibilidad de acondicionar el terreno (pendiente) según las propiedades del sedimento a través de la activación del Modelo de Avalancha

<span id='_Toc54965867'></span><span id='_Toc90987701'></span><span id='_Toc108335496'></span>

====6.2.2 Definición del número de estratos, clases y espesor de la capa activa====

La definición del número de estratos y clases del sedimento, así como el espesor de la capa activa, es primordial para la implementación de las propiedades del sedimento, ya que de estos valores depende la generación del resto de ventanas.

A esta ventana se accede desde el menú '''Datos >> Transporte de sedimentos >> Transporte de fondo (mezclas) >> Estratos&Clases'''. En la ventana que se muestra en la <span id='cite-_Ref107224079'></span>[[#_Ref107224079|Fig. 39]]a deben definirse el número de estratos y clases del sedimento, así como el espesor de la capa activa (en metros). Estos parámetros son fijos para todo el modelo.

<span id='_Ref107224079'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image39.png|400px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 39. (a) Definición del número de estratos, clases y espesor de la capa activa. (b) Definición del diámetro de cada clase.'''</span>


<span id='_Toc108335497'></span>

====6.2.3 Definición de las clases====

La definición del número y diámetro de cada clase se realiza mediante el menú '''Datos >> Transporte de sedimentos >> Transporte de fondo (mezclas) >> Definición de clases. '''Este menú da acceso a la ventana que se muestra en la <span id='cite-_Ref107224079'></span>[[#_Ref107224079|Fig. 39]]b, en la cual se debe introducir la relación de cada clase y el diámetro de partícula. El valor del diámetro de cada clase es constante, pero no así el porcentaje de cada uno de ellos en cada capa, ya que variará a lo largo del proceso de cálculo.

<span id='_Toc108335498'></span>

====6.2.4 Propiedades de las mezclas====

A continuación, se deben definir las propiedades del sedimento, es decir, introducir el espesor de cada capa y la relación de porcentajes que definen la granulometría en base al número de clases que se ha predefinido. Esto se realiza desde el menú '''Datos >> Transporte de sedimentos >> Transporte de fondo (mezclas) >> Propiedades de las mezclas''' (<span id='cite-_Ref107224997'></span>[[#_Ref107224997|Fig. 40]]a). En esta ventana se deben definir las propiedades de cada capa, indicando el identificador de la capa, el espesor de la capa (en metros) y el porcentaje (%) de cada clase, que está asociado al diámetro de la partícula previamente definido.

<span id='_Ref107224997'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image40.png|100%]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 40. (a) Definición de las propiedades de cada estructura del suelo. (b) Asignación de la estructura del suelo predefinida.'''</span>


En esta ventana también se debe especificar la distancia entre la capa superior del estrato 1 y el fondo, y que determina la posición de los estratos. Tiene que ser un número mayor que cero. Si es cero, es que el fondo coincide con este estrato. Si es mayor que cero, quiere decir que en el momento inicial el fondo está erosionado, o que puede haber sedimentación, por lo que pueden haber estratos por encima de la cota de fondo inicial.

Se pueden crear tantas “estructuras de suelo” de sedimentos como se desee, siendo el número de estratos y de clases fijo para todo el modelo.

<span id='_Toc108335499'></span>

====6.2.5 Asignación de las mezclas====

Una vez definidas las “mezclas”, es posible asignarlas mediante el menú '''Datos >> Transporte de sedimentos >> Transporte de fondo (mezclas) >> Asignación de mezclas'''. Esto se puede hacer tanto en la geometría como en la malla de cálculo (<span id='cite-_Ref107224997'></span>[[#_Ref107224997|Fig. 40]]b). Todas aquellas “mezclas” que se hayan creado en la ventana anterior, se mostrarán en una lista para poder asignarlas al modelo.

<span id='_Toc108335500'></span>

====6.2.6 Condiciones de contorno====

Una vez definidos los datos generales y definidas/asignadas las propiedades del sedimento, deben asignarse las condiciones de contorno. Dichas condiciones deben asignarse únicamente en los contornos de entrada de caudal. En caso de no asignarse, se considerará que no entra sedimento por dichos contornos.

Las condiciones de contorno se asignan desde el menú '''Datos >> Transporte de sedimentos >> Transporte de fondo (mezclas) >> Condiciones de contorno (mezclas)'''. Este menú da acceso a la ventana que se muestra en la <span id='cite-_Ref107227012'></span>[[#_Ref107227012|Fig. 41]], desde donde se puede elegir entre asignar las condiciones de contorno en función del tipo de entrada: no-sedimentos (o agua clara), por capacidad de arrastre y como caudal mediante un sedimentograma.

'''''No –Agua clara–'''''

Esta opción introducirá caudal de agua, pero no de sedimentos. Es decir, no habrá aporte de sedimentos por los contornos de entrada de caudal.

'''''Capacidad de arrastre'''''

Con esta opción, el modelo considera la entrada de sedimentos en los contornos de entrada que se indique, siendo el tipo de sedimento (mezcla) idéntico al que hay en el elemento interior de dicho contorno. La cantidad de sedimento entrante dependerá, por tanto, del tipo de sedimento impuesto en la zona inmediatamente adyacente a la condición de contorno, así como a las propiedades del mismo.

'''''Sedimentograma'''''

Finalmente, se proporciona la posibilidad de imponer una condición de contorno ad hoc. En este caso se debe definir el caudal de entrada de material sólido, es decir, un sedimentograma. El caudal sólido de entrada total será el especificado. La granulometría será la asignada en el elemento junto a la entrada.

<span id='_Ref107227012'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image41.png|350px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 41. Ventana de definición e implementación de las condiciones de contorno para mezclas.'''</span>


<span id='_Toc108335501'></span>

====6.2.7 Resultados====

Para el caso de sedimento con granulometría no-uniforme, los resultados específicos aparecen dentro del análisis con el nombre de '''Sedimentos'''. Con las herramientas habituales de consulta de resultados se pueden analizar los resultados globales o por clases (<span id='cite-_Ref107225581'></span>[[#_Ref107225581|Fig. 42]]). La particularidad de estos resultados, accesibles desde el análisis ''Sedimentos'', es que se han agrupado por estratos (incluyendo un análisis propio para la capa activa). En cada estrato se puede consultar la proporción de cada clase de sedimento, hallándose cada uno de ellos a su correspondiente profundidad.

<span id='_Ref107225581'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image42-c.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 42. Organización de los resultados para mezclas.'''</span>


Asimismo, se ha creado una opción nueva, accesible a través del menú '''Herramientas Iber >> Transporte de Sedimentos >> Granulometría''', que permite seleccionar unos elementos y ver la granulometría en ellos, tanto para la capa activa como para cada uno de los estratos existentes (<span id='cite-_Ref107225792'></span>[[#_Ref107225792|Fig. 43]]).

<span id='_Ref107225792'></span>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanz-Ramos_968687760-image43.png|600px]] </div>

<span style="text-align: center; font-size: 75%;">'''Fig. 43. Distribución granulométrica de una capa del sustrato en un modelo simulado con granulometría no uniforme.'''</span>


<span id='_Toc108335502'></span>

==Referencias==

<span id='_Bib001'></span>
[1] E. Bladé Castellet, L. Cea, G. Corestein, Numerical modelling of river inundations, Ing. Del Agua. 18 (2014) 68. doi: [https://doi.org/10.4995/ia.2014.3144 10.4995/ia.2014.3144].

<span id='_Bib002'></span>
[2] E. Bladé, L. Cea, G. Corestein, E. Escolano, J. Puertas, E. Vázquez-Cendón, J. Dolz, A. Coll, Iber: herramienta de simulación numérica del flujo en ríos, Rev. Int. Métodos Numéricos Para Cálculo y Diseño En Ing. 30 (2014) 1–10. doi: [https://doi.org/10.1016/j.rimni.2012.07.004 10.1016/j.rimni.2012.07.004].

<span id='_Bib003'></span>
[3] M. Sanz-Ramos, G. Olivares Cerpa, E. Bladé, Metodología para el análisis de rotura de presas con aterramiento mediante simulación con fondo móvil, Ribagua. 6 (2019) 138–147. doi: [https://doi.org/10.1080/23863781.2019.1705198 10.1080/23863781.2019.1705198].

<span id='_Bib004'></span>
[4] J. Anta Álvarez, M. Bermúdez, L. Cea, J. Suárez, P. Ures, J. Puertas, Modelización de los impactos por DSU en el río Miño (Lugo), Ing. Del Agua. 19 (2015) 105. doi: [https://doi.org/10.4995/ia.2015.3648 10.4995/ia.2015.3648].

<span id='_Bib005'></span>
[5] L. Cea, M. Bermúdez, J. Puertas, E. Bladé, G. Corestein, E. Escolano, A. Conde, B. Bockelmann-Evans, R. Ahmadian, IberWQ: new simulation tool for 2D water quality modelling in rivers and shallow estuaries, J. Hydroinformatics. 18 (2016) 816–830. doi: [https://doi.org/10.2166/hydro.2016.235 10.2166/hydro.2016.235].

<span id='_Bib006'></span>
[6] M. Sanz-Ramos, E. Bladé Castellet, A. Palau Ibars, D. Vericat Querol, A. Ramos-Fuertes, IberHABITAT: evaluación de la Idoneidad del Hábitat Físico y del Hábitat Potencial Útil para peces. Aplicación en el río Eume, Ribagua. 6 (2019) 158–167. doi: [https://doi.org/10.1080/23863781.2019.1664273 10.1080/23863781.2019.1664273].

<span id='_Bib007'></span>
[7] M. Sanz-Ramos, E. Bladé, L. Cea, Iberaula - Dissemination, (2022). [https://iberaula.es/56/iber-community/dissemination https://iberaula.es/56/iber-community/dissemination] (accessed July 5, 2022).

<span id='_Bib008'></span>
[8] A. Coll, R. Ribó, M. Pasenau, E. Escolano, J.S. Perez, A. Melendo, A. Monros, J. Gárate, GiD v.14 User Manual, (2018). [https://www.gidsimulation.com/ https://www.gidsimulation.com/]

<span id='_Bib009'></span>
[9] A. Coll, M. Pasenau, E. Escolano, J.S. Perez, A. Melendo, A. Monros, J. Gárate, [http://www.gidhome.com www.gidhome.com], (2018). [http://www.gidhome.com. www.gidhome.com.]

<span id='_Bib010'></span>
[10] A. Coll, R. Ribó, M. Pasenau, E. Escolano, J.S. Perez, A. Melendo, A. Monros, J. Gárate, GiD v.14 Reference Manual, (2018). [https://www.gidsimulation.com/ https://www.gidsimulation.com/]

<span id='_Bib011'></span>
[11] J.L. Aragón-Hernández, E. Bladé, Modelación numérica de flujo mixto en conductos cerrados con esquemas en volúmenes finitos, Tecnol. y Ciencias Del Agua. 08 (2017) 127–142. doi: [https://doi.org/10.24850/j-tyca-2017-03-08 10.24850/j-tyca-2017-03-08].

<span id='_Bib012'></span>
[12] E. Bladé, M. Sanz-Ramos, J. Dolz, J.M. Expósito-Pérez, M. Sánchez-Juny, Modelling flood propagation in the service galleries of a nuclear power plant, Nucl. Eng. Des. 352 (2019) 110180. doi: [https://doi.org/10.1016/j.nucengdes.2019.110180 10.1016/j.nucengdes.2019.110180].

<span id='_Bib013'></span>
[13] L. Cea, A. López‐Núñez, Extension of the two‐component pressure approach for modeling mixed free‐surface‐pressurized flows with the two‐dimensional shallow water equations, Int. J. Numer. Methods Fluids. 93 (2021) 628–652. doi: [https://doi.org/10.1002/fld.4902 10.1002/fld.4902].

<span id='_Bib014'></span>
[14] BOE-A-2008-755, Real Decreto 9/2008, de 11 de enero, por el que se modifica el Reglamento del Dominio Público Hidráulico, aprobado por el Real Decreto 849/1986, de 11 de abril, Boletín Of. Del Estado Núm. 14, 16 Enero 2008, Páginas 3141 a 3149. Minist. La Pres. (2008) 9. [https://www.boe.es/buscar/doc.php?id=BOE-A-2008-755 https://www.boe.es/buscar/doc.php?id=BOE-A-2008-755].

<span id='_Bib015'></span>
[15] ACA, Recomanacions tècniques per als estudis d’inundabilitat d’àmbit local. Guia Tècnica, Agència Catalana de l’Aigua. Generalitat de Catalunya. Març 2003, 2003. [http://www.gencat.net/aca. http://www.gencat.net/aca.]

<span id='_Bib016'></span>
[16] E. Bladé, M. Gómez-Valentín, J. Dolz, J.L.L. Aragón-Hernández, G. Corestein, M. Sánchez-Juny, Integration of 1D and 2D finite volume schemes for computations of water flow in natural channels, Adv. Water Resour. 42 (2012) 17–29. doi: [https://doi.org/10.1016/j.advwatres.2012.03.021 10.1016/j.advwatres.2012.03.021].

<span id='_Bib017'></span>
[17] E.F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, Springer: Berlin/Heidelberg, Germany, 2009. doi: [https://doi.org/10.1007/b79761 10.1007/b79761].

<span id='_Bib018'></span>
[18] M.H. Chaudhry, Open channel flow, Prentice Hall, New Jersey, USA, 1993.

<span id='_Bib019'></span>
[19] P.L. Roe, A basis for the upwind differencing of the two-dimensional unsteady Euler equations, Numer. Methods Fluid Dyn. 2. (1986) 55–80.

<span id='_Bib020'></span>
[20] E. Bladé, M. Gómez-Valentín, M. Sánchez-Juny, J. Dolz, Preserving steady-state in one-dimensional finite-volume computations of river flow, J. Hydraul. Eng. 134 (2008) 1343–1347. doi:  [https://doi.org/10.1061/(ASCE)0733-9429(2008)134:9(1343) 10.1061/(ASCE)0733-9429(2008)134:9(1343)].

<span id='_Bib021'></span>
[21] A. Preissmann, Propagation des intumescences dans les canaux et rivières (pp., in: 1st Congrèss Assoc. Fr. Calc. Grenoble, AFC, Paris, Fr. Sept., 1961: pp. 433–442.

<span id='_Bib022'></span>
[22] A.G. Barnett, Application of Godunov-type schemes to transient mixed flows, J. Hydraul. Res. 48 (2010) 686–688. doi: [https://doi.org/10.1080/00221686.2010.512793 10.1080/00221686.2010.512793].

<span id='_Bib023'></span>
[23] F. Kerger, P. Archambeau, S. Erpicum, B.J. Dewals, M. Pirotton, A fast universal solver for 1D continuous and discontinuous steady flows in rivers and pipes, Int. J. Numer. Methods Fluids. 66 (2011) 38–48. doi: [https://doi.org/10.1002/fld 10.1002/fld].

<span id='_Bib024'></span>
[24] J. Dolz, M. Gómez, Problemática del drenaje de aguas pluviales en zonas urbanas y del estudio hidráulico de las redes de colectores, Ing. Del Agua. 1 (1994) 55–66. doi: [https://doi.org/10.4995/ia.1994.2631  10.4995/ia.1994.2631].

<span id='_Bib025'></span>
[25] J. Despotovic, J. Plavsic, N. Stefanovic, D. Pavlovic, Inefficiency of storm water inlets as a source of urban floods, Water Sci. Technol. 51 (2005) 139–145.

<span id='_Bib026'></span>
[26] M. Gómez, B. Russo, Methodology to estimate hydraulic efficiency of drain inlets, Water Manag. 164 (2011) 81–90. [https://doi.org/10.1680/wama.900070  10.1680/wama.900070].

<span id='_Bib027'></span>
[27] J.Á. Aranda, C. Beneyto, M. Sánchez-Juny, E. Bladé, Efficient Design of Road Drainage Systems, Water. 13 (2021) 1661. doi: [https://doi.org/10.3390/w13121661  10.3390/w13121661].

<span id='_Bib028'></span>
[28] M. Gómez, B. Russo, Hydraulic Efficiency of Continuous Transverse Grates for Paved Areas, J. Irrig. Drain. Eng. 135 (2009) 225–230. doi: [https://doi.org/10.1061/(ASCE)0733-9437(2009)135:2(225) 10.1061/(ASCE)0733-9437(2009)135:2(225)].

<span id='_Bib029'></span>
[29] L. Cea, E. Bladé, A simple and efficient unstructured finite volume scheme for solving the shallow water equations in overland flow applications, Water Resour. Res. 51 (2015) 5464–5486. doi: [https://doi.org/10.1002/2014WR016547 10.1002/2014WR016547].

<span id='_Bib030'></span>
[30] V.T. Chow, D.R. Maidment, L.W. Mays, Applied Hydrology, MCGRAW-HIL, New York, USA, 1988.

<span id='_Bib031'></span>
[31] M.K. Kenneth, PART 630 HYDROLOGY, in: Natl. Eng. Handb., USDA, Natural resources Conservation Service, 2010: pp. 449–456.

<span id='_Bib032'></span>
[32] USDA-NRCS, Part 630 Hydrology - Chapter 10, in: Natl. Eng. Handb., US Department of Agriculture, Soil Conservation Service: Washington, DC, USA, 2004: p. 79.

<span id='_Bib033'></span>
[33] NCRS, Hydrology, in: Natl. Eng. Handb., US Department of Agriculture, Washington DC, USA, 1972.

<span id='_Bib034'></span>
[34] N. Kannan, C. Santhi, J.R. Williams, J.G. Arnold, Development of a continuous soil moisture accounting procedure for curve number methodology and its behaviour with different evapotranspiration methods, Hydrol. Process. 22 (2008) 2114–2121. doi: [https://doi.org/10.1002/hyp.6811 10.1002/hyp.6811].

<span id='_Bib035'></span>
[35] J. Doorenbos, W.O. Pruitt, Crop water requirements, FAO Irrig Drain. Paper 24. FAO of the United Nations, Rome, 1977.

<span id='_Bib036'></span>
[36] J.S. Famiglietti, E.F. Wood, Multiscale modeling of spatially variable water and energy balance processes, Water Resour. Res. 30 (1994) 3061–3078. doi: [https://doi.org/10.1029/94WR01498 10.1029/94WR01498].

<span id='_Bib037'></span>
[37] L.C. van Rijn, Sediment Transport, Part I: Bed Load Transport, J. Hydraul. Eng. 110 (1984) 1431–1456. doi: [https://doi.org/10.1061/(ASCE)0733-9429(1984)110:10(1431) 10.1061/(ASCE)0733-9429(1984)110:10(1431)].

<span id='_Bib038'></span>
[38] L.C. van Rijn, Mathematical modelling of morphological processes in the case of suspended sediment transport, (1987) 260.

<span id='_Bib039'></span>
[39] J.D. Smith, Modeling of sediment transport on continental shelves, United States, 1975. [https://www.osti.gov/biblio/7156268. https://www.osti.gov/biblio/7156268.]

<span id='_Bib040'></span>
[40] M. Garcia, G. Parker, Entrainment of Bed Sediment into Suspension, J. Hydraul. Eng. 117 (1991) 414–435. doi: [https://doi.org/10.1061/(ASCE)0733-9429(1991)117:4(414) 10.1061/(ASCE)0733-9429(1991)117:4(414)].

<span id='_Bib041'></span>
[41] J.A. Zyserman, J. Fredsøe, Data Analysis of Bed Concentration of Suspended Sediment, J. Hydraul. Eng. 120 (1994) 1021–1042. doi: [https://doi.org/10.1061/(ASCE)0733-9429(1994)120:9(1021) 10.1061/(ASCE)0733-9429(1994)120:9(1021)].

<span id='_Bib042'></span>
[42] G. Parker, G. Seminara, L. Solari, Bed load at low Shields stress on arbitrarily sloping beds: Alternative entrainment formulation, Water Resour. Res. 39 (2003). doi: [https://doi.org/10.1029/2001WR001253 10.1029/2001WR001253].

<span id='_Bib043'></span>
[43] E. Meyer-Peter, MüllerR., Formulas for Bed-Load transport, in: IAHSR 2nd Meet. Stock. Append. 2, IAHR, 1948.

<span id='_Bib044'></span>
[44] M. Wong, G. Parker, Reanalysis and Correction of Bed-Load Relation of Meyer-Peter and M?ller Using Their Own Database, J. Hydraul. Eng. 132 (2006) 1159–1168. doi: [https://doi.org/10.1061/(ASCE)0733-9429(2006)132:11(1159) 10.1061/(ASCE)0733-9429(2006)132:11(1159)].

<span id='_Bib045'></span>
[45] M.H. García, Sediment Transport and Morphodynamics, in: Sediment. Eng., American Society of Civil Engineers, Reston, VA, 2008: pp. 21–163. doi: [https://doi.org/10.1061/9780784408148.ch02 10.1061/9780784408148.ch02].

<span id='_Bib046'></span>
[46] G.J.C.M. Hoffmans, H.J. Verheij, Scour Manual, Routledge, 2017. doi: [https://doi.org/10.1201/9780203740132 10.1201/9780203740132].

<span id='_Bib047'></span>
[47] L. Cea, C. Legout, T. Grangeon, G. Nord, Impact of model simplifications on soil erosion predictions: application of the GLUE methodology to a distributed event-based model at the hillslope scale, Hydrol. Process. 30 (2016) 1096–1113. doi: [https://doi.org/10.1002/hyp.10697 10.1002/hyp.10697].

<span id='_Bib048'></span>
[48] J. Naves, J. Anta, J. Suárez, J. Puertas, Hydraulic, wash-off and sediment transport experiments in a full-scale urban drainage physical model, Sci. Data. 7 (2020) 44. doi: [https://doi.org/10.1038/s41597-020-0384-z 10.1038/s41597-020-0384-z].

<span id='_Bib049'></span>
[49] P.B. Hairsine, C.W. Rose, Modeling water erosion due to overland flow using physical principles: 2. Rill flow, Water Resour. Res. 28 (1992) 245–250. doi: [https://doi.org/10.1029/91WR02381 10.1029/91WR02381].

<span id='_Bib050'></span>
[50] P.B. Hairsine, C.W. Rose, Modeling water erosion due to overland flow using physical principles: 1. Sheet flow, Water Resour. Res. 28 (1992) 237–243. doi: [https://doi.org/10.1029/91WR02380 10.1029/91WR02380].

<span id='_Bib051'></span>
[51] J. Naves, J. Rieckermann, L. Cea, J. Puertas, J. Anta, Global and local sensitivity analysis to improve the understanding of physically-based urban wash-off models from high-resolution laboratory experiments, Sci. Total Environ. 709 (2020) 136152. doi: [https://doi.org/10.1016/j.scitotenv.2019.136152 10.1016/j.scitotenv.2019.136152].

<span id='_Bib052'></span>
[52] M. Hirano, River-bed degradation with armoring, Proc. Japan Soc. Civ. Eng. 1971 (1971) 55–65. doi: [https://doi.org/10.2208/jscej1969.1971.195_55 10.2208/jscej1969.1971.195_55].

<span id='_Bib053'></span>
[53] M. Church, J.K. Haschenburger, What is the “active layer”?, Water Resour. Res. 53 (2017) 5–10. doi: [https://doi.org/10.1002/2016WR019675 10.1002/2016WR019675].

<span id='_Bib054'></span>
[54] M. Spasojevic, F.M. Holly, Two- and Three-Dimensional Numerical Simulation of Mobile-Bed Hydrodynamics and Sedimentation, in: Sediment. Eng., American Society of Civil Engineers, Reston, VA, 2008: pp. 683–761. doi: [https://doi.org/10.1061/9780784408148.ch15 10.1061/9780784408148.ch15].

<span id='_Bib055'></span>
[55] K. Ashida, M. Michiue, Hydraulic Resistance of Flow in an Alluvia Bed and Bed Load Transport Rate, in: Proc. Japan Soc. Civ. Eng., Japan Society of Civil Engineers. No. 206 (in Japanese), 1972: pp. 59–69. [https://scirp.org/reference/referencespapers.aspx?referenceid=1312271. https://scirp.org/reference/referencespapers.aspx?referenceid=1312271.]

<span id='_Bib056'></span>
[56] I. V. Egiazaroff, Calculation of Nonuniform Sediment Concentrations, J. Hydraul. Div. 91 (1965) 225–247. doi: [https://doi.org/10.1061/JYCEAJ.0001277 10.1061/JYCEAJ.0001277].

<span id='_Bib057'></span>
[57] G. Parker, P.C. Klingeman, D.G. McLean, Bedload and size distribution in paved gravel-bed streams., J. Hydraul. Div. - ASCE. (1982). doi: [https://doi.org/10.1061/(ASCE)0733-9429(1983)109:5(793) 10.1061/(ASCE)0733-9429(1983)109:5(793)].

<span id='_Bib058'></span>
[58] P.R. Wilcock, J.C. Crowe, Surface-based Transport Model for Mixed-Size Sediment, J. Hydraul. Eng. 129 (2003) 120–128. doi:  [https://doi.org/10.1061/(ASCE)0733-9429(2003)129:2(120) 10.1061/(ASCE)0733-9429(2003)129:2(120)].