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

Sintesi Tecnica: Dispersione Longitudinale in Ambienti Fluviali

Definizione del Problema
Una valutazione affidabile del trasporto di inquinanti e dei processi di miscelazione è fondamentale per la gestione della qualità dell'acqua e la conservazione ambientale nei fiumi. Il coefficiente di dispersione longitudinale ($K$) è un parametole cruciale nell'equazione di advezione-dispersione unidimensionale, che governa l'accuratezza delle simulazioni del trasporto di inquinanti. Tuttavia, le formulazioni empiriche esistenti per $K$ si basano spesso su set di dati che combinano canali sia stretti che larghi, rischiando di trascurare i distinti meccanismi di dispersione inerenti alle diverse geometrie dei canali. Inoltre, nessuna formulazione esistente è stata dimostrata capace di operare con costanza in modo efficace sia su dataset sperimentali che su dati sul campo per diversi scenari di rilascio (istantaneo vs continuo). Questo studio affronta la necessità di formulazioni specifiche per regime e un'analisi dettagliata di come la geometria del canale influenzi la dinamica del trasporto di inquinanti.

Metodologia
La ricerca impiega un approccio su due fronti: lo sviluppo di nuove formulazioni empiriche e la simulazione numerica del trasporto di inquinanti.

1. Sviluppo della Formulazione Empirica:

   • Lo studio utilizza l'analisi dimensionale per esprimere il coefficiente di dispersione adimensionale ($\tilde{K} = K/uh$) in funzione del rapporto della velocità di taglio ($\Xi = u_*/u$) e del rapporto di aspetto del canale ($\beta = w/h$).
   • Viene adottato un approccio basato sul regime, classificando i canali come "stretti" ($\beta \leq 10$) e "larghi" ($\beta > 10$) sulla base dei trend osservati nei dataset derivati da rilievi sul campo.
   • L'analisi di regressione è stata eseguita su 105 grafici di dati sul campo provenienti da varie fonti della letteratura (ad es., Sooky, Godfrey e Frederick, Yotsukura et al.) per determinare i coefficienti e gli esponenti della relazione di legge di potenza $\tilde{K} = \mathcal{C} \Xi^a \beta^b$.
   • Parametri distinti sono stati derivati per canali stretti e larghi per catturare i loro unici comportamenti idraulici.

2. Simulazione Numerica:

   • L'equazione di advezione-dispersione unidimensionale è stata risolta numericamente utilizzando un metodo di differenza finita.
   • Il termine di advezione è stato discretizzato utilizzando lo schema esplicito di Lax-Wendroff (secondo ordine di accuratezza), mentre il termine di diffusione ha utilizzato la differenza centrale con avanzamento temporale di Euler forward (primo ordine di accuratezza).
   • La stabilità è stata garantita tramite l'analisi di von Neumann, aderendo ai criteri $\lambda_1 + 2\lambda_2 \leq 1$.
   • Le simulazioni sono state condotte per due scenari:
     • Rilascio Istantaneo: Modellato come una funzione gaussiana che rappresenta sversamenti accidentali.
     • Rilascio Continuo: Modellato con una concentrazione della sorgente sostenuta per un intervallo di tempo specifico.
   • L'analisi di sensibilità della griglia ha determinato un passo spaziale ($\Delta x$) di 1 m per l'indipendenza dalla dimensione della griglia.

Contributi Chiave

• Nuove Formulazioni Specifiche per Regime: L'articolo propone nuove equazioni empiriche per il coefficiente di dispersione longitudinale specificamente adattate per canali stretti e larghi, anziché una singola equazione unificata.
• Identificazione dei Parametri Dominanti: L'analisi rivela che la dispersione longitudinale nei canali stretti è governata principalmente dal rapporto di aspetto del canale ($\beta$), mentre nei canali larghi è governata dal rapporto della velocità di taglio ($\Xi$).
• Definizione della Soglia: Lo studio stabilisce un rapporto di aspetto del canale pari a 10 come soglia per demarcare i regimi di canale stretto e largo, un valore scelto per allinearsi meglio ai trend dei dati sul campo rispetto a precedenti soglie (ad es., 28 o 30,6).
• Dinamica Completa del Trasporto: Lo studio risolve numericamente l'equazione di advezione-dispersione per esaminare l'evoluzione spazio-temporale delle concentrazioni di inquinanti per rilasci sia istantanei che continui attraverso diverse geometrie di canale.

Risultati

• Comportamento del Coefficiente di Dispersione:
  • Nei canali stretti, il coefficiente di dispersione adimensionale ($\tilde{K}$) aumenta all'aumentare del rapporto della velocità di taglio ($\Xi$), guidato dalla dispersione indotta dal taglio e dal fenomeno del calo di velocità.
  • Nei canali larghi, $\tilde{K}$ diminuisce all'aumentare di $\Xi$. Ciò è attribuito alla dominanza della miscelazione turbolenta trasversale, che promuove una rapida omogeneizzazione della concentrazione dell'inquinante attraverso la sezione trasversale, riducendo così la diffusione longitudinale.
  • È stato riscontrato che le formulazioni esistenti non riescono a catturare questi comportamenti contrastanti, spesso sovrastimando o sottostimando $\tilde{K}$ a seconda del regime del canale.
• Dinamica del Trasporto di Inquinanti:
  • Canali Stretti (Istantaneo): I profili di concentrazione mostrano un picco distinto che si attenua man mano che il pennacchio si sposta a valle.
  • Canali Stretti (Continuo): La concentrazione rimane alla magnitudo massima per un intervallo temporale e spaziale finito, formando un ampio plateau senza decadimento significativo, indicando la dominanza dell'advezione sulla dispersione.
  • Canali Larghi (Istantaneo): La dispersione potenziata porta a un'attenuazione intensificata del picco di concentrazione rispetto ai canali stretti. I profoli mostrano un ramo ascendente rapido e un ramo discendente più lungo.
  • Canali Larghi (Continuo): A differenza dei canali stretti, i canali larghi mostrano un decadimento significativo della concentrazione di picco nella direzione a valle a causa della dispersione aumentata.
• Validazione:
  • Le formulazioni proposte hanno mostrato un forte accordo con le osservazioni sul campo, con la maggior parte dei punti dati che rientra nelle bande di errore di $\pm 20\%$.
  • Le previsioni numeriche per il rilascio istantaneo si sono allineate bene con le soluzioni analitiche classiche, sebbene siano state osservate concentrazioni di picco leggermente inferiori a causa della diffusione numerica.
  • I confronti con i dati dei traccianti sul campo del fiume Severn (Atkinson e Devis, 2000) hanno mostrato un accordo generale, sebbene le osservazioni sul campo abbiano mostrato picchi leggermente inferiori e un trasporto più veloce, probabilmente a causa di fattori non modellati come l'accumulo transitorio e gli apporti laterali.

Significatività e Rivendicazioni
Gli autori affermano che questo studio offre un quadro più affidabile per stimare il coefficiente di dispersione longitudinale, riconoscendo che un'unica equazione universale è insufficiente per tutte le geometrie fluviali. Distinguendo tra regimi stretti e larghi, le formulazioni proposte forniscono una rappresentazione più accurata della fisica della dispersione. Inoltre, lo studio elucida le distinte dinamiche di trasporto degli inquinanti sotto diversi scenari di rilascio e geometrie di canale, dimostrando che l'interazione tra advezione e dispersione varia significativamente tra canali stretti e larghi. Il lavoro funge da strumento robusto per simulare il trasporto di inquinanti, sebbene gli autori notino limitazioni riguardanti l'esclusione della curvatura del fiume, dello scambio con le pianure alluvionali, della resistenza della vegetazione e dei processi di trasporto reattivo, suggerendo questi come aree per indagini future.