Fractional Order Modeling and Control of a Flexible Satellite

Researchers have engrossed fractional-order modeling because of its ability to capture phenomena that are nearly impossible to describe owing to its long-term memory and inherited properties. Motivated by the research in fractional modeling, a fractional-order prototype for a flexible satellite whose dynamics are governed by fractional differential equations is proposed for the first time. These relations are derived using fractional attitude dynamic description of rigid body simultaneously coupled with the fractional Lagrange equation that governs the vibration of the appendages. Two attitude controls are designed in the presence of the faults and uncertainties of the system. The first is the fractional-order feedback linearization controller, in which the stability of the internal dynamics of the system is proved. The second is the fractional-order sliding mode control, whose asymptotic stability is demonstrated using the quadratic Lyapunov function. Several nonlinear simulations are implemented to analyze the performance of the proposed controllers.

N-Heterocyclic carbene complexes of silver, rhodium and iron: structures, dynamics and catalysis

The research performed in the framework of this Master Thesis has been directly inspired by the recent work of an organometallic research group led by Professor Maria Cristina Cassani on a topic related to the structures, dynamics and catalytic activity of N-heterocyclic carbene-amide rhodium(I) complexes1. A series of [BocNHCH2CH2ImR]X (R = Me, X = I, 1a’; R = Bz, X = Br, 1b’; R = trityl, X = Cl, 1c’) amide-functionalized imidazolium salts bearing increasingly bulky N-alkyl substituents were synthetized and characterized. Subsequently, these organic precursors were employed in the synthesis of silver(I) complexes as intermediate compounds on a way to rhodium(I) complexes [Rh(NBD)X(NHC)] (NHC = 1-(2-NHBoc-ethyl)-3-R-imidazolin-2-ylidene; X = Cl, R = Me (3a’), R = Bz (3b’), R = trityl (3c’); X = I, R = Me (4a’)). VT NMR studies of these complexes revealed a restricted rotation barriers about the metal-carbene bond. However, while the rotation barriers calculated for the complexes in which R = Me, Bz (3a’,b’ and 4a) matched the experimental values, this was not true in the trityl case 3c’, where the experimental value was very similar to that obtained for compound 3b’ and much smaller with respect to the calculated one. In addition, the energy barrier derived for 3c’ from line shape simulation showed a strong dependence on the temperature, while the barriers measured for 3a’,b’ did not show this effect. In view of these results and in order to establish the reasons for the previously found inconsistency between calculated and experimental thermodynamic data, the first objective of this master thesis was the preparation of a series of rhodium(I) complexes [Rh(NBD)X(NHC)] (NHC = 1-benzyl-3-R-imidazolin-2-ylidene; X = Cl, R = Me, Bz, trityl, tBu), containing the benzyl substituent as a chiral probe, followed by full characterization. The second objective of this work was to investigate the catalytic activity of the new rhodium compounds in the hydrosilylation of terminal alkynes for comparison purposes with the reported complexes. Another purpose of this work was to employ the prepared N-heterocyclic ligands in the synthesis of iron(II)-NHC complexes.

Polymer-liquid crystal interface: a molecular dynamics study

Liquid crystals (LCs) are an interesting class of soft condensed matter systems characterized by an unusual combination of fluidity and long-range order, mainly known for their applications in displays (LCDs). However, the interest in LC continues to grow pushed by their application in new technologies in medicine, optical imaging, micro and nano technologies etc. In LCDs uniaxial alignment of LCs is mainly achieved by a rubbing process. During this treatment, the surfaces of polymer coated display substrates are rubbed in one direction by a rotating cylinder covered with a rubbing cloth. Basically, LC alignment involves two possible aligning directions: uniaxial planar (homogeneous) and vertical (homeotropic) to the display substrate. An interesting unresolved question concerning LCs regards the origin of their alignment on rubbed surfaces, and in particular on the polymeric ones used in the display industry. Most studies have shown that LCs on the surface of the rubbed polymer film layer are lying parallel to the rubbing direction. In these systems, micrometric grooves are generated on the film surface along the rubbing direction and also the polymer chains are stretched in this direction. Both the parallel aligned microgrooves and the polymer chains at the film surface may play a role in the LC alignment and it is not easy to quantify the effect of each contribution. The work described in this thesis is an attempt to find new microscopic evidences on the origin of LC alignment on polymeric surfaces through molecular dynamics (MD) simulations, which allow the investigation of the phenomenon with atomic detail. The importance of the arrangement of the polymeric chains in LCs alignment was studied by performing MD simulations of a thin film of a typical nematic LC, 4-cyano-4’-pentylbiphenyl (5CB), in contact with two different polymers: poly(methyl methacrylate)(PMMA) and polystyrene (PS). At least four factors are believed to influence the LC alignment: 1. the interactions of LCs with the backbone vinyl chains; 2. the interactions of LCs with the oriented side groups; 3. the anisotropic interactions of LCs with nanometric grooves; 4. the presence of static surface charges. Here we exclude the effect of microgrooves and of static surface charges from our virtual experiment, by using flat and neutral polymer surfaces, with the aim of isolating the chemical driving factors influencing the alignment of LC phases on polymeric surfaces.

Dynamics of Bora wind over the Adriatic sea: atmospheric water balance and role of air-sea fluxes and orography

The Bora wind is a mesoscale phenomenon which typically affects the Adriatic Sea basin for several days each year, especially during winter. The Bora wind has been studied for its intense outbreak across the Dinaric Alps. The properties of the Bora wind are widely discussed in the literature and scientific papers usually focus on the eastern Adriatic coast where strong turbulence and severe gust intensity are more pronounced. However, the impact of the Bora wind can be significant also over Italy, not only in terms of wind speed instensity. Depending on the synoptic pressure pattern (cyclonic or anticyclonic Bora) and on the season, heavy snowfall, severe storms, storm surges and floods can occur along the Adriatic coast and on the windward flanks of the Apennines. In the present work five Bora cases that occurred in recent years have been selected and their evolution has been simulated with the BOLAM-MOLOCH model set, developed at ISAC-CNR in Bologna. Each case study has been addressed by a control run and by several sensitivity tests, performed with the purpose of better understanding the role played by air-sea latent and sensible heat fluxes. The tests show that the removal of the fluxes induces modifications in the wind approching the coast and a decrease of the total precipitation amount predicted over Italy. In order to assess the role of heat fluxes, further analysis has been carried out: column integrated water vapour fluxes have been computed along the Italian coastline and an atmospheric water balance has been evaluated inside a box volume over the Adriatic Sea. The balance computation shows that, although latent heat flux produces a significant impact on the precipitation field, its contribution to the balance is relatively minor. The most significant and lasting case study, that of February 2012, has been studied in more detail in order to explain the impressive drop in the total precipitation amount simulated in the sensitivity tests with removed heat fluxes with respect to the CNTRL run. In these experiments relative humidity and potential temperature distribution over different cross-sections have been examined. With respect to the CNTRL run a drier and more stable boundary layer, characterised by a more pronounced wind shear at the lower levels, has been observed to establish above the Adriatic Sea. Finally, in order to demonstrate that also the interaction of the Bora flow with the Apennines plays a crucial role, sensitivity tests varying the orography height have been considered. The results of such sensitivity tests indicate that the propagation of the Bora wind over the Adriatic Sea, and in turn its meteorological impact over Italy, is influenced by both the large air-sea heat fluxes and the interaction with the Apennines that decelerate the upstream flow.

Generic wind estimation and compensation based on residual generators and higher-order sliding mode schemes

This master thesis proposes a solution to the approach problem in case of unknown severe microburst wind shear for a fixed-wing aircraft, accounting for both longitudinal and lateral dynamics. The adaptive controller design for wind rejection is also addressed, exploiting the wind estimation provided by suitable estimators. It is able to successfully complete the final approach phase even in presence of wind shear, and at the same time aerodynamic envelope protection is retained. The adaptive controller for wind compensation has been designed by a backstepping approach and feedback linearization for time-varying systems. The wind shear components have been estimated by higher-order sliding mode schemes. At the end of this work the results are provided, an autonomous final approach in presence of microburst is discussed, performances are analyzed, and estimation of the microburst characteristics from telemetry data is examined.

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.

Characterization of Mars rotational dynamics from Doppler tracking of the InSight lander

In the last decades the evolution of radio science has made it possible to infer the atmosphere composition, the surface and the internal structure of the planets. Since the arrival of the first landers on Mars it was possible to make accurate measurements of the dynamics of this planet; in this thesis we will focus on InSight, considering the data disclosed by the JPL relative to the period from November 26th, 2018 to August 15th, 2021. In particular, the Doppler and Range measurements conducted by the RISE (Rotation and Interior Structure Experiment) will be analyzed. Since the accuracy of these measurements was improved significantly the effects due to the atmosphere of Mars might be measured so it should thus be possible to obtain a better estimate of the parameters characterizing the rotational dynamic of Mars. A large part of this study will therefore be dedicated to the study, modeling, implementation and analysis of the atmosphere of Mars, in both its components: troposphere and ionosphere. Once the complete model of Mars had been built, i.e. including the atmosphere, it was then possible to analyze the residuals, obtained between the data of the measurements carried out and the values predicted by the developed model, in order to obtain an estimate of the rotational dynamic of Mars.