Longitudinal dispersion in fluvial environments: new empirical formulations and applications to pollutant transport modeling

Resumen Técnico: Dispersión Longitudinal en Entornos Fluviales

Planteamiento del Problema
La evaluación confiable del transporte de contaminantes y de los procesos de mezcla es crítica para la gestión de la calidad del agua y la conservación ambiental en los ríos. El coeficiente de dispersión longitudinal ($K$) es un parámetro fundamental en la ecuación de advección-dispersión unidimensional, el cual rige la precisión de las simulaciones de transporte de contaminantes. Sin embargo, las formulaciones empíricas existentes para $K$ a menudo dependen de conjuntos de datos que combinan canales tanto estrechos como anchos, lo que podría pasar por alto los distintos mecanismos de dispersión inherentes a diferentes geometrías de canales. Además, ninguna formulación existente ha demostrado desempeñarse de manera consistente tanto en datos experimentales como de campo para diversos escenarios de liberación (instantánea frente a continua). Este estudio aborda la necesidad de formulaciones específicas por régimen y un análisis detallado de cómo la geometría del canal influye en la dinámica del transporte de contaminantes.

Metodología
La investigación emplea un enfoque de dos vertientes: el desarrollo de nuevas formulaciones empíricas y la simulación numérica del transporte de contaminantes.

1. Desarrollo de la Formulación Empírica:

   • El estudio utiliza el análisis dimensional para expresar el coeficiente de dispersión adimensional ($\tilde{K} = K/uh$) en términos de la relación de velocidad de corte ($\Xi = u_*/u$) y la relación de aspecto del canal ($\beta = w/h$).
   • Se adopta un enfoque basado en regímenes, clasificando los canales como "estrechos" ($\beta \leq 10$) y "anchos" ($\beta > 10$) basándose en las tendencias observadas en los conjuntos de datos derivados de campo.
   • Se realizó un análisis de regresión sobre 105 gráficos de datos de campo de diversas fuentes bibliográficas (p. ej., Sooky, Godfrey y Frederick, Yotsukura et al.) para determinar los coeficientes y exponentes de la relación de ley de potencia $\tilde{K} = \mathcal{C} \Xi^a \beta^b$.
   • Se derivaron parámetros distintos para canales estrechos y anchos para capturar sus comportamientos hidráulicos únicos.

2. Simulación Numérica:

   • La ecuación de advección-dispersión unidimensional se resolvió numéricamente utilizando un método de diferencias finitas.
   • El término de advección se discretizó utilizando el esquema Lax-Wendroff (de segundo orden de precisión), mientras que el término de difusión utilizó diferencias centrales con paso de tiempo de Euler hacia adelante (de primer orden de precisión).
   • La estabilidad se aseguró mediante el análisis de von Neumann, adhiriéndose a los criterios $\lambda_1 + 2\lambda_2 \leq 1$.
   • Las simulaciones se realizaron para dos escenarios:
     • Liberación Instantánea: Modelada como una función gaussiana que representa derrames accidentales.
     • Liberación Continua: Modelada con una concentración de fuente sostenida durante un intervalo de tiempo específico.
   • El análisis de sensibilidad de la malla determinó un tamaño de paso espacial ($\Delta x$) de 1 m para la independencia del tamaño de la malla.

Contribuciones Clave

• Nuevas Formulaciones Específicas por Régimen: El artículo propone nuevas ecuaciones empíricas para el coeficiente de dispersión longitudinal diseñadas específicamente para canales estrechos y anchos, en lugar de una ecuación unificada.
• Identificación de Parámetros Gobernantes: El análisis revela que la dispersión longitudinal en canales estrechos está gobernada principalmente por la relación de aspecto del canal ($\beta$), mientras que en canales anchos, está gobernada por la relación de velocidad de corte ($\Xi$).
• Definición de Umbral: El estudio establece una relación de aspecto de canal de 10 como el umbral para demarcar los regímenes de canales estrechos y anchos, un valor elegido para alinearse mejor con las tendencias de los datos de campo que los umbrales previos (p. ej., 28 o 30.6).
• Dinámica de Transporte Integral: El estudio resuelve numéricamente la ecuación de advección-dispersión para examinar la evolución espaciotemporal de las concentraciones de contaminantes tanto para liberaciones instantáneas como continuas a través de diferentes geometrías de canales.

Resultos

• Comportamiento del Coeficiente de Dispersión:
  • En canales estrechos, el coeficiente de dispersión adimensional ($\tilde{K}$) aumenta a medida que aumenta la relación de velocidad de corte ($\Xi$), impulsado por la dispersión inducida por el cizallamiento y el fenómeno de caída de velocidad.
    la canales anchos, $\tilde{K}$ disminuye a medida que $\Xi$ aumenta. Esto se atribuye al predominio de la mezcla turbulenta transversal, que promueve una homogeneización rápida de la concentración del contaminante a través de la sección transversal, reduciendo así la propagación longitudinal.
  • Se encontró que las formulaciones existentes fallan al capturar estos comportamientos contrastantes, sobreprediciendo o subprediciendo a menudo $\tilde{K}$ dependiendo del régimen del canal.
• Dinámica de Transporte de Contaminantes:
  • Canales Estrechos (Instantánea): Los perfiles de concentración exhiben un pico distintivo que se atenúa a medida que la pluma viaja río abajo.
  • Canales Estrechos (Continua): La concentración permanece en su magnitud máxima durante un intervalo temporal y espacial finito, formando una meseta amplia sin una degradación significativa, lo que indica un predominio de la advección sobre la dispersión.
  • Canales Anchos (Instantánea): La dispersión mejorada conduce a una atenuación intensificada del pico de concentración en comparación con los canales estrechos. Los perfiles muestran un flanco de ascenso rápido y un flanco de caída más largo.
  • Canales Anchos (Continua): A diferencia de los canales estrechos, los canales anchos exhiben una degradación significativa del pico de concentración en la dirección de la corriente debido a la dispersión aumentada.
• Validación:
  • Las formulaciones propuestas mostraron una fuerte concordancia con las observaciones de campo, con la mayoría de los puntos de datos dentro de las bandas de error de $\pm 20\%$.
  • Las predicciones numéricas para la liberación instantánea se alinearon bien con las soluciones analíticas clásicas, aunque se observaron concentraciones pico ligeramente menores debido a la difusión numérica.
  • Las comparaciones con los datos de trazadores de campo del río Severn (Atkinson y Devis, 2000) mostraron una concordancia general, aunque las observaciones de campo exhibieron picos ligeramente menores y un transporte más rápido, probablemente debido a factores no modelados como el almacenamiento transitorio y las entradas laterales.

Significancia y Reivindicaciones
Los autores afirman que este estudio ofrece un marco más confiable para estimar el coeficiente de dispersión longitudinal al reconocer que una única ecuación universal es insuficiente para todas las geometrías de ríos. Al distinguir entre regímenes estrechos y anchos, las formulaciones propuestas proporcionan una representación más precisa de la física de la dispersión. Además, el estudio elucida las distintas dinámicas de transporte de contaminantes bajo diferentes escenarios de liberación y geometrías de canales, demostrando que la interacción entre la advección y la dispersión varía significamente entre canales estrechos y anchos. El trabajo sirve como una herramienta robusta para simular el transporte de contaminantes, aunque los autores señalan las limitaciones respecto a la exclusión de la curvatura del río, el intercambio con llanuras de inundación, la resistencia de la vegetación y los procesos de transporte reactivo, sugiriéndolos como áreas para futuras investigaciones.