Electrical spectroscopy of rubrene bulk and thin film crystals

One of the most diffused electronic device is the field effect transistor (FET), contained in number of billions in each electronic device. Organic optoelectronics is an emerging field that exploits the unique properties of conjugated organic materials to develop new applications that require a combination of performance, low cost and processability. Organic single crystals are the material with best performances and purity among the variety of different form of organic semiconductors. This thesis is focused on electrical and optical characterization of Rubrene single crystal bulk and thin films. Rubrene bulk is well known but for the first time we studied thin films. The first Current-voltage characterization has been performed for the first time on three Rubrene thin films with three different thickness to extract the charge carriers mobility and to assess its crystalline structure. As results we see that mobility increase with thickness. Field effect transistor based on Rubrene thin films on $SiO_2$ have been characterize by current-voltage (I-V) analyses (at several temperatures) and reveals a hopping conduction. Hopping behavior probably is due to the lattice mismatch with the substrate or intrinsic defectivity of the thin films. To understand effects of contact resistance we tested thin films with the Transmission Line Method (TLM) method. The TLM method revealeds that contact resistance is negligible but evidenced a Schottky behavior in a limited but well determined range of T. To avoid this effect we carried out annealing treatment after the electrode evaporation iswe performed a compete I-V characterization as a function of in temperature to extract the electronic density of states (DOS) distribution through the Space Charge Limited Current (SCLC) method. The results show a DOS with an exponential trenddistribution, as expected. The measured mobility of thin films is about 0.1cm^2/Vs and it increases with the film thickness. Further studies are necessary to investigate the reason and improve performances. From photocurrent spectrum we calculated an Eg of about 2.2eV and both thin films and bulk have a good crystal order. Further measurement are necessary to solve some open problems

Numerical and semi-analytical models of sliding masses: application to the 1783 Scilla tsunamigenic landslide

The Scilla rock avalanche occurred on 6 February 1783 along the coast of the Calabria region (southern Italy), close to the Messina Strait. It was triggered by a mainshock of the Terremoto delle Calabrie seismic sequence, and it induced a tsunami wave responsible for more than 1500 casualties along the neighboring Marina Grande beach. The main goal of this work is the application of semi-analtycal and numerical models to simulate this event. The first one is a MATLAB code expressly created for this work that solves the equations of motion for sliding particles on a two-dimensional surface through a fourth-order Runge-Kutta method. The second one is a code developed by the Tsunami Research Team of the Department of Physics and Astronomy (DIFA) of the Bologna University that describes a slide as a chain of blocks able to interact while sliding down over a slope and adopts a Lagrangian point of view. A wide description of landslide phenomena and in particular of landslides induced by earthquakes and with tsunamigenic potential is proposed in the first part of the work. Subsequently, the physical and mathematical background is presented; in particular, a detailed study on derivatives discratization is provided. Later on, a description of the dynamics of a point-mass sliding on a surface is proposed together with several applications of numerical and analytical models over ideal topographies. In the last part, the dynamics of points sliding on a surface and interacting with each other is proposed. Similarly, different application on an ideal topography are shown. Finally, the applications on the 1783 Scilla event are shown and discussed.

Testing of subgrid scale (SGS) models for large-eddy simulation (LES) of turbulent channel flow

Sub-grid scale (SGS) models are required in order to model the influence of the unresolved small scales on the resolved scales in large-eddy simulations (LES), the flow at the smallest scales of turbulence. In the following work two SGS models are presented and deeply analyzed in terms of accuracy through several LESs with different spatial resolutions, i.e. grid spacings. The first part of this thesis focuses on the basic theory of turbulence, the governing equations of fluid dynamics and their adaptation to LES. Furthermore, two important SGS models are presented: one is the Dynamic eddy-viscosity model (DEVM), developed by \cite{germano1991dynamic}, while the other is the Explicit Algebraic SGS model (EASSM), by \cite{marstorp2009explicit}. In addition, some details about the implementation of the EASSM in a Pseudo-Spectral Navier-Stokes code \cite{chevalier2007simson} are presented. The performance of the two aforementioned models will be investigated in the following chapters, by means of LES of a channel flow, with friction Reynolds numbers $Re_\tau=590$ up to $Re_\tau=5200$, with relatively coarse resolutions. Data from each simulation will be compared to baseline DNS data. Results have shown that, in contrast to the DEVM, the EASSM has promising potentials for flow predictions at high friction Reynolds numbers: the higher the friction Reynolds number is the better the EASSM will behave and the worse the performances of the DEVM will be. The better performance of the EASSM is contributed to the ability to capture flow anisotropy at the small scales through a correct formulation for the SGS stresses. Moreover, a considerable reduction in the required computational resources can be achieved using the EASSM compared to DEVM. Therefore, the EASSM combines accuracy and computational efficiency, implying that it has a clear potential for industrial CFD usage.

Numerical simulations of the effects of residual stresses in adhesively bonded composite specimen

In this work the problem of performing a numerical simulation of quasi-static crack propagation within an adhesive layer of a bonded joint under Mode I loading affected by stress field changes due to thermal-chemical shrinkage induced by cure process is addressed. Secondly, a parametric study on fracture critical energy, cohesive strength and Young's modulus is performed. Finally, a particular case of adhesive layer stiffening is simulated in order to verify qualitatively the major effect.

In-situ detection of defect formation in organic flexible electronics by Kelvin Probe Force Microscopy

Organic semiconductor technology has attracted considerable research interest in view of its great promise for large area, lightweight, and flexible electronics applications. Owing to their advantages in processing and unique physical properties, organic semiconductors can bring exciting new opportunities for broad-impact applications requiring large area coverage, mechanical flexibility, low-temperature processing, and low cost. In order to achieve highly flexible device architecture it is crucial to understand on a microscopic scale how mechanical deformation affects the electrical performance of organic thin film devices. Towards this aim, I established in this thesis the experimental technique of Kelvin Probe Force Microscopy (KPFM) as a tool to investigate the morphology and the surface potential of organic semiconducting thin films under mechanical strain. KPFM has been employed to investigate the strain response of two different Organic Thin Film Transistor with active layer made by 6,13-bis(triisopropylsilylethynyl)-pentacene (TIPS-Pentacene), and Poly(3-hexylthiophene-2,5-diyl) (P3HT). The results show that this technique allows to investigate on a microscopic scale failure of flexible TFT with this kind of materials during bending. I find that the abrupt reduction of TIPS-pentacene device performance at critical bending radii is related to the formation of nano-cracks in the microcrystal morphology, easily identified due to the abrupt variation in surface potential caused by local increase in resistance. Numerical simulation of the bending mechanics of the transistor structure further identifies the mechanical strain exerted on the TIPS-pentacene micro-crystals as the fundamental origin of fracture. Instead for P3HT based transistors no significant reduction in electrical performance is observed during bending. This finding is attributed to the amorphous nature of the polymer giving rise to an elastic response without the occurrence of crack formation.

Estimation of the ephemerides and gravity fields of the Galilean moons through orbit determination of the JUICE mission

Jupiter and its moons are a complex dynamical system that include several phenomenon like tides interactions, moon's librations and resonances. One of the most interesting characteristics of the Jovian system is the presence of the Laplace resonance, where the orbital periods of Ganymede, Europa and Io maintain a 4:2:1 ratio respectively. It is interesting to study the role of the Laplace Resonance in the dynamic of the system, especially regarding the dissipative nature of the tidal interaction between Jupiter and its closest moon, Io. Numerous theories have been proposed regarding the orbital evolution of the Galilean satellites, but they disagree about the amount of dissipation of the system, therefore about the magnitude and the direction of the evolution of the system, mainly because of the lack of experimental data. The future JUICE space mission is a great opportunity to solve this dispute. JUICE is an ESA (European Space Agency) L-class mission (the largest category of missions in the ESA Cosmic Vision) that, at the beginning of 2030, will be inserted in the Jovian system and that will perform several flybys of the Galilean satellites, with the exception of Io. Subsequently, during the last part of the mission, it will orbit around Ganymede for nine months, with a possible extension of the mission. The data that JUICE will collect during the mission will have an exceptional accuracy, allowing to investigate several aspects of the dynamics the system, especially, the evolution of Laplace Resonance of the Galilean moons and its stability. This thesis will focus on the JUICE mission, in particular in the gravity estimation and orbit reconstruction of the Galilean satellites during the Jovian orbital phase using radiometric data. This is accomplished through an orbit determination technique called multi-arc approach, using the JPL's orbit determination software MONTE (Mission-analysis, Operations and Navigation Tool-kit Environment).

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.

Numerical prediction of the lateral displacement capacity of flexural reinforced concrete shear walls

Previous earthquakes showed that shear wall damage could lead to catastrophic failures of the reinforced concrete building. The lateral load capacity of shear walls needs to be estimated to minimize associated losses during catastrophic events; hence it is necessary to develop and validate reliable and stable numerical methods able to converge to reasonable estimations with minimum computational effort. The beam-column 1-D line element with fiber-type cross-section model is a practical option that yields results in agreement with experimental data. However, shortcomings of using this model to predict the local damage response may come from the fact that the model requires fine calibration of material properties to overcome regularization and size effects. To reduce the mesh-dependency of the numerical model, a regularization method based on the concept of post-yield energy is applied in this work to both the concrete and the steel material constitutive laws to predict the nonlinear cyclic response and failure mechanism of concrete shear walls. Different categories of wall specimens known to produce a different response under in plane cyclic loading for their varied geometric and detailing characteristics are considered in this study, namely: 1) scaled wall specimens designed according to the European seismic design code and 2) unique full-scale wall specimens detailed according to the U.S. design code to develop a ductile behavior under cyclic loading. To test the boundaries of application of the proposed method, two full-scale walls with a mixed shear-flexure response and different values of applied axial load are also considered. The results of this study show that the use of regularized constitutive models considerably enhances the response predictions capabilities of the model with regards to global force-drift response and failure mode. The simulations presented in this thesis demonstrate the proposed model to be a valuable tool for researchers and engineers.

Experimental analysis of the aerodynamic and mixing performance of a ventilation device

Numerous types of acute respiratory failure are routinely treated using non-invasive ventilatory support (NIV). Its efficacy is well documented: NIV lowers intubation and death rates in various respiratory disorders. It can be delivered by means of face masks or head helmets. Currently the scientific community’s interest about NIV helmets is mostly focused on optimising the mixing between CO2 and clean air and on improving patient comfort. To this end, fluid dynamic analysis plays a particularly important role and a two- pronged approach is frequently employed. While on one hand numerical simulations provide information about the entire flow field and different geometries, they exhibit require huge temporal and computational resources. Experiments on the other hand help to validate simulations and provide results with a much smaller time investment and thus remain at the core of research in fluid dynamics. The aim of this thesis work was to develop a flow bench and to utilise it for the analysis of NIV helmets. A flow test bench and an instrumented mannequin were successfully designed, produced and put into use. Experiments were performed to characterise the helmet interface in terms of pressure drop and flow rate drop over different inlet flow rates and outlet pressure set points. Velocity measurements by means of Particle Image Velocimetry were performed. Pressure drop and flow rate characteristics from experiments were contrasted with CFD data and sufficient agreement was observed between both numerical and experimental results. PIV studies permitted qualitative and quantitative comparisons with numerical simulation data and offered a clear picture of the internal flow behaviour, aiding the identification of coherent flow features.

Simulation and Implementation of a PMSM Sensorless Control Algorithm for Submarine Applications

In this thesis, the study and the simulation of two advanced sensorless speed control techniques for a surface PMSM are presented. The aim is to implement a sensorless control algorithm for a submarine auxiliary propulsion system. This experimental activity is the result of a project collaboration with L3Harris Calzoni, a leader company in A&D systems for naval handling in military field. A Simulink model of the whole electric drive has been developed. Due to the satisfactory results of the simulations, the sensorless control system has been implemented in C code for STM32 environment. Finally, several tests on a real brushless machine have been carried out while the motor was connected to a mechanical load to simulate the real scenario of the final application. All the experimental results have been recorded through a graphical interface software developed at Calzoni.

Numerical modeling on the short-term response of fiber-reinforced concrete round panels

Fiber-reinforced concrete is a composite material consisting of discrete, discontinuous, and uniformly distributed fibers in plain concrete primarily used to enhance the tensile properties of the concrete. FRC performance depends upon the fiber, interface, and matrix properties. The use of fiber-reinforced concrete has been increasing substantially in the past few years in different fields of the construction industry such as ground-level application in sidewalks and building floors, tunnel lining, aircraft parking, runways, slope stabilization, etc. Many experiments have been performed to observe the short-term and long-term mechanical behavior of fiber-reinforced concrete in the last decade and numerous numerical models have been formulated to accurately capture the response of fiber-reinforced concrete. The main purpose of this dissertation is to numerically calibrate the short-term response of the concrete and fiber parameters in mesoscale for the three-point bending test and cube compression test in the MARS framework which is based on the lattice discrete particle model (LDPM) and later validate the same parameters for the round panels. LDPM is the most validated theory in mesoscale theories for concrete. Different seeds representing the different orientations of concrete and fiber particles are simulated to produce the mean numerical response. The result of numerical simulation shows that the lattice discrete particle model for fiber-reinforced concrete can capture results of experimental tests on the behavior of fiber-reinforced concrete to a great extent.