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.

Ultra-Low Temperature Partial Oxidation: Detailed Kinetics, Safety, and Environmental Issues

The innovation in several industrial sectors has been recently characterized by the need for reducing the operative temperature either for economic or environmental related aspects. Promising technological solutions require the acquisition of fundamental-based knowledge to produce safe and robust systems. In this sense, reactive systems often represent the bottleneck. For these reasons, this work was focused on the integration of chemical (i.e., detailed kinetic mechanism) and physical (i.e., computational fluid dynamics) models. A theoretical-based kinetic mechanism mimicking the behaviour of oxygenated fuels and their intermediates under oxidative conditions in a wide range of temperature and pressure was developed. Its validity was tested against experimental data collected in this work by using the heat flux burner, as well as measurements retrieved from the current literature. Besides, estimations deriving from existing models considered as the benchmark in the combustion field were compared with the newly generated mechanism. The latter was found to be the most accurate for the investigated conditions and fuels. Most influential species and reactions on the combustion of butyl acetate were identified. The corresponding thermodynamic parameter and rate coefficients were quantified through ab initio calculations. A reduced detailed kinetic mechanism was produced and implemented in an open-source computational fluid dynamics model to characterize pool fires caused by the accidental release of aviation fuel and liquefied natural gas, at first. Eventually, partial oxidation processes involving light alkenes were optimized following the quick, fair, and smoot (QFS) paradigm. The proposed procedure represents a comprehensive and multidisciplinary approach for the construction and validation of accurate models, allowing for the characterization of developing industrial sectors and techniques.

Development of a patient-specific model for stroke risk assessment in atrial fibrillation patients

Atrial fibrillation is associated with a five-fold increase in the risk of cerebrovascular events,being responsible of 15-18% of all strokes.The morphological and functional remodelling of the left atrium caused by atrial fibrillation favours blood stasis and, consequently, stroke risk. In this context, several clinical studies suggest that stroke risk stratification could be improved by using haemodynamic information on the left atrium (LA) and the left atrial appendage (LAA). The goal of this study was to develop a personalized computational fluid-dynamics (CFD) model of the left atrium which could clarify the haemodynamic implications of atrial fibrillation on a patient specific basis. The developed CFD model was first applied to better understand the role of LAA in stroke risk. Infact, the interplay of the LAA geometric parameters such as LAA length, tortuosity, surface area and volume with the fluid-dynamics parameters and the effects of the LAA closure have not been investigated. Results demonstrated the capabilities of the CFD model to reproduce the real physiological behaviour of the blood flow dynamics inside the LA and the LAA. Finally, we determined that the fluid-dynamics parameters enhanced in this research project could be used as new quantitative indexes to describe the different types of AF and open new scenarios for the patient-specific stroke risk stratification.

Advanced voxel-based CAD modelling for FSI simulations for automotive structures design

Additive Manufacturing (AM) is nowadays considered an important alternative to traditional manufacturing processes. AM technology shows several advantages in literature as design flexibility, and its use increases in automotive, aerospace and biomedical applications. As a systematic literature review suggests, AM is sometimes coupled with voxelization, mainly for representation and simulation purposes. Voxelization can be defined as a volumetric representation technique based on the model’s discretization with hexahedral elements, as occurs with pixels in the 2D image. Voxels are used to simplify geometric representation, store intricated details of the interior and speed-up geometric and algebraic manipulation. Compared to boundary representation used in common CAD software, voxel’s inherent advantages are magnified in specific applications such as lattice or topologically structures for visualization or simulation purposes. Those structures can only be manufactured with AM employment due to their complex topology. After an accurate review of the existent literature, this project aims to exploit the potential of the voxelization algorithm to develop optimized Design for Additive Manufacturing (DfAM) tools. The final aim is to manipulate and support mechanical simulations of lightweight and optimized structures that should be ready to be manufactured with AM with particular attention to automotive applications. A voxel-based methodology is developed for efficient structural simulation of lattice structures. Moreover, thanks to an optimized smoothing algorithm specific for voxel-based geometries, a topological optimized and voxelized structure can be transformed into a surface triangulated mesh file ready for the AM process. Moreover, a modified panel code is developed for simple CFD simulations using the voxels as a discretization unit to understand the fluid-dynamics performances of industrial components for preliminary aerodynamic performance evaluation. The developed design tools and methodologies perfectly fit the automotive industry’s needs to accelerate and increase the efficiency of the design workflow from the conceptual idea to the final product.

Characterization of pyrolysis and oxidation processes of biomaterials, biomasses and biofuels

The urgent need for alternative solutions mitigating the impacts of human activities on the environment has strongly opened new challenges and opportunities in view of the energy transition. Indeed, the automotive industry is going through a revolutionary moment in its quest to reduce its carbon footprint, with biofuels being one of the viable alternatives. The use of different classes of biofuels as fuel additives/standalone components has attracted the attention of many researchers. Despite their beneficial effects, biofuel’s combustion can also result in the production of undesirable pollutants, requiring complete characterization of the phenomena occurring during their production and consumption. Industrial scale-up of biomass conversion is challenging owing to the complexity of its chemistry and transport phenomena involved in the process. In this view, the role of solid-phase and gas-phase chemistry is paramount. Thus, this study is devoted to detailed analysis of physical-chemical phenomena characterizing biomass pyrolysis and biofuel oxidation. The pyrolysis mechanism has been represented by 20 reactions whereas, the gas-phase kinetic models; manually upgraded model (KiBo_MU) and automated model (KiBo_AG), comprises 141 species and 453 reactions, and 631 species and 28329 reactions, respectively. The accuracy of the kinetic models was tested against experimental data and the models captured experimental trends very well. While the development and validation of detailed kinetic mechanisms is the main deliverable of this project, the realized procedure integrating schematic classifications with methodologies for the identification of common decomposition pathways and intermediates represents an additional source of novelty. Besides, the fundamentally oriented nature of the adopted method allows the identification of most relevant reactions and species under the operating conditions different industrial applications, paving the way for reduced kinetic mechanisms. Ultimately, the resulting detailed mechanisms can be used to integrate with more complex fluid dynamics model to accurately reproduce the behavior of real systems and reactors.

Decoding complex flow phenomena via Dynamic Mode Decomposition

In this thesis, the viability of the Dynamic Mode Decomposition (DMD) as a technique to analyze and model complex dynamic real-world systems is presented. This method derives, directly from data, computationally efficient reduced-order models (ROMs) which can replace too onerous or unavailable high-fidelity physics-based models. Optimizations and extensions to the standard implementation of the methodology are proposed, investigating diverse case studies related to the decoding of complex flow phenomena. The flexibility of this data-driven technique allows its application to high-fidelity fluid dynamics simulations, as well as time series of real systems observations. The resulting ROMs are tested against two tasks: (i) reduction of the storage requirements of high-fidelity simulations or observations; (ii) interpolation and extrapolation of missing data. The capabilities of DMD can also be exploited to alleviate the cost of onerous studies that require many simulations, such as uncertainty quantification analysis, especially when dealing with complex high-dimensional systems. In this context, a novel approach to address parameter variability issues when modeling systems with space and time-variant response is proposed. Specifically, DMD is merged with another model-reduction technique, namely the Polynomial Chaos Expansion, for uncertainty quantification purposes. Useful guidelines for DMD deployment result from the study, together with the demonstration of its potential to ease diagnosis and scenario analysis when complex flow processes are involved.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.

Linear and nonlinear thermal instability of Newtonian and non-Newtonian fluid saturated porous media

The present work aims to investigate the influence of different aspects, such as non-standard steady solutions, complex fluid rheologies and non-standard porous-channel geometries, on the stability of a Darcy-Bénard system. In order to do so, both linear and nonlinear stability theories are considered. A linear analysis focuses on studying the dynamics of the single disturbance wave present in the system, while its nonlinear counterpart takes into consideration the interactions among the single modes. The scope of the stability analysis is to obtain information regarding the transition from an equilibrium solution to another one, and also information regarding the transition nature and the emergent solution after the transition. The disturbance governing equations are solved analytically, whenever possible, and numerical by considering different approaches. Among other important results, it is found that a cylinder cross-section does not affect the thermal instability threshold, but just the linear pattern selection for dilatant and pseudoplastic fluid saturated porous media. A new rheological model is proposed as a solution for singular issues involving the power-law model. Also, a generalised class of one parameter basic solutions is proposed as an alternative description of the isoflux Darcy--Bénard problem. Its stability is investigated.