Un modello 1D-2D di acque basse per il calcolo della propagazione di piene naturali

Abstract

Le equazioni 2D del Saint-Venant (SV), o delle acque basse (shallow water), si utilizzano ampiamente per le simulazioni idrodinamiche in fiumi, laghi, zone di estuari e zone golenali. Come ampiamente descritto in letteratura, in base ai termini che si trascurano nell'equazione del momento, si possono distinguere modelli completi, cinematici, diffusivi e a zero inerzia convettiva. Escludendo eventi estremi di shock, nella normalità dei casi relativi a fenomeni di inondazione, si possono trascurare i termini inerziali nell'equazione del momento, ottenendo il modello diffusivo. La scelta del modello da utilizzare dipende anche dalla disponibilità e dall'accuratezza dei dati di input (Smith et al., 2004), nonché dall'accuratezza richiesta ai risultati dal problema gestionale che motiva l'uso del modello. I dati di input sono per esempio la topografia, le proprietà idrauliche dei tratti di fiume e delle zone golenali, la conoscenza degli idrogrammi di ingresso. La soluzione del modello completo è molto sensibile all'errore commesso nella stima delle quote topografiche (Aricò et al., 2011), mentre il modello cinematico non è in grado di simulare processi di backwater. Nel presente lavoro si propone un modello diffusivo che è particolarmente adatto alla simulazione degli eventi di piena, quando il dominio di calcolo è costituito dalla parte centrale del canale, in cui si realizza una propagazione di tipo sostanzialmente monodimensionale, nonché dalle fasce laterali, dove la corrente può divagare secondo uno schema tipicamente bidimensionale. Il modello proposto è di tipo predizione-correzione e garantisce la conservazione della massa, sia localmente che globalmente. Le equazioni di governo sono risolte su mesh miste, composte da domini 2D, discretizzati con elementi triangolari, e da domini 1D discretizzati da elementi quadrilateri. Nei domini 1D si assume una quota piezometrica costante nella direzione normale al flusso, nonché la relazione fra il gradiente piezometrico e la portata data dalla legge di resistenza propria delle sezioni di forma complessa. I flussi scambiati lateralmente con il dominio 2D sono calcolati invece in base alla discretizzazione spaziale degli elementi triangolari. Ciò consente di calcolare il flusso principale, nella direzione dell'alveo, tenendo conto mediante la legge di resistenza anche degli sforzi turbolenti causati dai vortici con asse verticale all'interno dell'alveo, sforzi che altrimenti verrebbero trascurati anche applicando un tradizionale modello completo 2D di acque basse. Utilizzando la schematizzazione suddetta è altresì possibile evitare un drastico aumento della densità della mesh all'interno dell'alveo, altrimenti necessario per una corretta rappresentazione del fondo. Un ulteriore vantaggio dell'approccio diffusivo consiste nel fatto che l'onere computazionale risulta poco più che proporzionale al numero di elementi di calcolo e non vi è per la stabilità delle soluzione alcuna restrizione del passo temporale, in base al massimo numero di Courant calcolato nel dominio di calcolo. L'approccio diffusivo consente infine di trattare facilmente le discontinuità esistenti lungo gli elementi 1D a causa di briglie o salti, ovvero a causa di pareti laterali sub-verticali. Il sistema di equazioni di governo alle derivate parziali si risolve ad ogni livello temporale frazionandolo in un sistema di predizione ed uno di correzione. Si dimostra che il sistema di predizione ha le caratteristiche funzionali di un problema puramente convettivo, mentre quello di correzione ha le caratteristiche funzionali di un problema puramente diffusivo. L'input del modello è costituito, per il dominio 1D, dalla geometria delle sezioni trasversali degli elementi 1D, assegnate mediante serie di punti quotati, nonché dalle linee di livello e/o da singoli punti quotati per quanto concerne il dominio 2D. Il software denominato NETGEN (Schöberl, 1997), genera quind