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

Sintesi Tecnica: Soluzione dell'Equazione di Advezione-Diffusione in Due Dimensioni utilizzando la Velocità di Deposizione

Enunciato del Problema
Questo studio affronta la soluzione analitica dell'equazione di advezione-diffusione (ADE) stazionaria in due dimensioni sotto condizioni atmosferiche instabili. L'obiettivo primario è modellare la dispersione degli inquinanti, concentrandosi specificamente sulle concentrazioni normalizzate integrate trasversalmente al vento lungo la linea centrale. Il modello incorpora una condizione al contorno di velocità di deposizione al livello del suolo e utilizza funzioni di legge di potenza sia per la diffusività di eddy verticale che per la velocità del vento in funzione dell'altezza verticale. Lo studio mira a validare questo approccio analitico rispetto ai dati osservati di concentrazione di esanofluoruro di zolfo (SF6) provenienti dall'esperimento sul campo di Copenhagen, in Danimarca.

Metodologia
Gli autori utilizzano la tecnica della separazione delle variabili per risolvere l'ADE governante. Il quadro matematico è definito come segue:

• Equazione Governante: L'ADE bidimensionale è formulata dove la velocità del vento $u(z)$ e la diffusività di eddy verticale $k_z(z)$ sono espresse come leggi di potenza dell'altezza verticale ($z^p$ e $z^q$, rispettivamente).
• Condizioni al Contorno:
  1. Confine Superiore (Altezza di Miscelamento $h$): Condizione di flusso nullo ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Confine al Suolo ($z=0$): Una condizione di flusso di deposizione dove il flusso è uguale al prodotto della velocità di deposizione ($V_d$) e della concentrazione ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Termine Sorgente: Una sorgente puntiforme continua all'altezza del camino $h_s$ rappresentata da una funzione delta di Dirac.
• Derivazione Analitica:
  • L'equazione viene separata in componenti spaziali $X(x)$ e $Z(z)$.
  • La componente verticale porta a un'equazione differenziale di tipo Sturm-Liouville con coefficienti variabili.
  • Assumendo valori specifici per gli esponenti (ad esempio, $q=0.85$) e utilizzando una trasformazione di variabile, l'equazione viene ridotta a un'equazione di Bessel generalizzata.
  • La soluzione per il profilo di concentrazione verticale è derivata utilizzando le funzioni di Bessel di prima e seconda specie, che si semplificano in forme trigonometriche (seno e coseno) per le funzioni di Bessel di ordine semintero.
  • Autovalori e autofunzioni sono determinati applicando le condizioni al contorno, portando a una soluzione in serie per la concentrazione $C_y(x, z)$.
  • I coefficienti sono calcolati utilizzando le condizioni di ortogonalità derivate dal profilo della velocità del vento.

Contributi Chiave

• Quadro Analitico: L'articolo presenta una soluzione analitica a forma chiusa per l'ADE 2D in condizioni instabili, incorporando esplicitamente una condizione al contorno di velocità di deposizione, che è spesso un fattore complesso nella modellazione della dispersione.
• Parametrizzazione per Legge di Potenza: Il modello utilizza rappresentazioni per legge di potenza per la velocità del vento e la diffusività di eddy, permettendo flessibilità attraverso diverse classi di stabilità (definite dagli esponenti $p$ e $q$).
• Validazione: Il modello analitico è rigorosamente testato contro il noto dataset SF6 di Copenhagen, coprendo varie classi di stabilità (A, B, C, D) e distanze a valle.

Risultati
Lo studio confronta le concentrazioni normalizzate integrate trasversalmente al vento previste con le concentrazioni di SF6 osservate nell'esperimento di Copenhagen.

• Confronto dei Dati: La Tabella 2 e la Figura 1 presentano un confronto diretto tra le concentrazioni osservate e quelle previste attraverso molteplici run (Run 1 attraverso Run 9) a diverse distanze a valle.
• Performance Statistica: La performance del modello è valutata utilizzando metriche statistiche standard (Hanna, 1989):
  • Errore Quadratico Medio Normalizzato (NMSE): 0.07 (che indica un errore basso).
  • Bias Frazionario (FB): 0.003 (che indica un bias trascurabile).
  • Coefficiente di Correlazione (COR): 0.82.
  • Valore entro un Fattore di Due (FAC2): 0.997.
• Accordo Visivo: La Figura 1 illustra un alto grado di accordo tra i profili di concentrazione osservati e previsti. La Figura 2 conferma che i valori previsti rientrano in un fattore di due rispetto ai dati osservati.

Significatività e Rivendicazioni
Gli autori affermano che il proposto modello analitico raggiunge un "alto grado di accordo chiuso" con i dati misurati sul campo. Nello specifico, l'articolo dichiara che le concentrazioni previste hanno ottenuto un accordo del 100% con le concentrazioni osservate sulla linea centrale, come evidenziato dal valore FAC2 di 0.997 e dal quasi nullo bias. Lo studio conclude che il metodo di separazione, combinato con le parametrizzazioni per legge di potenza per vento e diffusività, fornisce uno strumento analitico robusto per prevedere la dispersione degli inquinanti in condizioni atmosferiche instabili senza la necessità di complesse simulazioni numeriche. Il lavoro rinforza l'utilità delle soluzioni analitiche basate su trasformate e separazione nel catturare la fisica essenziale del trasporto di inquinanti negli strati limite instabili.