Large Eddy Simulations applied to wind loading and pollutant dispersion

The application of Computational Fluid Dynamics based on the Reynolds-Averaged Navier-Stokes equations to the simulation of bluff body aerodynamics has been thoroughly investigated in the past. Although a satisfactory accuracy can be obtained for some urban physics problems their predictive capability is limited to the mean flow properties, while the ability to accurately predict turbulent fluctuations is recognized to be of fundamental importance when dealing with wind loading and pollution dispersion problems. The need to correctly take into account the flow dynamics when such problems are faced has led researchers to move towards scale-resolving turbulence models such as Large Eddy Simulations (LES). The development and assessment of LES as a tool for the analysis of these problems is nowadays an active research field and represents a demanding engineering challenge. This research work has two objectives. The first one is focused on wind loads assessment and aims to study the capabilities of LES in reproducing wind load effects in terms of internal forces on structural members. This differs from the majority of the existing research, where performance of LES is evaluated only in terms of surface pressures, and is done with a view of adopting LES as a complementary design tools alongside wind tunnel tests. The second objective is the study of LES capabilities in calculating pollutant dispersion in the built environment. The validation of LES in this field is considered to be of the utmost importance in order to conceive healthier and more sustainable cities. In order to validate the numerical setup adopted, a systematic comparison between numerical and experimental data is performed. The obtained results are intended to be used in the drafting of best practice guidelines for the application of LES in the urban physics field with a particular attention to wind load assessment and pollution dispersion problems.

Experimental and Numerical Analysis of Gas Flows in Microchannels and Micro Heat Exchanger

Due to increased interest in miniaturization, great attention has been given in the recent decade to the micro heat exchanging systems. Literature survey suggests that there is still a limited understanding of gas flows in micro heat exchanging systems. The aim of the current thesis is to further the understanding of fluid flow and heat transfer phenomenon inside such geometries when a compressible working fluid is utilized. A combined experimental and numerical approach has been utilized in order to overcome the lack of employable sensors for micro dimensional channels. After conducting a detailed comparison between various data reduction methodologies employed in the literature, the best suited methodology for gas microflow experimentalists is proposed. A transitional turbulence model is extensively validated against the experimental results of the microtubes and microchannels under adiabatic wall conditions. Heat transfer analysis of single microtubes showed that when the compressible working fluid is used, Nusselt number results are in partial disagreement with the conventional theory at highly turbulent flow regime for microtubes having a hydraulic diameter less than 250 microns. Experimental and numerical analysis on a prototype double layer microchannel heat exchanger showed that compressibility is detrimental to the thermal performance. It has been found that compressibility effects for micro heat exchangers are significant when the average Mach number at the outlet of the microchannel is greater than 0.1 compared to the adiabatic limit of 0.3. Lastly, to avoid a staggering amount of the computational power needed to simulate the micro heat exchanging systems with hundreds of microchannels, a reduced order model based on the porous medium has been developed that considers the compressibility of the gas inside microchannels. The validation of the proposed model against experimental results of average thermal effectiveness and the pressure loss showed an excellent match between the two.

Experimental determination of burning velocity in metal dust explosions

The modeling of metal dust explosion phenomenon is important in order to safeguard industries from potential accidents. A key parameter of these models is the burning velocity, which represents the consumption rate of the reactants by the flame front, during the combustion process. This work is focused on the experimental determination of aluminium burning velocity, through an alternative method, called "Direct method". The study of the methods used and the results obtained is preceded by a general analysis on dust explosion phenomenon, flame propagation phenomenon, characteristics of the metals combustion process and standard methods for determining the burning velocity. The “Direct method” requires a flame propagating through a tube recorded by high-speed cameras. Thus, the flame propagation test is carried out inside a vertical prototype made of glass. The study considers two optical technique: the direct visualization of the light emitted by the flame and the Particle Image Velocimetry (PIV) technique. These techniques were used simultaneously and allow the determination of two velocities: the flame propagation velocity and the flow velocity of the unburnt mixture. Since the burning velocity is defined by these two quantities, its direct determination is done by substracting the flow velocity of the fresh mixture from the flame propagation velocity. The results obtained by this direct determination, are approximated by a linear curve and different non-linear curves, which show a fluctuating behaviour of burning velocity. Furthermore, the burning velocity is strongly affected by turbulence. Turbulence intensity can be evaluated from PIV technique data. A comparison between burning velocity and turbulence intensity highlighted that both have a similar trend.

Development of a numerical platform for the modeling and optimal control of liquid metal flows

The main purpose of this work is to develop a numerical platform for the turbulence modeling and optimal control of liquid metal flows. Thanks to their interesting thermal properties, liquid metals are widely studied as coolants for heat transfer applications in the nuclear context. However, due to their low Prandtl numbers, the standard turbulence models commonly used for coolants as air or water are inadequate. Advanced turbulence models able to capture the anisotropy in the flow and heat transfer are then necessary. In this thesis, a new anisotropic four-parameter turbulence model is presented and validated. The proposed model is based on explicit algebraic models and solves four additional transport equations for dynamical and thermal turbulent variables. For the validation of the model, several flow configurations are considered for different Reynolds and Prandtl numbers, namely fully developed flows in a plane channel and cylindrical pipe, and forced and mixed convection in a backward-facing step geometry. Since buoyancy effects cannot be neglected in liquid metals-cooled fast reactors, the second aim of this work is to provide mathematical and numerical tools for the simulation and optimization of liquid metals in mixed and natural convection. Optimal control problems for turbulent buoyant flows are studied and analyzed with the Lagrange multipliers method. Numerical algorithms for optimal control problems are integrated into the numerical platform and several simulations are performed to show the robustness, consistency, and feasibility of the method.

Non-thermal phenomena in galaxy clusters: the LOFAR revolution

Turbulence introduced into the intra-cluster medium (ICM) through cluster merger events transfers energy to non-thermal components (relativistic particles and magnetic fields) and can trigger the formation of diffuse synchrotron radio sources. Owing to their steep synchrotron spectral index, such diffuse sources can be better studied at low radio frequencies. In this respect, the LOw Frequency ARray (LOFAR) is revolutionizing our knowledge thanks to its unprecedented resolution and sensitivity below 200 MHz. In this Thesis we focus on the study of radio halos (RHs) by using LOFAR data. In the first part of this work we analyzed the largest-ever sample of galaxy clusters observed at radio frequencies. This includes 309 Planck clusters from the Second Data Release of the LOFAR Two Metre Sky Survey (LoTSS-DR2), which span previously unexplored ranges of mass and redshift. We detected 83 RHs, half of which being new discoveries. In 140 clusters we lack a detected RH; for this sub-sample we developed new techniques to derive upper limits to their radio powers. By comparing detections and upper limits, we carried out the first statistical analysis of populations of clusters observed at low frequencies and tested theoretical formation models. In the second part of this Thesis we focused on ultra-steep spectrum radio halos. These sources are almost undetected at GHz frequencies, but are thought to be common at low frequencies. We presented LOFAR observations of two interesting clusters hosting ultra-steep spectrum radio halos. With complementary radio and X-ray observations we constrained the properties and origin of these targets.

Gyrokinetic modelling of plasma transport and investigation on test particle dynamics

Turbulent plasmas inside tokamaks are modeled and studied using guiding center theory, applied to charged test particles, in a Hamiltonian framework. The equations of motion for the guiding center dynamics, under the conditions of a constant and uniform magnetic field and turbulent electrostatic field are derived by averaging over the fast gyroangle, for the first and second order in the guiding center potential, using invertible changes of coordinates such as Lie transforms. The equations of motion are then made dimensionless, exploiting temporal and spatial periodicities of the model chosen for the electrostatic potential. They are implemented numerically in Python. Fast Fourier Transform and its inverse are used. Improvements to the original Python scripts are made, notably the introduction of a power-law curve fitting to account for anomalous diffusion, the possibility to integrate the equations in two steps to save computational time by removing trapped trajectories, and the implementation of multicolored stroboscopic plots to distinguish between trapped and untrapped guiding centers. The post-processing of the results is made in MATLAB. The values and ranges of the parameters chosen for the simulations are selected based on numerous simulations used as feedback tools. In particular, a recurring value for the threshold to detect trapped trajectories is evidenced. Effects of the Larmor radius, the amplitude of the guiding center potential and the intensity of its second order term are studied by analyzing their diffusive regimes, their stroboscopic plots and the shape of guiding center potentials. The main result is the identification of cases anomalous diffusion depending on the values of the parameters (mostly the Larmor radius). The transitions between diffusive regimes are identified. The presence of highways for the super-diffusive trajectories are unveiled. The influence of the charge on these transitions from diffusive to ballistic behaviors is analyzed.

Experimental characterization of a planar jet by hot-wire anemometry

The study of turbulence is also nowadays a problem that does not have solution from the mathematical point of view due to the lack of solution to link the mean part of the flow with the fluctuating one. To solve this problem, in the CICLoPE laboratory of Predappio, experiments on different type of jets are performed in order to derive a closure model able to close our mathematical model. One of the most interesting type of jet that could be studied is the planar turbulent free jet which is a two dimensional canonical jet characterized by the self-similarity condition of the velocity profiles. To study this particular jet, a new facility was built. The aim of this project is to characterize the jet at different distances from the nozzle exit, for different values of Reynolds number, to demonstrate that the self-similarity condition is respected. To do that, the evaluation of quantities such as spreading rate, centerline velocity decay and relation between fluctuations and mean part of the flow has to be obtain. All these parameters could be detected thanks to the use of single and X hot-wire anemometry with which it is possible to analyzed the fluctuating behaviour of the flow by associating to an electric signal a physical variable expressed in terms of velocity. To justify the data obtain by the measures, a comparison with results coming from the literature has to be shown.

A comparison between vertical mixing parameterizations in the Levantine Sea

The goal of this thesis was the study of an optimal vertical mixing parameterization scheme in a mesoscale dominated field characterized from a strong vorticity and the presence of a layer of colder, less saline water at about 100 m depth (Atlantic Waters); in these conditions we compared six different experiments, that differ by the turbulent closure schemes, the presence or not of an enhanced diffusion parameterization and the presence or not of a double diffusion mixing parameterization. To evaluate the performance of the experiments and the model we compared the simulations with the ARGO observations of temperature and salinity available in our domain, in our period of interest. The conclusions were the following: • the increase of the resolution gives better results in terms of temperature in all the considered cases, and in terms of salinity. • The comparisons between the Pacanovski-Philander and the TKE turbulent closure schemes don’t show significant differences when the simulations are compared to the observations. • The removing of the enhanced diffusion parameterization in presence of the TKE turbulent closure submodel doesn’t give positive results, and show limitations in the resolving of gravitational instabilities near the surface • The k-ϵ turbulent closure model utilized in all the GLS experiments, is the best performing closure model among the three considered, with positive results in all the salinity comparison with the in situ observation and in most of the temperature comparisons. • The double mixing parameterization utilized in the k-ϵ closure submodel improves the results of the experiments improving both the temperature and salinity in comparison with the ARGO data.

Temperature measurements in turbulent pipe flow with cold wire anemometry

Passive scalars measurements in turbulent pipe flows are difficult to perform and only few experimental data are available in literature. The present thesis deals with the experimental acquisition and study of the first turbulent temperature profile inside the CICLoPE wind tunnel through cold wire anemometry technique at Reτ = 6000 and Reτ = 9500. This type of measurements provides not only useful data on temperature (and passive scalars) behaviour and statistics in turbulent pipe flows, but could be used also for temperature correction of turbulent velocity profiles. In the present work, subsequent acquisitions of temperature and velocity profiles has been performed at the same Reynolds number and in the same points, through cold wire and hot wire techniques respectively. Taking as reference data from both DNS and experimental campaigns, the activity has been carried out obtaining satisfactory results. We have verified the presence of turbulent temperature profile inside the CICLoPE wind tunnel and then studied its statistical and spectral behaviours obtaining results in agreement with existing data from Hishida, Nagano, and Ferro. Cold wire temperature data were then used to correct hot wire velocity data, obtaining a slightly improvement in the near wall region.

Impact of convective motions on pollutant dispersion in urban canyons: an analysis through RANS simulations

Since the majority of the population of the world lives in cities and that this number is expected to increase in the next years, one of the biggest challenges of the research is the determination of the risk deriving from high temperatures experienced in urban areas, together with improving responses to climate-related disasters, for example by introducing in the urban context vegetation or built infrastructures that can improve the air quality. In this work, we will investigate how different setups of the boundary and initial conditions set on an urban canyon generate different patterns of the dispersion of a pollutant. To do so we will exploit the low computational cost of Reynolds-Averaged Navier-Stokes (RANS) simulations to reproduce the dynamics of an infinite array of two-dimensional square urban canyons. A pollutant is released at the street level to mimic the presence of traffic. RANS simulations are run using the k-ɛ closure model and vertical profiles of significant variables of the urban canyon, namely the velocity, the turbulent kinetic energy, and the concentration, are represented. This is done using the open-source software OpenFOAM and modifying the standard solver simpleFoam to include the concentration equation and the temperature by introducing a buoyancy term in the governing equations. The results of the simulation are validated with experimental results and products of Large-Eddy Simulations (LES) from previous works showing that the simulation is able to reproduce all the quantities under examination with satisfactory accuracy. Moreover, this comparison shows that despite LES are known to be more accurate albeit more expensive, RANS simulations represent a reliable tool if a smaller computational cost is needed. Overall, this work exploits the low computational cost of RANS simulations to produce multiple scenarios useful to evaluate how the dispersion of a pollutant changes by a modification of key variables, such as the temperature.