Solution of the Advection-Diffusion Equation in Two Dimensions using the Deposition Velocity

Resumen Técnico: Solución de la Ecuación de Advección-Difusión en Dos Dimensiones utilizando la Velocidad de Deposición

Planteamiento del Problema
Este estudio aborda la solución analítica de la ecuación de advección-difusión (ADE) estacionaria en dos dimensiones bajo condiciones atmosféricas inestables. El objetivo principal es modelar la dispersión de contaminantes, centrándose específicamente en las concentraciones normalizadas integradas en el sentido transversal al viento en la línea central. El modelo incorpora una condición de contorno de velocidad de deposición en el nivel del suelo y utiliza funciones de ley de potencia tanto para la difusividad turbulenta vertical como para la velocidad del viento en función de la altura vertical. El estudio busca validar este enfoque analítico frente a los datos observados de concentración de hexafluoruro de azufre (SF6) del experimento de campo de Copenhague en Dinamarca.

Metodología
Los autores emplean la técnica de separación de variables para resolver la ADE gobernante. El marco matemático se define de la siguiente manera:

• Ecuación Gobernante: Se formula la ADE en dos dimensiones donde la velocidad del viento $u(z)$ y la difusividad turbulenta vertical $k_z(z)$ se expresan como leyes de potencia de la altura vertical ($z^p$ y $z^q$, respectivamente).
• Condiciones de Contorno:
  1. Límite Superior (Altura de Mezcla $h$): Condición de flujo cero ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Límite del Suelo ($z=0$): Una condición de flujo de deposición donde el flujo es igual al producto de la velocidad de deposición ($V_d$) y la concentración ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Término de Fuente: Una fuente puntual continua a una altura de chimenea $h_s$ representada por una función delta de Dirac.
• Derivación Analítica:
  • La ecuación se separa en componentes espaciales $X(x)$ y $Z(z)$.
  • El componente vertical conduce a una ecuación diferencial de tipo Sturm-Liouville con coeficientes variables.
  • Al asumir valores específicos para los exponentes (por ejemplo, $q=0.85$) y utilizar una transformación de variable, la ecuación se reduce a una ecuación de Bessel generalizada.
  • La solución para el perfil de concentración vertical se deriva utilizando funciones de Bessel de primera y segunda especie, que se simplifican a formas trigonométricas (seno y coseno) para funciones de Bessel de orden semientero.
  • Los autovalores y autofunciones se determinan aplicando las condiciones de contorno, lo que conduce a una solución en serie para la concentración $C_y(x, z)$.
  • Los coeficientes se calculan utilizando las condiciones de ortogonalidad derivadas del perfil de velocidad del viento.

Contribuciones Clave

• Marco Analítico: El artículo presenta una solución analítica de forma cerrada para la ADE en 2D bajo condiciones inestables, incorporando explícitamente una condición de contorno de velocidad de deposición, lo cual es un factor que suele ser complejo en la modelización de la dispersión.
• Parametrización por Ley de Potencia: El modelo utiliza representaciones de ley de potencia para la velocidad del viento y la difusividad turbulenta, lo que permite flexibilidad a través de diferentes clases de estabilidad (definidas por los exponentes $p$ y $q$).
• Validación: El modelo analítico es rigurosamente probado contra el conocido conjunto de datos de SF6 de Copenhague, cubriendo diversas clases de estabilidad (A, B, C, D) y distancias a sotavento.

Resultados
El estudio compara las concentraciones normalizadas integradas en el sentido transversal al viento predichas con las concentraciones de SF6 observadas en el experimento de Copenhague.

• Comparación de Datos: La Tabla 2 y la Figura 1 presentan una comparación directa de las concentraciones observadas frente a las predichas en múltiples ejecuciones (Run 1 a Run 9) a distintas distancias a sotavento.
• Desempeño Estadístico: El desempeño del modelo se evalúa mediante métricas estadísticas estándar (Hanna, 1989):
  • Error Cuadrático Medio Normalizado (NMSE): 0.07 (indicando un error bajo).
  • Sesgo Fraccional (FB): 0.003 (indicando un sesgo insignificante).
  • Coeficiente de Correlación (COR): 0.82.
  • Fracción dentro de un Factor de Dos (FAC2): 0.997.
• Acuerdo Visual: La Figura 1 ilustra un alto grado de acuerdo entre los perfiles de concentración observados y predichos. La Figura 2 confirma que los valores predichos caen dentro de un factor de dos de los datos observados.

Significancia y Reivindicaciones
Los autores afirman que el modelo analítico propuesto logra un "alto grado de acuerdo cerrado" con los datos de campo observados. Específicamente, el artículo establece que las concentraciones predichas alcanzaron un 100% de acuerdo con las concentraciones observadas en la línea central, como lo evidencia el valor de FAC2 de 0.997 y el sesgo cercano a cero. El estudio concluye que el método de separación, combinado con las parametrizaciones de ley de potencia para el viento y la difusividad, proporciona una herramienta analítica robusta para predecir la dispersión de contaminantes bajo condiciones atmosféricas inestables sin la necesidad de simulaciones numéricas complejas. El trabajo refuerza la utilidad de las soluciones analíticas basadas en transformaciones y separación para capturar la física esencial del transporte de contaminantes en capas límite inestables.