Dissertations / Theses on the topic 'SCUOLA DI DOTTORATO DI ENVIRONMENTAL AND INDUSTRIAL FLUID MECHANICS'

2010/2011
In the past two decades, the growing interest in the study of fluid flows involving discontinuities, such as shocks or high gradients, where a quadratic-convergent method may not provide a satisfactory solution, gave a notable impulse to the employment of high-order techniques. The present dissertation comprises the analysis and numerical testing of two high-order methods. The first one, belonging to the discontinuous finite-element class, is the discontinuous control-volume/finite-element method (DCVFEM) for the advection/ diffusion equation. The second method refers to the high-order finite-difference class, and is the mixed weighted non-oscillatory scheme (MWCS) for the solution of the compressible Euler equations. The methods are described from a formal point of view, a Fourier analysis is used to assess the dispersion and dissipation errors, and numerical simulations are conducted to confirm the theoretical results.
XXIV Ciclo
1980

2009/2010
Fluid structure interaction (FSI) is one of fundamental phenomena encountered everywhere in nature or in industrial systems as well as one of the most studied and the most challenging topics in the fluid mechanics. Its research presents the core objective of this dissertation, along with experimental study of artificial heart devices. Better understanding of FSI could turn the still unexploited phenomenon into a powerful tool for resolving wealthy of multi-physics problems. Recently computational fluid dynamics community has been putting enormous efforts to uncover, make clear and answer yet numerous issues related to this developing topic. In addition, the FSI is often followed by the vortex formation, one more phenomena that could be both powerful driving force as well as distracting, disturbing occurrence. Consequently, this dissertation will begin with addressing some open issues related to the fluid-structure interaction associated with the simple system made of movable rigid leaflet and an unsteady viscous fluid flow. Such two-dimensional model, even if it appears extremely simple, is able to produce fairly rich flow features which deserve careful analytical and accurate numerical solution. Thus, we have performed a significant number of numerical experiments with the objective to uncover the role of the structure inertia in the overall behavior of the fluid-leaflet system, under the different flow recurrences. For that purpose, we have constructed a strong-coupling code and resolved the fluid and structure dynamics simultaneously, paying particular care of solution accuracy around the moving boundary. The complex problem of large fluid deformation in response to the rapid structure movements has been resolved by the time-dependent conformal mapping, exclusively developed for this specific physical arrangement. The numerical findings, even if theoretical in nature, allowed for the classification and characterization of body’s and fluid dynamics in functionality of different structure inertia and Strouhal numbers, which have been used as free parameters in all numerical experiments. The study is completed by a brief analysis of the more realistic system of actual prosthetic heart valves. Besides many problems that follow the performance of mechanical heart valve prosthesis, the complications related to the complex blood-leaflet interaction are a key factor. The intraventricular flow is characterized by large vortical structures, without significant turbulence, in a smooth circulatory pattern that, in presence of pathological conditions or mechanical devices, could be disturbed. Thus, among the criteria for the assessment of mitral valve functionality and mechanical valve design are the proper vortical features inside the left ventricle. Until nowadays the standard mechanical valves, designed originally for the aortic replacement and without exceptions symmetrical, have never satisfied the regularity of natural vortical dynamics. Thus, we have been motivated to investigate the flow features downstream of asymmetrical prototypes, exclusively designed for the mitral replacement with attempt to better mimic the natural intraventricular flow. Experimental outcomes allowed for preliminary conclusions that the break of symmetry in the novel prosthesis creates the asymmetrical vortical flow in the left ventricle, which is more similar to the natural one, although the concept introduced by this prototype has to undergo deeper testing and careful improvements before querying in the real hearts.
XXIII Ciclo
1982

2011/2012
Evaporation and condensation of thin liquid films on solid surfaces are common elements of industrial processes. In many cases they have a significant impact on the physics of the studied case. At the same time, experimental studies can prove to be troublesome, mostly because of the amount of possible setups, complex geometries of interest, numerous materials being used and cost. For that reason it is reasonable to study this phenomena using numerical methods. Having the advantage in speed and cost of performance, computational studies become a valuable tool. For evaporation and condensation process, one has to deal with buoyancy driven fluid flows, conjugated heat transfer between gaseous and solid phases, film thickness modeling, vapor phase behavior, and phase transition of the thin fluid film into vapor phase. The strong conjunction and mutual interaction of mentioned effects is the main focus of presented work. The gas phase behavior is being calculated using incompressible Navier-Stokes equations under Boussinesq approximation. The solutions of the partial differential equations are obtained with numerical methods using Eulerian finite volume discretization (Kundu and Cohen [2002]). Time advancement is being treated with second order implicit discretization. For cases with high Reynolds number, large eddy simulation (LES) techniques are used. Due to the complexity of the geometries of interest a dynamic computation of the Smagorinsky constant is preferred, applying the lagrangian dynamic model proposed by Meneveau et al. [1996]. The liquid film present on the surface of the solids is modeled following Petronio[2010]. Since the film is thin, it is assumed that it can be represented only by its thickness. This also leads to assumption that the heat transfer through the film is instantaneous. The vapor is represented by concentration of this phase in the volume of gas. The concentration is transported by convection and diffusion. The phenomena of evaporation and condensation of the thin films are driven by the presence of concentration gradients next to the surfaces. Phase transition of vapor to fluid, or other way around, acts on the energy balance, id. est latent heat is released into the gas when condensation occurs or the solid is cooled during evaporation. The heat transport is modeled in both solid and fluid domains. The case is split into separate regions with different material properties. These regions are solved one by one in a serial way using numerical techniques consistent with domain decomposition methods described by Quarteroni and Valli [1999]. The energy transport among the regions is performed by applying a heat coupling boundary conditions. The main focus of this work is to provide a reliable model for simulation system with complex physics involving fluid motion, heat transport in multi region domains (fluid-solid), vapor transport, thin film evolution and evaporation and condensation effects on energy balance. Proposed model is validated on simple geometries and later applied to problem of evaporation in vertical channel flow. The reference to the channel case is work of Laaroussi et at. [2009]. Presented study aims in providing comprehensive insights into physical effects that appear when the solid wall is being directly modeled and when latent heat transformations are taken into account. The final test is performed on a vertical channel with forced turbulent flow, directly modeled solid walls and evaporation or condensation happening on the boundary. Having the model working within such complex frame allows for its future usage in elaborate industrial applications.
XXV Ciclo
1985

2009/2010
L'obiettivo principale del progetto di ricerca sviluppato nella presente tesi di dottorato è quello di comprendere meglio le problematiche riguardanti le prestazioni di asciugatura della lavastoviglie, con speciale riferimento alla modellazione dei fenomeni di evaporazione e condensazione che avvengono nella vasca. Tipicamente l'acqua per il risciacquo finale viene portata ad una temperatura di 70°C che scalda le stoviglie permettendo a queste di immagazzinare energia termica. La lavastoviglie si raffredda dall'esterno, cosicché la vasca risulta essere più fredda dei piatti posti all'interno. In questo contesto l'acqua può evaporare dalle superfici delle stoviglie e condensare sulle pareti della vasca stessa. Il sistema fisico può essere descritto come un flusso in presenza di cambiamenti di fase. Tali tipologia di flussi ha un ruolo cruciale in molti processi naturali e tecnologici, in particolare in quelli in cui si hanno asciugatura o formazione di condensa sulle superfici solide. Tuttavia, pur essendo così comuni nelle applicazioni ingegneristiche, la loro comprensione è lontana dall'essere completa. Il complesso problema fisico può essere suddiviso in tre sotto-problemi: la trasmissione del calore tra il corpo bagnato ed il liquido sulla sua superficie; il trasferimento di calore e massa tra la fase liquida e quella gassosa; il flusso della fase gassosa che risulta essere molto influenzato dalle forze di galleggiamento dovute alle variazioni di densità causate dalla diffusione di temperatura e concentrazione di vapor acqueo. Dallo studio della letteratura risulta che tale problema non sia stato ancora investigato completamente. In particolare non è mai stato proposto un modello adatto a scopi ingegneristici, cioè per problemi di larga scala con geometrie complesse, che consideri l'evoluzione del film liquido durante processi di asciugatura. Questo progetto vuole contribuire allo sviluppo della ricerca in questo settore. Il modello matematico del flusso d'aria in presenza di evaporazione e condensazione è stato implementato numericamente nell'ambiente open-source OpenFoam. Il modello consiste nella formulazione delle equazioni di Navier-Stokes per flussi incomprimibili più le equazioni del trasporto per la temperatura e la concentrazione di vapore. Entrambi gli scalari sono considerati attivi e le variazioni di densità sono state incorporate sotto l'approssimazione di Boussinesq. Si assume inoltre l'approssimazione a film sottile, per cui si è inteso che film liquidi, gocce ed, in genere, le zone bagnate di un solido possano essere considerate come un film liquido continuo. Tale film sottile è stato interpretato come una condizione al contorno per il flusso d'aria, prescrivendo una condizione di Dirichlet per la temperatura e per il vapore. Quest'ultimo all'interfaccia del liquido è considerato in condizione di saturazione. Il calcolo della velocità di evaporazione all'interfaccia, imposta anche come condizione al contorno per il campo di velocità, ha permesso la quantificazione del processo di evaporazione/condensazione consentendo il calcolo della massa d'acqua evaporata/condensata. Il modello numerico è stato validato con i dati di letteratura per poi essere applicato nello studio del flusso su un cilindro bagnato, tra due piani paralleli. In questo lavoro è stato evidenziato l'effetto sul flusso di evaporazione attorno al cilindro delle condizioni alle pareti, considerate come bagnate o asciutte ed adiabatiche. Inoltre è stato valutato anche l'effetto della distanza del cilindro stesso dalle pareti. Successivamente il modello è stato applicato ad una geometria 2D della lavastoviglie. I risultati mostrano che il flusso evolve secondo uno schema preciso: le forze di galleggiamento danno luogo ad un moto convettivo che si alza dalle stoviglie più calde ed umide, che poi scende lungo alle pareti più fredde e meno umide. Un'ulteriore analisi è stata fatta simulando il processo fino all'asciugatura completa di un film uniformemente distribuito su tutte le stoviglie. Negli stadi intermedi del processo è stato osservato che, attorno alle porzioni di stoviglie già asciutte, il galleggiamento risulta essere ridotto e la velocità dell'aria minore, per il mancato rilascio di vapor acqueo. Un passo ulteriore verso la modellazione della lavastoviglie è stato condotto considerando una geometria 3D semplificata per testare il modello e verificare le caratteristiche richieste alla griglia computazionale. Anche per questa configurazione è stato osservato l'instaurarsi del moto convettivo e l'effetto dell'asciugatura sul flusso. Infine si è iniziato a studiare il caso della lavastoviglie 3D. La simulazione è potuta durare pochi secondi fisici, nei quali hanno iniziato a svilupparsi sopra le stoviglie i caratteristici plume. Successivamente delle instabilità numeriche hanno dato luogo a valori di pressione non fisici determinando l'interruzione del programma. Tale comportamento è stato spiegato dalla mancata dissipazione turbolenta nel flusso. L'attivazione del modello LES di Smagorinsky con l'analogia di Reynolds per la determinazione delle diffusività turbolente dei due scalari ha dato luogo ad una soluzione numericamente stabile. Tuttavia la eccessiva viscosità di sotto-griglia ha sovrastimato la diffusione degli scalari, inficiando l'accuratezza della simulazione. Per includere nel modello l'accoppiamento termico che caratterizza il processo di asciugatura è stata scelta la tecnica di decomposizione di domini detta di Dirichlet-Neumann in quanto è risultata essere efficace e semplice da implementare. Essa impone la continuità della temperatura e il bilancio dei flussi di calore attraverso le interfacce. Inoltre è stato proposto un modello opportuno per la distribuzione della temperatura nel film liquido, per riprodurre nella maniera corretta il trasferimento di calore attraverso il film stesso. Ciascuna di queste due parti è stata implementata e testata individualmente cosicché la loro inclusione nel modello potrà avvenire in un successivo sviluppo della presente ricerca. Inoltre è in fase di sviluppo l'implementazione del modello di sotto-griglia LES dinamico lagrangiano che permetterà di superare i limiti riscontrati nel utilizzo del modello di Smagorinsky.
XXIII Ciclo
1980

2010/2011
La presente tesi concerne la modellizzazione e simulazione numerica, attraverso l'equazione di Gross-Pitaevskii (chiamata anche equazione di Schroedinger non lineare), della dinamica dei vortici quantistici nell'elio superfluido e in particolare del fenomeno della riconnessione. La riconnessione si verifica qualora due vortici approssimativamente antiparalleli, si intersecano e si scambiano le estremità. Questo fenomeno è stato osservato sperimentalmente e risulta essere una caratteristica essenziale della turbolenza quantistica.
XXIII Ciclo
1983

2011/2012
In questo lavoro l'analisi dello spazio delle fasi viene applicata a tre sistemi di equazioni differenziali, due riguardanti i fluidi geofisici e uno riguardante le onde termoelastiche. Il primo risultato mostra l'esistenza e l'unicità di soluzioni "mild" per l'equazione di Navier-Stokes-Coriolis, nel caso i dati iniziali siano piccoli nella norma dello spazio ibrido $B_{\dot{H}^{\frac{1}{2}},\dot{B}^{\frac{3}{p}-1}_{p,\infty}}$, con $3XXIV Ciclo
1984

2013/2014
Renewable energy resources, such as wind, are available worldwide. Locating areas with high and continual wind sources are crucial in pre-planning of wind farms. Vast offshore areas are characterized with higher and more reliable wind resources in comparison with continental areas. However, offshore wind energy production is in a quite preliminary phase. Elaborating the potential productivity of wind farms over such areas is challenging due to sparse in situ observations. Mediterranean basin is not an exception. The overall aim of this thesis is to perform analysis in model efficiency in estimation of wind energy from regional to local scale. First, we are proposing numerical simulations of near-surface wind fields from regional climate models (RCMs) in order to obtain and fill the gaps in observations over the Mediterranean basin. Four simulations produced with two regional climate models are examined. Remote sensing observations (QuikSCAT satellite) are used to assess the skill of the simulated fields. A technique in estimation the potential energy from the wind fields over the region is introduced locating the three potentially interesting sub-regions for wind farms. Then, we use local-scale model (large-eddy simulation) with implemented parameterization of wind turbine in order to simulate real case flow in theoretical wind farm. Information reported with regional climate model would be used to create inflow conditions for the selected sub-region of the Mediterranean Sea for simulating theoretical offshore wind farm. Finally, we would compare the estimation of wind power potential obtained by regional climate model and power production of theoretical wind farm obtained with large-eddy simulations for chosen sub-region. Within this multi-scale approach, we would present different numerical computational efficiency in application of wind energy and justification in usage of both regional and local scale models. The novelty of this multi-model methodological approach could be considered in offering significant information for wind industry.
XXVII Ciclo
1983

2012/2013
Turbulent Rayleigh-Bénard convection (RBC) occurs when a shallow layer of fluid is heated from below. It is a challenging subject in non-linear physics, with many important applications in natural and engineering systems. Because of the complexity of the governing equations, analytical progress in understanding convection has been slow, and laboratory experiments and numerical simulations have assumed increased importance. In regard to numerical work, Large-Eddy Simulation (LES) techniques have proved to be reliable and powerful tool to understand the physics since it provides better coverage for measurements, that are not as easily obtained in physical experiments or the other numerical approaches. This thesis addresses different aspects of Rayleigh-Bénard convection in fully developed turbulent regime through Large Eddy Simulation (LES) to shed light on some important aspect of the geometrical shape of the convection cell. The layout of the thesis is as follows: In Chapter 1, we first introduce Rayleigh-Bénard convection and the equations and parameters that govern it. This is followed by a discussion on different types of boundary conditions used in numerical and theoretical studies of RBC. Subsequently we present various convection states that are observed analytically and experimentally in RBC as a function of Ra and Ʈ. To this end we present a brief survey of the analytical, experimental and numerical works on confined thermal convection. We introduce different regimes and related scaling according to Grossman and Lohse theory. We also present the experimental and numerical results related to the Large Scale Circulation (LSC) within different geometries. In Chapter 2, we present the details of the numerical methods used to solve the governing non-linear equations . In the second part, we provide the details of the solver and the algorithm used to solve the RBC dynamical equations in a Cartesian geometry together with boundary conditions. In Chapter 3, we demonstrate that our numerical method and solver give results consistent with earlier numerical results. Results from the Direct Numerical Simulations (DNS) and Large Eddy Simulation (LES) with constant and dynamic subgrid scale Prandtl number (P_sgs) are presented and compared. We observe close agreement with Lagrangian dynamic approaches. In the first part of Chapter 4 we analyse the local fluctuations of turbulent Rayleigh-Bénard convection in a cubic confinement with aspect ratio one for Prandtl number Pr = 0.7 and Rayleigh numbers (Ra) up to 10^9 by means of LES methodology on coarse grids. Our results reveal that the scaling of the root-mean-square density and velocity fluctuations measured in the cell center are in excellent agreement with the unexpected scaling measured in the laboratory experiments of Daya and Ecke (2001) in their square cross-section cell. Moreover we find that the time-averaged spatial distributions of density fluctuations show a fixed inhomogeneity that maintains its own structure when the flow switches from one diagonal to the other. The largest level of rms density fluctuations corresponds to the diagonal opposite that of the Large Scale Circulation (LSC) where we observed strong counter-rotating vortex structures. In the second part we extended our simulations and Ra up to 1011, in order to identify the time periods in which the orientation of LSC is constant. Surprisingly we find that the LSC switches stochastically from one diagonal to the other. In Chapter 5, we study the effect of 3D-roughness on scaling of Nu(Ra) and consequently on the fluctuations of density. Moreover we present the effect of roughness shape when the tip has a wide angle and the other one is smooth. We study two types of elements, one of which is a pyramid and the other is a sinusoidal function spread over the bottom (heated) and top (cooled) plates, in a cubic confinement. However preliminary results suggest that the effect of roughness appears evident at high Ra numbers when the thermal boundary layer is thin enough to shape around the obstacles.
XXVI Ciclo
1983

2013/2014
Mostriamo un modello allo stato dell’arte, che considera i principali processi fisici che governano il greggio in mare nelle prime ore dopo il rilascio, (Zanier, et al., 2014). Le particelle e i tar sono trattati come particelle lagrangiane, ognuna con la propria densità e il proprio diametro; consideriamo le forze principali che agiscono su di esse ossia: galleggiamento, trascinamento e la forza di Coriolis. Il greggio in forma di film sottile è modellato tramite le equazioni proposte da Nihoul (Nihoul 1983/84). Il modello originale di Nihoul considera le forze principali (ossia gravità, stress indotto da vento e correnti marine) che agiscono sulla macchia e governano il suo trasporto e diffusione, sulla superficie del mare, nelle prime 24 ore dopo il rilascio. Il nostro miglioramento al modello consiste nell’introduzione della forza di Coriolis evitando di utilizzare formulazioni empiriche (Zanier, et al., 2015). Infine i principali processi di weathering che agiscono sulla macchia nelle prime 12-24 ore dopo il rilascio (ossia emulsificazione ed evaporazione) sono considerate in accordo con i modelli presenti in letteratura (Mackay, Peterson, et al., 1980 e Mackay, Buist, et al., 1980, rispettivamente). Per preservare un’accuratezza del secondo ordine del metodo numerico, i termini convettivi, nel modello Euleriano, sono discretizzati usando SMART uno schema numerico upwind del terzo ordine (Gaskell and Lau 1988). Il modello è validato con dei casi test standard. Le correnti marine sono risolte con il modello LES-COAST (IEFLUIDS Università di Trieste), un modello numerico ad alta definizione, adatto per simulare flussi in aree costiere e portuali. Il modello LES-COAST risolve la forma filtrata delle equazioni di Navier-Stokes tridimensionali e non-idrostatiche, assumendo che valga l’approssimazione di Boussinesq; e l’equazione di trasporto degli scalari, salinità e temperatura. Il modello usa l’approccio della large eddy simulation per parametrizzare la turbolenza, le variabili sono filtrate con una funzione filtro, rappresentante la grandezza delle celle. I flussi di sottogriglia (SGS), che appaiono dopo l’operazione di filtraggio delle equazioni, sono parametrizzati con un modello di Smagorinsky anisotropo con due eddy viscosity, per adattare il modello a simulare flussi costieri dove le lunghezze scala orizzontali sono molto più grandi di quelle verticali (Roman et al., 2010 ). Le diffusività di sotto griglia della temperatura e salinità, cioè i numeri di Prandtl e Schmidt, sono imposti come $Pr_{sgs}=Sc_{sgs}=0.8$, assumendo che l’analogia di Reynolds sia valida per entrambi gli scalari. La complessità geometrica che caratterizza le aree costiere, è trattata con una combinazione di griglie curvilinee e il metodo dei contorni immersi (IBM) (Roman, Napoli, et al., 2009). L’azione del vento sulla superficie libera del mare è imposta tramite una formula proposta da Wu (Wu, 1982), nella quale lo stress del vento sul mare è calcolato dalla velocità del vento a 10 m sopra il livello del mare. Allo stress aggiungiamo una varianza del 20% per agevolare la generazione di turbolenza e per tener conto che l’azione del vento non è costante nel tempo e nello spazio. Inoltre vicino agli ostacoli, come moli, navi e frangiflutti, lo stress del vento è ridotto linearmente, per considerare la riduzione del vento che si ha nelle zone di ricircolo. Sui contorni aperti le velocità e le quantità scalari sono ottenute innestando il modello LES-COAST con modelli di larga scala (Petronio, et al., 2013) oppure sono impostati secondo dati rilevati. Vicino ai bordi solidi le velocità sono modellate tramite funzioni parete (Roman, Armenio, et al., 2009). Il modello di rilascio di petrolio e il modello idrodinamico sono stati applicati assieme per simulare degli ipotetici scenari di trasporto e diffusione del greggio in mare nel porto di Barcellona (Mar Mediterraneo Nord-Ovest, Spagna, Galea, et al. 2014) e nella baia di Panzano (Mar Adriatico, Nord, Italia).
We present a novel, state of the art model, which accounts for the relevant short-term physical processes governing oil spill at sea, (Zanier, et al., 2014). Particles and tars are modelled as Lagrangian phase having its own density and diameter; taking into account the main forces acting on them, namely: buoyancy, drag and Coriolis forces. Oil transported in form of thin-film is treated by means of an improved Nihoul’s model (Nihoul 1983/84). The latter considers the main forces (gravity, wind and sea currents stresses) governing oil slick spreading and transport in the first hours after spilling, up to 24h for large spill. Our main improvement to the classical model consists in the introduction of Coriolis effect, avoiding using empirical formulations (Zanier, et al., 2015). Finally the relevant short-term (12-24 hours) weathering processes (mainly emulsification and evaporation) are taken into account through established literature models (Mackay, Peterson, et al., 1980 and Mackay, Buist, et al., 1980, respectively). To preserve second-order accuracy of the overall numerical method, convective terms, in the Eulerian model, are discretized using SMART a third order accurate upwind numerical scheme (Gaskell and Lau 1988). We validate the model on standard test cases. The underground hydrodynamics is resolved using LES-COAST (IEFLUIDS University of Trieste), a high definition numerical model suited for coastal or harbour areas. LES-COAST model solves the filtered form of three dimensional, non-hydrostatic Navier-Stokes equations under Boussinesq approximation and the transport equation for salinity and temperature. It makes use of Large Eddy Simulation approach to parametrize turbulence, the variables are filtered by way of a top-hat filter function represented by the size of the cells. The subgrid-scale fluxes (SGS), which appear after filtering operations, are parametrized by a two-eddy viscosity anisotropic Smagorinsky model, to better adapt to coastal flow in which horizontal length scale is larger than vertical one (Roman et al., 2010). The subgrid-scale eddy diffusivities of temperature and salinity, Prandtl and Schmidt numbers, are set $Pr_{sgs}=Sc_{sgs}=0.8$, by assuming that Reynolds analogy holds also for both scalars. Complex geometry that characterizes coastal flow is treated by a combination of curvilinear grid and Immersed Boundary Method (IBM) (Roman, Napoli, et al., 2009). Wind action on the free surface is taken into account by means of the formula proposed by Wu (Wu, 1982), in which the wind stress on the sea surface is computed from the wind velocities at 10 m above the surface. A 20% of variance is added to the stress to ease the generation of turbulence and to take into account of wind stress variations in time and space. Moreover near obstacles such as docks, ships and breakwaters, the wind stress is linearly reduced considering the relevant reduction of stress in recirculation regions. On the open boundaries the velocities and scalars quantities are obtained by nesting LES-COAST within Large Circulation Models (Petronio, et al., 2013) or are imposed from in-situ measurements. Near the wall velocities are modelled using wall functions (Roman, Armenio, et al., 2009). We apply the coupled oil spill model and hydrodynamical one to simulate hypothetical oil spill events in real case scenarios in Barcelona harbour (North-west Mediterranean Sea, Spain, Galea, et al. 2014) and in Panzano bay (North Adriatic Sea, Italy).
XXVII Ciclo
1986

2012/2013
The subject of modelling flow near wall is still open in turbulent wall bounded flows, since there is no wall layer model which works perfectly. Most of the present models work well in attached flows, specially for very simple geometries like plane channel flows. Weakness of the models appears in complex geometries, and many of them do not capture flow separation accurately in detached flows, specially when the slope of wall changes gradually. In many engineering applications, we deal with complex geometries. A possible way to simulate flows influenced by complex geometry using a structured grid, is to consider the geometry as immersed boundary for the simulation. Current wall layer models for the immersed boundaries are more complex and less accurate than the body-fitted cases (cases without immersed boundaries). In this project the accuracy of wall layer model in high Reynolds number flows is targeted, using LES for attached flows as well as detached flows (flows with separation). In addition to the body fitted cases, wall layer model in the presence of immersed boundaries which is treated totally different also regarded. A single solver LES-COAST (IE-Fluids, University of Trieste) is used for the flow simulations, and the aim is to improve wall layer model in the cases with uniform coarse grid. This is in fact novelty of the thesis to introduce a wall layer model applied on the first off-wall computational node of a uniform coarse grid, and merely use LES on the whole domain. This work for the immersed boundaries is in continuation of the methodology proposed by Roman et al. (2009) in which velocities at the cells next to immersed boundaries are reconstructed analytically from law of the wall. In body-fitted cases, since smaller Smagorinsky constant is required close to the walls than the other points, wall layer model in dynamic Smagorinsky sub-grid scale model using dynamic k (instead of Von Karman constant) is applied to optimize wall function in separated flows. In the presence of immersed boundaries, the present wall layer model is calibrated, and then improved in attached and also detached flows with two different approaches. The results are also compared to experiment and resolved LES. Consequently the optimized wall layer models show an acceptable accuracy, and are more reliable. In the last part of this thesis, LES is applied to model the wave and wind driven sea water circulation in Kaneohe bay, which is a bay with a massive coral reef. This is the first time that LES-COAST is applied on a reef-lagoon system which is very challenging since the bathymetry changes very steeply. For example the water depth differs from less than 1 meter over the reef to more than 10 meters in vicinity of the reef, in lagoon. Since a static grid is implemented, the effect of wave is imposed as the velocity of current over the reef, which is used on the boundary of our computational domain. Two eddies Smagorinsky SGS model is used for this simulation.
XXVI Ciclo
1983

2010/2011
The surface transport in the Northeastern Adriatic Sea has been investigated by evaluating, for the first time, the Finite-Size Lyapunov Exponent (FSLE) on the current field detected by the High Frequency (HF) coastal radar network active in the period August 2007 - August 2008. A similar analysis has been carried out on the MITgcm (Massachusetts Institute of Technology general circulation model) current field in order to have a perspective on the various results given by the application of the same FSLE evaluation procedure to different velocity fields. This work includes also the development, testing and calibration of the FSLE evaluation algorithm for the specific Adriatic area considered. The aim of this thesis is to study the surface dynamics of the Northeastern Adriatic current field associated with the typical wind regimes, namely Bora and Sirocco. The strongest and most persistent wind events, coinciding with the radar network activity, have been identified in the first instance from ALADIN (Aire Limitée Adaptation dynamique Développement InterNational) model meteorological data and then confirmed by the in situ meteo-mareographic time series in Trieste and Venice. For comparison purposes, the transport structures developing throughout the longest calm wind periods have also been investigated. In this thesis, the combination of the FSLE resulting from model and radar current fields contributes to interpret the surface transport dynamics in the studied area. In fact, it allows: i) to determine the strongest dynamical features, ii) to extend the transport information beyond the radar coverage and iii) to bring to the light the possible ambiguity of some structures originated from model currents. The FSLE analysis evidenced along the northern Adriatic margin an attractive transport structure with a filament-shaped conformation. The dynamics of this structure is driven by the water exchange between the Gulf of Trieste (GoT) and the North Adriatic Sea. The spatial location of this filament and the advective direction associated with it, vary according to the wind patterns. In fact, in calm wind periods this attractive filament is found right in front of the gulf entrance indicating the pattern of the GoT outflow (east-west direction). During Bora episodes this attractive filament is advected westward and it elongates following the northern Adriatic coast. Differently, Sirocco wind pushes this filament closer to the northeastern coastline reversing the transport direction along it (from west to east). Previous studies evidenced that the current signal in the southern part of the domain is less correlated to the wind pattern with respect to what observed in the northern area of the current field. Therefore also from the FSLE analysis no recurrent transport dynamics is observed in the southern area except for the Bora cases, when a repulsive structure originates from the Istrian coast and it looses strength while it is advected northwestward.
Il trasporto superficiale nella parte nord-orientale del Mar Adriatico è stato studiato calcolando, per la prima volta, i Finite-Size Lyapunov Exponent (FSLE) sui campi di corrente misurati dalla rete di radar ad alta frequenza attiva durante il periodo Agosto 2007 - Agosto 2008. Uno studio analogo è stato effettuato sui campi di corrente del modello MITgcm (Massachusetts Institute of Technology general circulation model), per avere una panoramica sulla varietà dei risultati prodotti dalla valutazione degli FSLE su diversi campi di velocità. Il lavoro include anche lo sviluppo, i test e la calibrazione dell'algoritmo di calcolo degli FSLE per la specifica area dell'Adriatico considerata. L'obiettivo è analizzare la dinamica superficiale del campo di corrente associata ai regimi di vento tipici dell'area nord adriatica, ovvero Bora e Scirocco. Gli episodi di vento più intensi e persistenti, verificatisi durante il periodo di attività dei radar, sono stati identificati in primo luogo dai dati meteorologici prodotti dal modello ALADIN (Aire Limitée Adaptation dynamique Développement InterNational) e successivamente sono stati confermati dalle serie temporali dei dati meteo-mareografici delle stazioni di Trieste e Venezia. A scopo comparativo, sono state analizzate anche le strutture di trasporto sviluppatesi durante i prolungati periodi di calma di vento. In questa tesi, la combinazione degli FSLE risultanti dai campi di corrente dei radar e del modello contribuisce ad una più completa interpretazione della dinamica del trasporto superficiale nell'area di interesse. I risultati ottenuti permettono: i) di determinare le strutture dinamiche più intense, ii) di estendere l'informazione sul trasporto al di fuori della copertura della rete dei radar, nonché iii) di portare alla luce eventuali ambiguità di alcune strutture originatesi dalle correnti del modello. L'analisi degli FSLE ha evidenziato una struttura di trasporto attrattiva che assume l'aspetto di un filamento lungo il margine settentrionale del Mar Adriatico. La dinamica di questa struttura è regolata dagli scambi di acqua tra il Golfo di Trieste ed il Nord Adriatico. La posizione del filamento e la direzione avvettiva associata ad esso dipendono dal regime di vento. Infatti, nei periodi di calma di vento il filamento attrattivo si trova di fronte all'entrata del golfo, ad indicare la direzione del flusso d'acqua in uscita dal golfo stesso (da est ad ovest). Durante gli episodi di Bora il filamento attrattivo è trasportato verso ovest e si allunga seguendo la costa nord adriatica. Invece, lo Scirocco spinge questo filamento verso la costa nord-orientale invertendo la direzione del trasporto associata ad esso (da ovest ad est). Precedenti lavori hanno evidenziato che le correnti nella parte meridionale del dominio sono meno correlate al segnale di vento, rispetto a quanto osservato nell'area settentrionale del campo delle correnti. Quindi anche dall'analisi degli FSLE non si osserva una dinamica di trasporto ricorrente ad eccezione dei casi di Bora, durante i quali si sviluppa una struttura repulsiva lungo la costa Istriana, la cui intensità diminuisce mentre viene trasportata verso nordovest.
XXIV Ciclo
1984

2012/2013
We investigate by different techniques, the solvability of a class of capillarity-type problems, in a bounded N-dimensional domain. Since our approach is variational, the natural context where this problem has to be settled is the space of bounded variation functions. Solutions of our equation are defined as subcritical points of the associated action functional.
We first introduce a lower and upper solution method in the space of bounded variation functions. We prove the existence of solutions in the case where the lower solution is smaller than the upper solution. A solution, bracketed by the given lower and upper solutions, is obtained as a local minimizer of the associated functional without any assumption on the boundedness of the right-hand side of the equation. In this context we also prove order stability results for the minimum and the maximum solution lying between the given lower and upper solutions. Next we develop an asymmetric version of the Poincaré inequality in the space of bounded variation functions. Several properties of the curve C are then derived and basically relying on these results, we discuss the solvability of the capillarity-type problem, assuming a suitable control on the interaction of the supremum and the infimum of the function at the right-hand side with the curve C. Non-existence and multiplicity results are investigated as well. The one-dimensional case, which sometimes presents a different behaviour, is also discussed. In particular, we provide an existence result which recovers the case of non-ordered lower and upper solutions.
XXV Ciclo
1985

2012/2013
The Regional Climate Model RegCM4 (Giorgi et al., 2012) treats nonconvective clouds and precipitation following the Subgrid Explicit SUBEX param- eterization (Pal et al., 2000). This scheme includes a simple representation for the formation of raindrops and solves diagnostically the precipitation: rain forms when the cloud water content exceeds the autoconversion threshold, that is an increasing function of the temperature and assumes different values over the land and over the ocean to account for the difference in number of the cloud condensation nuclei over continental and oceanic regions. The SUBEX scheme does not account for the presence of clouds ice, and the fraction of ice is diagnosed as a function of temperature in the radiation scheme. Due to the increasing emphasis on cloud representations in the climate community and the forthcoming increasing resolution due to the inclusion, in the close future, of a non-hydrostatic compressible core, the treatment of the ice microphysics and a prognostic representation of the precipitation is required in RegCM4. This thesis presents the new parameterization for stratiform cloud microphysics and precipitation implemented in RegCM4. The approach of the new parameterization is based on an implicit numerical framework recently developed and implemented into the ECMWF operational forecasting model (Tiedtke, 1993). The new parameterization solves 5 prognostic equations for the water vapour, the liquid water, the rain, the ice and the snow mixing ratios. It allows a proper treatment of mixed-phase clouds and a more physically realistic representation of the precipitation as it is no more an instantaneous response to the microphysical processes occurring in clouds and is subjected to the horizontal advection. A first discussion of the results contains an evaluation of the vertical distributions of the main microphysical quantities, such as the liquid and ice water mixing ratios and the relative fractions. It also presents a series of sensitivity tests to understand how the moisture and radiation quantities respond to the variation of the microphysical parameters used in the scheme, such as the fall speeds of the falling categories, the autoconversion scheme and the evaporation coefficient. Cloud properties are afterwards evaluated through the implementation for RegCM4 of the new cloud evaluation COSP tool (Bodas-Salcedo et al., 2011), developed by the Cloud Feedback Model In- tercomparison Project (CFMIP), that facilitates the comparison of simulated clouds with observations from passive and active remote sensing by diagnosing from model outputs the quantities that would be observed from satellites if they were flying above an atmosphere similar to that predicted by the model. Different hypothesis are presented to explain the reasons for RegCM4 biases in representing different types of clouds over the tropical band and new prospectives for the future investigations designed to answer to the open questions are outlined.
XXVI Ciclo
1983

2012/2013
A novel high-resolution, eddy-resolving numerical model (LES-COAST) is used to investigate currents, mixing and water renewal in Barcelona harbour and Taranto bay. These environmental sites are of particular importance due to the interplay between touristic and commercial activities, requiring detailed and high-definition studies of water quality within the harbour. We use Large Eddy Simulation (LES) which directly resolves the anisotropic and energetic large scales of motion and parametrizes the small, dissipative, ones. Small-scale turbulence is modelled by the Anisotropic Smagorinsky Model (ASM) which is employed in presence of large cell anisotropy. The complexity of the harbour is modelled using a combination of curvilinear, structured, non-staggered grid and the Immersed Boundary Method (IBM). Both computation grids and harbour structures are purposely constructed for these applications by appropriate programs. Boundary conditions for wind forcing at the free surface and currents at the inlets of the port are obtained from in-situ measurements (for the case of Barcelona harbour) or by nesting this numerical model into a coastal model (Taranto bay). In this dissertation thesis an important modification to LES-COAST is implemented and is proposed as a prototype scheme, namely the possibility to consider the effect of surface waves in coastal semi-closed areas. Particularly, a linear formulation of the free surface boundary condition is considered, which would be able to reproduce the presence of seiches and tides on the dynamics of the area under investigation. The methodology is validated against analytical solution for a stationary oscillating surface wave in a simple computational grid.In both harbours considered, first- and second-order statistics, such as the mean velocity field, turbulent kinetic energy, and horizontal and vertical eddy viscosities are calculated and their spatial distribution is assessed. Water residence time is also considered for the two coastal semi-closed areas examined. Finally, the LES solution is validated against available field data.The study shows the presence of sub-surface elongated rolling structures (with a time scale of a few hours), contributing to the vertical water mixing. The time-averaged velocity field reveals intense upwelling and downwelling zones along the walls of the harbours. The analysis of second-order statistics in these harbours shows strong inhomogeneity of turbulent kinetic energy and horizontal and vertical eddy viscosities in the horizontal plane, with larger values in the regions characterized by stronger currents. The water renewal within the port is quantified for particular sub-domain regions, showing that the complexity of the harbour is such that certain inner basins of Barcelona harbour have a water renewal of over five days, including its yacht marina area, and over seven days for Taranto bay. For the Barcelona simulation, the LES solution compares favourably with available current-meter data; it is also compared with a RANS solution obtained in literature for the same site under the same forcing conditions, the comparison demonstrating a large sensitivity of properties to model resolution and frictional parametrization.
XXVI Ciclo
1986

2013/2014
Questa tesi è dedicata allo studio di alcuni modelli differenziali che nascono nell'ambito della fluidodinamica o della relatività generale e che coinvolgono gli operatori di curvatura media nello spazio $N$-dimensionale euclideo o di Minkowski. Entrambi sono operatori ellittici quasi-lineari che non soddisfano la proprietà di uniforme ellitticità, essendo l'operatore di curvatura media euclidea degenere, mentre quello di curvatura media nello spazio di Minkowski singolare. Il lavoro è suddiviso in tre parti. La prima riguarda lo studio delle soluzioni periodiche dell'equazione di curvatura prescritta unidimensionale nello spazio euclideo, equazione che modellizza fenomeni di tipo capillarità. In accordo con la struttura dell'operatore di curvatura e imponendo un opportuno comportamento in 0, o all'infinito, della curvatura prescritta, si dimostra l'esistenza di infinite soluzioni subarmoniche classiche arbitrariamente piccole aventi opportune proprietà nodali, oppure di infinite soluzioni subarmoniche a variazione limitata con oscillazioni arbitrariamente grandi. La tecnica per la ricerca delle soluzioni classiche è topologica e si basa sull'uso del numero di rotazione e su una generalizzazione del teorema di Poincaré-Birkhoff; d'altro lato l'approccio per lo studio delle soluzioni non classiche poggia sulla teoria dei punti critici per funzionali non lisci, in particolare su un lemma di passo di montagna nello spazio delle funzioni a variazione limitata. La seconda parte della tesi è dedicata allo studio del problema di Dirichlet omogeneo associato a un'equazione della curvatura media prescritta anisotropa nello spazio euclideo, il quale fornisce un modello di descrizione della geometria della cornea umana. Il problema è ambientato in un dominio regolare in $\mathbb{R}^N$ con frontiera lipschitziana. Il capitolo è suddiviso a sua volta in tre sezioni, che sono rispettivamente focalizzate sui casi unidimensionale, radiale e $N$-dimensionale. Nel caso unidimensionale e nel caso radiale in una palla, si dimostrano l'esistenza e l'unicità di una soluzione classica, che presenta alcune proprietà qualitative aggiuntive. Le tecniche usate in questo contesto sono di natura topologica. Infine, nel caso $N$-dimensionale in un dominio generale, si provano l'esistenza, l'unicità e la regolarità di una soluzione di tipo forte del problema. In relazione ai possibili fenomeni di scoppio del gradiente, l'approccio è variazionale nello spazio delle funzioni a variazione limitata. Si enunciano e si dimostrano prima di tutto alcuni risultati preliminari riguardo al comportamento del funzionale associato al problema; tra questi, si sottolinea l'importanza di una proprietà di approssimazione. Successivamente si provano l'esistenza e l'unicità del minimizzante globale del funzionale, che è regolare all'interno ma non necessariamente sulla frontiera, e soddisfa il problema secondo un'opportuna definizione. Infine si mostra l'unicità della soluzione del problema. Sotto alcune ipotesi rafforzate sulla geometria del dominio, la soluzione ottenuta è classica. La terza parte della tesi riguarda il problema di Dirichlet associato a un'equazione della curvatura media prescritta nello spazio di Minkowski, che è di interesse in relatività generale. Il problema è ambientato in un dominio limitato regolare in $\mathbb{R}^N$ e un modello di curvatura media prescritta è dato da una funzione $f(x,s)$ che può avere comportamento sublineare, lineare, superlineare o sub-superlineare in $s=0$. L'attenzione è rivolta all'esistenza e alla molteplicità di soluzioni positive del problema. Come il precedente, anche questo capitolo è suddiviso in tre sezioni, che trattano rispettivamente i casi unidimensionale, radiale e $N$-dimensionale in un dominio generale. Nel caso unidimensionale, viene impiegato un approccio di tipo mappa-tempo per studiare una semplice situazione autonoma. Nel caso radiale in una palla, la tecnica è variazionale e lo studio del funzionale associato al problema evidenzia l'esistenza di un punto critico (casi sublineare o lineare), o di due (caso superlineare), o di tre punti critici (caso sub-superlineare): ciascuno di questi è una soluzione positiva del problema. Infine, nel caso generale in dimensione $N$, si adotta un approccio topologico che permette di studiare il problema non variazionale, in cui la funzione $f$ può dipendere dal gradiente della soluzione. Più nel dettaglio, con un metodo di sotto- e sopra-soluzioni specificamente sviluppato per questo problema, proviamo vari risultati di esistenza, molteplicità e localizzazione, in relazione alla presenza di una singola sotto-soluzione, o di una singola sopra-soluzione, o di una coppia di sotto- e sopra-soluzione ordinate o non ordinate. L'Appendice chiude la tesi: qui sono raccolti vari strumenti matematici utilizzati nel corso del lavoro.
This thesis is devoted to the study of some differential models arising in fluid mechanics or general relativity and involving the mean curvature operators in the $N$-dimensional Euclidean or Minkowski spaces. In both cases the operators are quasilinear elliptic operators which do not satisfy the property of uniform ellipticity, the Euclidean mean curvature operator being degenerate, whereas the Minkowski mean curvature operator being singular. This work is subdivided into three parts. The first one concerns the study of the periodic solutions of the one-dimensional prescribed curvature equation in the Euclidean space, which models capillarity-type phenomena. According to the structure of the curvature operator and imposing a suitable behaviour at zero, or at infinity, of the prescribed curvature, we prove the existence of infinitely many arbitrarily small classical subharmonic solutions with suitable nodal properties, or bounded variation subharmonic solutions with arbitrarily large oscillations. The technique for the search of classical solutions is topological and relies on the use of the rotation number and on a generalization of the Poincaré-Birkhoff theorem; whereas the approach for the study of non-classical solutions is based on non-smooth critical point theory, namely on a mountain pass lemma set in the space of bounded variation functions. The second part of the thesis is devoted to the study of the homogeneous Dirichlet problem associated with an anisotropic prescribed mean curvature equation in the Euclidean space, which provides a model for describing the geometry of the human cornea. The problem is set in a bounded domain in $\mathbb{R}^N$ with Lipschitz boundary. This chapter is subdivided into three sections, which are focused on the one-dimensional, the radial and the general $N$-dimensional case, respectively. In the one-dimensional and in the radial case in a ball, we prove an existence and uniqueness result of classical solution, which also displays some additional qualitative properties. Here the techniques used are topological in nature. Finally, in the $N$-dimensional case, we prove the existence, the uniqueness and the regularity of a strong-type solution of the problem. In order to tackle the possible gradient blow-up phenomena, the approach is variational and the framework is the space of bounded variation functions. We first collect some preliminary results about the behaviour of the action functional associated with the problem; among them, we remark the importance of an approximation property. We then prove the existence and uniqueness of the global minimizer of the action functional, which is smooth in the interior but non necessarily on the boundary, and satisfies the problem in a suitable sense. We finally prove the uniqueness of solution. Under some strengthened assumptions on the geometry of the domain, the solution obtained is classical. The third part of the thesis deals with the Dirichlet problem associated with a prescribed mean curvature equation in the Minkowski space, which is of interest in general relativity. The problem is set in a bounded regular domain in $\mathbb{R}^N$ and a model prescribed curvature is given by a function $f(x,s)$ whose behaviour is sublinear, linear, superlinear or sub-superlinear at $s=0$. The attention is addressed towards the existence and the multiplicity of positive solutions of the problem. In parallel to the second part of the thesis, this chapter is subdivided into three sections, which are focused on the one-dimensional, the radial and the general $N$-dimensional case, respectively. In the one-dimensional case, a time-map approach is employed for treating a simple autonomous situation. In the radial case in a ball, the technique is variational and the study of the action functional associated with the problem evidences the existence of either one (sublinear or linear cases), or two (superlinear case), or three (sub-superlinear case) non-trivial critical points of the action functional: each of them is a positive solution of the problem. Finally, in the general $N$-dimensional case, we adopt a topological approach which allows to study the non-variational problem, where the function $f$ may also depend on the gradient of the solution. Namely, by a lower and upper solution method specifically developed for this problem, we prove several existence, multiplicity and localization results, in relation to the presence of a single lower solution, or a single upper solution, or a couple of ordered or non-ordered lower and upper solutions of the problem. The Appendix completes this thesis: here several mathematical tools that have been used to prove the results are collected.
XXVI Ciclo
1986

2011/2012
La presente tesi ha come principale obiettivo lo studio della variabilità di proprietà idrologiche in mezzi porosi, con particolare attenzione alla permeabilità. A tal fine, ci si avvale di un approccio che combina l'analisi di proprietà statistiche e di scaling applicata a dataset di permeabilità, con lo studio di risultati numerici di simulazioni di flusso alla microscala in mezzi porosi. Con la prima analisi è possibile caratterizzare variazioni di permeabilità alla scala di misura (tipicamente dell'ordine del centimetro), mentre la seconda analisi dà una descrizione dell'eterogeneità di permeabilità ad una scala inferiore (nell'ordine del millimetro), ottenuta risolvendo processi fisici alla scala dei pori e derivando le quantità integrali di interesse. L'analisi statistica e di scaling, effettuata sia su distribuzioni di permeabilità sintetiche, sia su dataset raccolti su campioni reali, avvalora la validità dei modelli truncated fractional Brownian motion (tfBm) e truncated fractional Gaussian noise (tfGn), o di processi random sub-Gaussiani ad essi subordinati, per l'interpretazione della variabilità di proprietà idrologiche. Soluzioni numeriche di campi di flusso (i.e. velocità e pressione) alla scala dei pori sono ottenute sia per campioni sintetici, sia per campioni reali, la cui geometria è ricostruita mediante micro-tomografia a raggi X. Diverse metodologie di applicazione delle condizioni al contorno in corrispondenza dell'interfaccia liquido-solido forniscono risultati qualitativamente simili sia in termini di quantità microscopiche, sia in termini di quantità medie.
The work is aimed at providing some insights on the variability of hydrological properties in porous media, focusing in particular on permeability. We consider an approach which combines scaling and statistical analyses of air-permeability datasets with pore-scale numerical simulations of flow through porous media. The former investigation allows to characterize permeability heterogeneity at the centimeter observation scale; the latter provides a description of heterogeneity on a millimeter scale by resolving physical processes occurring at the microscopic scale and deriving up-scaled quantities. Scaling and statistical analyses performed on synthetic permeability distributions as well as on datasets collected on real media support the identification of truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn) and of sub-Gaussian random processes subordinated to tfBm (or tfGn) as viable models for the interpretation of hydrological properties variability. Pore-scale numerical solutions of flow (i.e., in terms of velocity and pressure distributions) are performed on both randomly generated samples and real porous media reconstructed via X-ray Micro-Tomography. Different approaches for the enforcement of boundary conditions at the fluid-solid interface provide qualitatively similar results in terms of both microscopic and averaged quantities.
XXV Ciclo
1984

2013/2014
Il sistema Adriatico Ionio è un importante componente nella circolazione termoalina del bacino orientale mediterraneo. Il mar Adriatico è la più importante sorgente di acqua profonde per il bacino. Il mar Ionio è una sorta di punto di incrocio per acque con caratteristiche diverse : le acque atlantiche in superficie , quelle levantine negli strati intermedi , le acque profonde adriatiche in quelli sul fondo . La variabilità osservata nella circolazione Ionica , con reversal periodici da stato ciclonico a stato anticiclonico e e viceversa, negli ultimi 20 anni è stata oggetto di dibattito in quanto è stata attribuita rispettivamente all'influenza della variazione delle proprietà termoaline delle acque profonde prodotte nel sud Adriatico sul bilancio di vorticità dello Ionio o a variazioni nel rotore del wind stress sul bacino ionico.Questo lavoro di tesi si prefigge di esaminare tutte le ipotesi relative alla dinamica del sistema adriatico-ionio da un punto di vista modellistico e la sua risposta in termini di bilancio di vorticità e di energia all'azione di forzanti esterne come il wind stress e i flussi termoalini. I risultati finali di questo approccio modellistico hanno mostrato la maggiore importanza del bacino est (mar cretese e passaggio di creta) rispetto all'Adriatico nel determinare la variabilità , a livello di bilancio di vorticità e di energia , della circolazione ionica .
The Adriatic-Ionian system represents an important component in the Mediterranean thermohaline circulation.The Adriatic Sea is the main source of deep waters for the Eastern Basin.The Ionian Sea is a sort of cross-road points for different kind of waters: Atlantic water moving on its surface layers , Levantine waters moving in its intermediate layers, Adriatic deep waters moving on its bottom layers.The variability observed in the Ionian upper layers circulation , which revers periodically from cyclonic to anticyclonic pattern and viceversa, in the last 20 years has been object of debate as it has been attributed to the effects of variation of thermohaline properties of deep waters produced in the Adriatic Sea on the vorticity balance of the Ionian or to variation in the wind stress curl over the Ionian basin.This thesis work aims to challenge all the hypotheses related the Adriatic-Ionian System dynamics using a modeling approach and its behavior in terms of vorticity and energy balance to some forcing as wind stress and thermohaline fluxes.The final results of this analysis have shown the predominant role , with respect to the Adriatic Sea , of Eastern basin in determining the vorticity and energy balance of the Ionian and shaping the variability of its circulation.
XXVII Ciclo
1982