Relatori: Prof. Cristiana Di Cristo Prof. Rui Ferreira (IST, Lisbona)

1. UNIVERSITÀ DEGLI STUDI DI CASSINO E DEL LAZIO MERIDIONALE Dipartimento Di Ingegneria Civile e Meccanica Tesi di Laurea Magistrale in Ingegneria Civile VORTEX DYNAMICS IN FLOWS WITH EMERGENT VEGETATION Candidato: Simona Di Carlo Matr. 0033453 Relatori: Prof. Cristiana Di Cristo Prof. Rui Ferreira (IST, Lisbona)

2. La vegetazione emergente e le Wetlands La vegetazione acquatica in generale riveste un ruolo di primaria importanza negli ambienti costieri e fluviali. In particolare, la vegetazione a fusto emergente rappresenta uno dei principali costituenti delle zone umide(wetlands).

3. Le Wetlands Le wetland sono ecosistemi di transizione e di collegamento tra terra ed acqua. Tali ecosistemi rivestono notevole importanza sotto diversi punti di vista: • Idrogeologico • Chimico e fisico • Biologico

4. Obiettivi Lo scopo di questo lavoro di tesi è quello di fornire un contributo allo studio sulla gestione e ripristino delle wetland, accrescendo la conoscenza dei processi fisici inerenti ad esse. E’ stata condotta una campagna sperimentale per investigare gli aspetti cinematici relativi al fenomeno di distacco dei vortici da cilindri rigidi e circolari che simulano la vegetazione emergente. La presenza dei vortici all’interno dei corpi idrici influenza i processi di sedimentazione e dispersione.

5. Apparato sperimentale PIV System (Particle Image Velocimetry) Fonte di luce laser pulsato Generatore Telecamera Unità di sincronizzazione Software di controllo ed acquisizione

6. Prove sperimentali effettuate Test A e A* Test B Test D . Test C .

7. Prove sperimentali effettuate Caratteristiche idrauliche dei test effettuati

8. Prove sperimentali effettuate Velocità medie temporali Test B Test A Test C Test D Test A*

9. Metodologia Il presente lavoro di tesi introduce un originale metodo, implementato in linguaggio Matlab, per la rilevazione di vortici in campi di moto bidimensionali. 1) Identificazione dei vortici 2) Monitoraggio dei vortici d 3d 4d

10. Metodologia 1) Identificazione dei vortici Criterio delle linee isovelocità

11. Metodologia 1) Identificazione dei vortici Criterio delle linee isovelocità Criterio delle linee di flusso chiuse

12. Metodologia 1) Identificazione dei vortici Criterio delle linee isovelocità Criterio delle linee di flusso chiuse Vortici positivi Vortici negativi

13. Metodologia 2) Monitoraggio dei vortici Vortici positivi Vortici negativi • La traiettoriadeivortici è definita da almeno due vorticicorrelatitraloro. • Le traiettorie positive sonoformate da vorticipositivi e le traiettorie negative sonoformate da vorticinegativi.

14. Risultati dell’analisi statistica dei dati L’analisi dei dati effettuata è stata di tipo statistico. • Frequenza media di distacco dei vortici • Numero di Strouhal • Numero normalizzato di vortici • Distribuzione longitudinale • Survival rate • Distribuzione laterale • Traiettorie medie e relativa deviazione standard

15. Risultati dell’analisi statistica dei dati • Frequenza media di distacco dei vortici • Numero di Strouhal fs : frequenza di distacco dei vortici L : lunghezza caratteristica V : velocità media di portata Il meccanismo di distacco dei vortici è molto simile in tutti i test indagati.

16. Risultati dell’analisi statistica dei dati • Numero normalizzato di vortici ND : numero di vortici identificati NS: numero di vortici che si formano nell’intervallo temporale di acquisizione NV : numero ideale di visualizzazioni del singolo vortice Il massimo numero normalizzato di vortici è stato riscontrato nel test D.

17. Risultati dell’analisi statistica dei dati • Distribuzione longitudinale dei vortici Un’ ulteriore analisi è stata effettuata al fine di individuare la durata e la collocazione spaziale dei vortici lungo la direzione del moto. Test C Vortici positivi Vortici negativi La distribuzione mostra un andamento approsimativamente lineare

18. Risultati dell’analisi statistica dei dati • Distribuzione longitudinale dei vortici Un’ ulteriore analisi è stata effettuata al fine di individuare la durata e la collocazione spaziale dei vortici lungo la direzione del moto. Test A Vortici positivi Vortici negativi Test statistico di ‘bontà di adattamento’ di Kolmogorov – Smirnov La distribuzione longitudinale dei vortici ben si adatta alla distribuzione teorica di probabilità di Burr.

19. Risultati dell’analisi statistica dei dati • Survival rate Il survival rate indica il tasso di sopravvivenza dei vortici lungo la direzione del moto. Vortici positivi Vortici negativi I risultati di tale analisi trovano spiegazione nell’interazione con la turbolenza di fondo.

20. Risultati dell’analisi statistica dei dati • Distribuzione laterale dei vortici Test B Area A1 Area A2 La distribuzione dei vortici lungo la direzione trasversale al moto è stata indagata per comprendere meglio la collocazione spaziale dei vortici. Vortici positivi Vortici negativi Le distribuzioni laterali seguono la distribuzione di probabilità di Burr.

21. Risultati dell’analisi statistica dei dati • Traiettorie medie dei vortici e deviazione standard

22. Risultati dell’analisi statistica dei dati • Traiettorie medie dei vortici e deviazione standard Test B Test A Test A* Test C Vortici positivi Vortici negativi Test D Vortici positivi Vortici negativi

23. Risultati dell’analisi statistica dei dati • Traiettorie medie dei vortici e deviazione standard La maggiore lunghezza media delle traiettorie è stata individuata nel test D.

24. Conclusioni • Il meccanismo di distacco dei vortici non è influenzato dalla particolare configurazione geometrica dei test oggetti di studio: ciò è evidenziato dai valori del Numero di Strouhal, che risultano essere mediamente uguali per tutti i test effettuati; • Il valore medio del numero di Strouhal è di St≈0.26, poco più alto del valore St=0.21 noto in letteratura per il caso del cilindro isolato (numero di Reynolds Rep=1000). Questa differenza potrebbe essere dovuta a macroscabrezze presenti sui cilindri; • Il numero di vortici rilevati positivi e negativi è mediamente lo stesso; • L’analisi relativa al survival rate ha dimostrato che nei test con distribuzione densa di cilindri i vortici perdono la loro forma molto più rapidamente rispetto ai test caratterizzati da densità minori. Il test D, ovvero il caso del cilindro isolato situato in prossimità di un gruppo di cilindri, mostra il maggiore valore del survival rate;

25. Conclusioni • Il test D presenta la maggiore lunghezza media delle traiettorie; • I test A ed A* presentano delle distribuzioni longitudinali e trasversali entrambe pienamente adattabili alla distribuzione di probabilità teorica di Burr; • In definitiva è quindi possibile affermare che tutti i risultati ottenuti sono perfettamente in accordo tra di essi. L’ assenza di cilindri a monte del singolo cilindro indagato determina un maggiore numero di vortici e una maggiore durata spaziale e temporale di essi. • Questo comportamento è stato interpretato con una ulteriore indagine sull’interazione tra i vortici e la turbolenza di fondomediante l’analisi delle mappe relative all’ intensità di turbolenza.

26. GRAZIE PER L’ATTENZIONE SIMONA DI CARLO • Il lavoro sperimentale di questa tesi è stato anche finanziato dalla ‘Portuguese Foundation for Science and Technology’ nell’ ambito del progetto di ricerca PTDC/ECM/11760/2010.

27. Kineticturbulentenergy production Una possibile spiegazione relativa all’analisi condotta sul survival rate può essere trovata in termini di produzione di energia cinetica turbolenta. Una successiva analisi ha infatti mostrato che il termine PTKE risulta essere maggiore di zero nei test A ed A* ed invece risulta essere praticamente nullo nei test B, C e D. Si può quindi assumere che è un problema di turbolenza di fondo. Questo si vede in maniera molto chiara nei grafici mostrati: Where i and j are indices running from 1 to 3; xj stands for three cartesian spatial directions; is the product of the time averaged velocity fluctuations in the directions i and j ; is the time averaged velocity in the direction j. Maps of (on the left hand side) and (on the right hand side). The colormaps are represented in SI units (m2s-2) and the arrow indicates the flow direction

28. Kineticturbulentenergy production Test A Test C Test B