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.

WEC-Sim Application for the Assessment of Wave Energy Generation in the Tyrrhenian Sea and the North Sea Through the 2-Body Point Absorber (2BPA) and the Oscillating Surge Wave Energy Converter (OSWEC)

In recent years, developed countries have turned their attention to clean and renewable energy, such as wind energy and wave energy that can be converted to electrical power. Companies and academic groups worldwide are investigating several wave energy ideas today. Accordingly, this thesis studies the numerical simulation of the dynamic response of the wave energy converters (WECs) subjected to the ocean waves. This study considers a two-body point absorber (2BPA) and an oscillating surge wave energy converter (OSWEC). The first aim is to mesh the bodies of the earlier mentioned WECs to calculate their hydrostatic properties using axiMesh.m and Mesh.m functions provided by NEMOH. The second aim is to calculate the first-order hydrodynamic coefficients of the WECs using the NEMOH BEM solver and to study the ability of this method to eliminate irregular frequencies. The third is to generate a *.h5 file for 2BPA and OSWEC devices, in which all the hydrodynamic data are included. The BEMIO, a pre-and post-processing tool developed by WEC-Sim, is used in this study to create *.h5 files. The primary and final goal is to run the wave energy converter Simulator (WEC-Sim) to simulate the dynamic responses of WECs studied in this thesis and estimate their power performance at different sites located in the Mediterranean Sea and the North Sea. The hydrodynamic data obtained by the NEMOH BEM solver for the 2BPA and OSWEC devices studied in this thesis is imported to WEC-Sim using BEMIO. Lastly, the power matrices and annual energy production (AEP) of WECs are estimated for different sites located in the Sea of Sicily, Sea of Sardinia, Adriatic Sea, Tyrrhenian Sea, and the North Sea. To this end, the NEMOH and WEC-Sim are still the most practical tools to estimate the power generation of WECs numerically.

On the characterization of the mechanical behaviour of bonded steel-CFRP adhesive joints

The use of adhesively bonded carbon fiber reinforced polymers (CFRP) is well established to repair metallic structural elements in the aerospace industry for more than three decades. Despite a few exceptions, this technology has yet not been exploited for the steel construction industry where there is a great need to rehabilitate old metallic bridges. For instance, in Europe more than 30% of the railway bridge stock operated for more than 100 years. These bridges are made of old mild steel or puddle iron that exhibits poor behaviour due to the quality of the material itself and degradation caused by the long-term loading or environmental effects. The modest results for Steel/CFRP joints obtained may be due to the type of adhesive used. In fact, most of the previous studies utilized brittle adhesives specially developed for concrete structures. Recent ductile adhesives that made for the automotive industry for metallic joints should be more appropriate. In this study, an experimental investigation on the behaviour of CFRP/steel adhesively bonded joints is presented. A comparison between brittle adhesives and ductile adhesives is conducted. The results show that the ductile adhesives achieve much higher performance than the brittle ones. The brittle adhesives provide more stiffness to the adhesive joint. In the specimens with the ductile adhesives, the failure pattern started by yielding the steel bars first then the adhesive joint which is promising since it can facilitate the design significantly if the steel yielding can be used as a design criterion. The main disadvantage of ductile adhesives is they are usually more expensive than brittle ones. In order to solve this issue, bi-adhesive joints, in which the joint is mainly made of (low cost) brittle adhesive and ductile adhesive in the stress concentration region, are proposed. The results revealed very high improvement up to the yielding strength of the steel bars and with a balanced stiffness.

Root cause analysis applied to a finite element model's refinement of a negative stiffness structure

Negative Stiffness Structures are mechanical systems that require a decrease in the applied force to generate an increase in displacement. They are structures that possess special characteristics such as snap-through and bi-stability. All of these features make them particularly suitable for different applications, such as shock-absorption, vibration isolation and damping. From this point of view, they have risen awareness of their characteristics and, in order to match them to the application needed, a numerical simulation is of great interest. In this regard, this thesis is a continuation of previous studies in a circular negative stiffness structure and aims at refine the numerical model by presenting a new solution. To that end, an investigation procedure is needed. Amongst all of the methods available, root cause analysis was the chosen one to perform the investigation since it provides a clear view of the problem under analysis and a categorization of all the causes behind it. As a result of the cause-effect analysis, the main causes that have influence on the numerical results were obtained. Once all of the causes were listed, solutions to them were proposed and it led to a new numerical model. The numerical model proposed was of nonlinear type of analysis with hexagonal elements and a hyperelastic material model. The results were analyzed through force-displacement curves, allowing for the visualization of the structure’s energy recovery. When compared to the results obtained from the experimental part, it is evident that the trend is similar and the negative stiffness behaviour is present.

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.

Analysis of the Mediterranean tropical-like cyclone Numa from a numerical model simulation

I Medicane sono rari cicloni che si sviluppano sul Mar Mediterraneo e presentano caratteristiche dei cicloni tropicali, come la forma a spirale delle bande di nubi, un occhio privo di venti e nubi, venti intensi nella banda che circonda l’occhio e la presenza di un nucleo caldo. Nel presente lavoro, è stata compiuta un’analisi del Medicane Numa, verificatosi nel novembre del 2017, utilizzando gli output della simulazione del modello RAMS-ISAC. L’obiettivo di questa tesi è l’identificazione e la descrizione delle caratteristiche tropicali di Numa, focalizzandosi sulla descrizione dei diversi stadi di sviluppo del ciclone. Il sistema di bassa pressione è stato identificato utilizzando la pressione sul livello del mare, mentre l’anomalia termica e il vento orizzontale hanno permesso una descrizione della struttura del Medicane e del suo nucleo caldo. Un’illustrazione della struttura a spirale delle bande di nubi e dell’occhio è stata ottenuta con il grafico dei rapporti di mescolanza delle idrometeore di nubi e pioggia. Questi parametri hanno consentito di ricavare il diametro dell’occhio, pari a 75 km, mentre il diametro del Medicane è risultato 230 km. Numa ha registrato una velocità massima del vento in superficie di 20 m/s nella banda adiacente all’occhio del ciclone. Il diagramma di Hart dello spazio delle fasi del ciclone ha confermato la natura simil-tropicale di Numa e ne ha descritto l’evoluzione, identificando la transizione da sistema a caratteri tropicali a sistema ibrido. La traiettoria nello spazio delle fasi ha consentito l’identificazione delle sottofasi dell’evoluzione di Numa, confermate dai grafici dell’evoluzione temporale dei parametri menzionati in precedenza. L’analisi ha mostrato il ruolo cruciale della presenza di una struttura organizzata nel determinare l’intensità e la durata delle caratteristiche tropicali. Tutti i parametri hanno evidenziato la simmetria della struttura durante la persistente fase matura di Numa.