Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Design and development of new pressure sensors for aerodynamic applications

This artwork reports on two different projects that were carried out during the three years of Doctor of the Philosophy course. In the first years a project regarding Capacitive Pressure Sensors Array for Aerodynamic Applications was developed in the Applied Aerodynamic research team of the Second Faculty of Engineering, University of Bologna, Forlì, Italy, and in collaboration with the ARCES laboratories of the same university. Capacitive pressure sensors were designed and fabricated, investigating theoretically and experimentally the sensor’s mechanical and electrical behaviours by means of finite elements method simulations and by means of wind tunnel tests. During the design phase, the sensor figures of merit are considered and evaluated for specific aerodynamic applications. The aim of this work is the production of low cost MEMS-alternative devices suitable for a sensor network to be implemented in air data system. The last two year was dedicated to a project regarding Wireless Pressure Sensor Network for Nautical Applications. Aim of the developed sensor network is to sense the weak pressure field acting on the sail plan of a full batten sail by means of instrumented battens, providing a real time differential pressure map over the entire sail surface. The wireless sensor network and the sensing unit were designed, fabricated and tested in the faculty laboratories. A static non-linear coupled mechanical-electrostatic simulation, has been developed to predict the pressure versus capacitance static characteristic suitable for the transduction process and to tune the geometry of the transducer to reach the required resolution, sensitivity and time response in the appropriate full scale pressure input A time dependent viscoelastic error model has been inferred and developed by means of experimental data in order to model, predict and reduce the inaccuracy bound due to the viscolelastic phenomena affecting the Mylar® polyester film used for the sensor diaphragm. The development of the two above mentioned subjects are strictly related but presently separately in this artwork.

Creep and damage models for the service behaviour of structural members

Deformability is often a crucial to the conception of many civil-engineering structural elements. Also, design is all the more burdensome if both long- and short-term deformability has to be considered. In this thesis, long- and short-term deformability has been studied from the material and the structural modelling point of view. Moreover, two materials have been handled: pultruded composites and concrete. A new finite element model for thin-walled beams has been introduced. As a main assumption, cross-sections rigid are considered rigid in their plane; this hypothesis replaces that of the classical beam theory of plane cross-sections in the deformed state. That also allows reducing the total number of degrees of freedom, and therefore making analysis faster compared with twodimensional finite elements. Longitudinal direction warping is left free, allowing describing phenomena such as the shear lag. The new finite-element model has been first applied to concrete thin-walled beams (such as roof high span girders or bridge girders) subject to instantaneous service loadings. Concrete in his cracked state has been considered through a smeared crack model for beams under bending. At a second stage, the FE-model has been extended to the viscoelastic field and applied to pultruded composite beams under sustained loadings. The generalized Maxwell model has been adopted. As far as materials are concerned, long-term creep tests have been carried out on pultruded specimens. Both tension and shear tests have been executed. Some specimen has been strengthened with carbon fibre plies to reduce short- and long- term deformability. Tests have been done in a climate room and specimens kept 2 years under constant load in time. As for concrete, a model for tertiary creep has been proposed. The basic idea is to couple the UMLV linear creep model with a damage model in order to describe nonlinearity. An effective strain tensor, weighting the total and the elasto-damaged strain tensors, controls damage evolution through the damage loading function. Creep strains are related to the effective stresses (defined by damage models) and so associated to the intact material.

Kinematic description of the rupture from strong motion data: strategies for a robust inversion

We present a non linear technique to invert strong motion records with the aim of obtaining the final slip and rupture velocity distributions on the fault plane. In this thesis, the ground motion simulation is obtained evaluating the representation integral in the frequency. The Green’s tractions are computed using the discrete wave-number integration technique that provides the full wave-field in a 1D layered propagation medium. The representation integral is computed through a finite elements technique, based on a Delaunay’s triangulation on the fault plane. The rupture velocity is defined on a coarser regular grid and rupture times are computed by integration of the eikonal equation. For the inversion, the slip distribution is parameterized by 2D overlapping Gaussian functions, which can easily relate the spectrum of the possible solutions with the minimum resolvable wavelength, related to source-station distribution and data processing. The inverse problem is solved by a two-step procedure aimed at separating the computation of the rupture velocity from the evaluation of the slip distribution, the latter being a linear problem, when the rupture velocity is fixed. The non-linear step is solved by optimization of an L2 misfit function between synthetic and real seismograms, and solution is searched by the use of the Neighbourhood Algorithm. The conjugate gradient method is used to solve the linear step instead. The developed methodology has been applied to the M7.2, Iwate Nairiku Miyagi, Japan, earthquake. The estimated magnitude seismic moment is 2.6326 dyne∙cm that corresponds to a moment magnitude MW 6.9 while the mean the rupture velocity is 2.0 km/s. A large slip patch extends from the hypocenter to the southern shallow part of the fault plane. A second relatively large slip patch is found in the northern shallow part. Finally, we gave a quantitative estimation of errors associates with the parameters.

Evoluzione Morfotettonica delle Aree Alpine "Sempione" e "Brennero" attraverso Studi Termocronologici di bassa temperatura

E’ mostrata l’analisi e la modellazione di dati termocronologici di bassa temperatura da due regioni Alpine: il Sempione ed il Brennero. Le faglie distensive presenti bordano settori crostali profondi appartenenti al dominio penninico: il duomo metamorfico Lepontino al Sempione e la finestra dei Tauri al Brennero. I dati utilizzati sono FT e (U-Th)/He su apatite. Per il Sempione i dati provengono dalla bibliografia; per il Brennero si è provveduto ad un nuovo campionamento, sia in superficie che in sotterraneo. Gli attuali lavori per la galleria di base del Brennero (BBT), hanno consentito, per la prima volta, di raccogliere dati di FT e (U-Th)/He in apatite in sottosuolo per la finestra dei Tauri occidentale. Le analisi sono state effettuate tramite un codice a elementi finiti, Pecube, risolvente l’equazione di diffusione del calore per una topografia evolvente nel tempo. Il codice è stato modificato per tener conto dei dati sotterranei. L’inversione dei dati è stata effettuata usando il Neighbourhood Algorithm (NA), per ottenere il più plausibile scenario di evoluzione morfotettonico. I risultati ottenuti per il Sempione mostrano: ipotetica evoluzione dello stile tettonico della faglia del Sempione da rolling hinge a low angle detachment a 6.5 Ma e la cessazione dell’attività a 3 Ma; costruzione del rilievo fino a 5.5 Ma, smantellamento da 5.5 Ma ad oggi, in coincidenza dei cambiamenti climatici Messiniani e relativi all’inizio delle maggiori glaciazioni; incremento dell’esumazione da 0–0.6 mm/anno a 0.6–1.2 mm/anno a 2.4 Ma nell’emisfero settentrionale. I risultati al Brennero mostrano: maggiore attività tettonica della faglia del Brennero (1.3 mm/anno), maggiore attività esumativa (1–2 mm/anno) prima dei 10 Ma; crollo dell’attività della faglia del Brennero fra 10 Ma e oggi (0.1 mm/anno) e dell’attività esumativa nello stesso periodo (0.1–0.3 mm/anno); nessun aumento del tasso esumativo o variazioni topografiche negli ultimi 5 Ma.

Behaviour and applications of elastic waves in structures and metamaterials

The present thesis focuses on elastic waves behaviour in ordinary structures as well as in acousto-elastic metamaterials via numerical and experimental applications. After a brief introduction on the behaviour of elastic guided waves in the framework of non-destructive evaluation (NDE) and structural health monitoring (SHM) and on the study of elastic waves propagation in acousto-elastic metamaterials, dispersion curves for thin-walled beams and arbitrary cross-section waveguides are extracted via Semi-Analytical Finite Element (SAFE) methods. Thus, a novel strategy tackling signal dispersion to locate defects in irregular waveguides is proposed and numerically validated. Finally, a time-reversal and laser-vibrometry based procedure for impact location is numerically and experimentally tested. In the second part, an introduction and a brief review of the basic definitions necessary to describe acousto-elastic metamaterials is provided. A numerical approach to extract dispersion properties in such structures is highlighted. Afterwards, solid-solid and solid-fluid phononic systems are discussed via numerical applications. In particular, band structures and transmission power spectra are predicted for 1P-2D, 2P-2D and 2P-3D phononic systems. In addition, attenuation bands in the ultrasonic as well as in the sonic frequency regimes are experimentally investigated. In the experimental validation, PZTs in a pitch-catch configuration and laser vibrometric measurements are performed on a PVC phononic plate in the ultrasonic frequency range and sound insulation index is computed for a 2P-3D phononic barrier in the sonic frequency range. In both cases the numerical-experimental results comparison confirms the existence of the numerical predicted band-gaps. Finally, the feasibility of an innovative passive isolation strategy based on giant elastic metamaterials is numerically proved to be practical for civil structures. In particular, attenuation of seismic waves is demonstrated via finite elements analyses. Further, a parametric study shows that depending on the soil properties, such an earthquake-proof barrier could lead to significant reduction of the superstructure displacement.

Modelling the dynamic response of rockfall protection barriers

Mountainous areas are prone to natural hazards like rockfalls. Among the many countermeasures, rockfall protection barriers represent an effective solution to mitigate the risk. They are metallic structures designed to intercept rocks falling from unstable slopes, thus dissipating the energy deriving from the impact. This study aims at providing a better understanding of the response of several rockfall barrier types, through the development of rather sophisticated three-dimensional numerical finite elements models which take into account for the highly dynamic and non-linear conditions of such events. The models are built considering the actual geometrical and mechanical properties of real systems. Particular attention is given to the connecting details between the structural components and to their interactions. The importance of the work lies in being able to support a wide experimental activity with appropriate numerical modelling. The data of several full-scale tests carried out on barrier prototypes, as well as on their structural components, are combined with results of numerical simulations. Though the models are designed with relatively simple solutions in order to obtain a low computational cost of the simulations, they are able to reproduce with great accuracy the test results, thus validating the reliability of the numerical strategy proposed for the design of these structures. The developed models have shown to be readily applied to predict the barrier performance under different possible scenarios, by varying the initial configuration of the structures and/or of the impact conditions. Furthermore, the numerical models enable to optimize the design of these structures and to evaluate the benefit of possible solutions. Finally it is shown they can be also used as a valuable supporting tool for the operators within a rockfall risk assessment procedure, to gain crucial understanding of the performance of existing barriers in working conditions.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

Rafael Moneo: la complessità del vuoto. Dalle matrici formali dell'opera di Oteiza e Chillida

The experience of void, essential to the production of forms and to make use them, can be considered as the base of the activities that attend to the formative processes. If void and matter constitutes the basic substances of architecture. Their role in the definition of form, the symbolic value and the constructive methods of it defines the quality of the space. This job inquires the character of space in the architecture of Moneo interpreting the meaning of the void in the Basque culture through the reading of the form matrices in the work of Jorge Oteiza and Eduardo Chillida. In the tie with the Basque culture a reading key is characterized by concurring to put in relation some of the theoretical principles expressed by Moneo on the relationship between place and time, in an unique and specific vision of the space. In the analysis of the process that determines the genesis of the architecture of Moneo emerges a trajectory whose direction is constructed on two pivos: on the one hand architecture like instrument of appropriation of the place, gushed from an acquaintance process who leans itself to the reading of the relations that define the place and of the resonances through which measuring it, on the other hand the architecture whose character is able to represent and to extend the time in which he is conceived, through the autonomy that is conferred to them from values. Following the trace characterized from this hypothesis, that is supported on the theories elaborated from Moneo, surveying deepens the reading of the principles that construct the sculptural work of Oteiza and Chillida, features from a search around the topic of the void and to its expression through the form. It is instrumental to the definition of a specific area that concurs to interpret the character of the space subtended to a vision of the place and the time, affine to the sensibility of Moneo and in some way not stranger to its cultural formation. The years of the academic formation, during which Moneo enters in contact with the Basque artistic culture, seem to be an important period in the birth of that knowledge that will leads him to the formulation of theories tied to the relationship between time, place and architecture. The values expressed through the experimental work of Oteiza and Chillida during years '50 are valid bases to the understanding of such relationships. In tracing a profile of the figures of Oteiza and Chillida, without the pretension that it is exhaustive for the reading of the complex historical period in which they are placed, but with the needs to put the work in a context, I want to be evidenced the important role carried out from the two artists from the Basque cultural area within which Moneo moves its first steps. The tie that approaches Moneo to the Basque culture following the personal trajectory of the formative experience interlaces to that one of important figures of the art and the Spanish architecture. One of the more meaningful relationships is born just during the years of his academic formation, from 1958 to the 1961, when he works like student in the professional office of the architect Francisco Sáenz de Oiza, who was teaching architectural design at the ETSAM. In these years many figures of Basque artists alternated at the professional office of Oiza that enjoys the important support of the manufacturer and maecenas Juan Huarte Beaumont, introduced to he from Oteiza. The tie between Huarte and Oteiza is solid and continuous in the years and it realizes in a contribution to many of the initiatives that makes of Oteiza a forwarder of the Basque culture. In the four years of collaboration with Oiza, Moneo has the opportunity to keep in contact with an atmosphere permeated by a constant search in the field of the plastic art and with figures directly connected to such atmosphere. It’s of a period of great intensity as in the production like in the promotion of the Basque art. The collective “Blanco y Negro”, than is held in 1959 at the Galería Darro to Madrid, is only one of the many times of an exhibition of the work of Oteiza and Chillida. The end of the Fifties is a period of international acknowledgment for Chillida that for Oteiza. The decade of the Fifties consecrates the hypotheses of a mythical past of the Basque people through the spread of the studies carried out in the antecedent years. The archaeological discoveries that join to a context already rich of signs of the prehistoric era, consolidate the knowledge of a strong cultural identity. Oteiza, like Chillida and other contemporary artists, believe in a cosmogonist conception belonging to the Basques, connected to their matriarchal mythological past. The void in its meaning of absence, in the Basque culture, thus as in various archaic and oriental religions, is equivalent to the spiritual fullness as essential condition to the revealing of essence. Retracing the archaic origins of the Basque culture emerges the deep meaning that the void assumes as key element in the religious interpretation of the passage from the life to the death. The symbology becomes rich of meaningful characters who derive from the fact that it is a chthonic cult. A representation of earth like place in which divine manifest itself but also like connection between divine and human, and this manipulation of the matter of which the earth it is composed is the tangible projection of the continuous search of the man towards God. The search of equilibrium between empty and full, that characterizes also the development of the form in architecture, in the Basque culture assumes therefore a peculiar value that returns like constant in great part of the plastic expressions, than in this context seem to be privileged regarding the other expressive forms. Oteiza and Chillida develop two original points of view in the representation of the void through the form. Both use of rigorous systems of rules sensitive to the physics principles and the characters of the matter. The last aim of the Oteiza’s construction is the void like limit of the knowledge, like border between known and unknown. It doesn’t means to reduce the sculptural object to an only allusive dimension because the void as physical and spiritual power is an active void, that possesses that value able to reveal the being through the trace of un-being. The void in its transcendental manifestation acts at the same time from universal and from particular, like in the atomic structure of the matter, in which on one side it constitutes the inner structure of every atom and on the other one it is necessary condition to the interaction between all the atoms. The void can be seen therefore as the action field that concurs the relations between the forms but is also the necessary condition to the same existence of the form. In the construction of Chillida the void represents that counterpart structuring the matter, inborn in it, the element in absence of which wouldn’t be variations neither distinctive characters to define the phenomenal variety of the world. The physics laws become the subject of the sculptural representation, the void are the instrument that concurs to catch up the equilibrium. Chillida dedicate himself to experience the space through the senses, to perceive of the qualities, to tell the physics laws which forge the matter in the form and the form arranges the places. From the artistic experience of the two sculptors they can be transposed, to the architectonic work of Moneo, those matrices on which they have constructed their original lyric expressions, where the void is absolute protagonist. An ambit is defined thus within which the matrices form them drafts from the work of Oteiza and Chillida can be traced in the definition of the process of birth and construction of the architecture of Moneo, but also in the relation that the architecture establishes with the place and in the time. The void becomes instrument to read the space constructed in its relationships that determine the proportions, rhythms, and relations. In this way the void concurs to interpret the architectonic space and to read the value of it, the quality of the spaces constructing it. This because it’s like an instrument of the composition, whose role is to maintain to the separation between the elements putting in evidence the field of relations. The void is that instrument that serves to characterize the elements that are with in the composition, related between each other, but distinguished. The meaning of the void therefore pushes the interpretation of the architectonic composition on the game of the relations between the elements that, independent and distinguished, strengthen themselves in their identity. On the one hand if void, as measurable reality, concurs all the dimensional changes quantifying the relationships between the parts, on the other hand its dialectic connotation concurs to search the equilibrium that regulated such variations. Equilibrium that therefore does not represent an obtained state applying criteria setting up from arbitrary rules but that depends from the intimate nature of the matter and its embodiment in the form. The production of a form, or a formal system that can be finalized to the construction of a building, is indissolubly tied to the technique that is based on the acquaintance of the formal vocation of the matter, and what it also can representing, meaning, expresses itself in characterizing the site. For Moneo, in fact, the space defined from the architecture is above all a site, because the essence of the site is based on the construction. When Moneo speaks about “birth of the idea of plan” like essential moment in the construction process of the architecture, it refers to a process whose complexity cannot be born other than from a deepened acquaintance of the site that leads to the comprehension of its specificity. Specificity arise from the infinite sum of relations, than for Moneo is the story of the oneness of a site, of its history, of the cultural identity and of the dimensional characters that that they are tied to it beyond that to the physical characteristics of the site. This vision is leaned to a solid made physical structure of perceptions, of distances, guideline and references that then make that the process is first of all acquaintance, appropriation. Appropriation that however does not happen for directed consequence because does not exist a relationship of cause and effect between place and architecture, thus as an univocal and exclusive way does not exist to arrive to a representation of an idea. An approach that, through the construction of the place where the architecture acquires its being, searches an expression of its sense of the truth. The proposal of a distinction for areas like space, matter, spirit and time, answering to the issues that scan the topics of the planning search of Moneo, concurs a more immediate reading of the systems subtended to the composition principles, through which is related the recurrent architectonic elements in its planning dictionary. From the dialectic between the opposites that is expressed in the duality of the form, through the definition of a complex element that can mediate between inside and outside as a real system of exchange, Moneo experiences the form development of the building deepening the relations that the volume establishes in the site. From time to time the invention of a system used to answer to the needs of the program and to resolve the dual character of the construction in an only gesture, involves a deep acquaintance of the professional practice. The technical aspect is the essential support to which the construction of the system is indissolubly tied. What therefore arouses interest is the search of the criteria and the way to construct that can reveal essential aspects of the being of the things. The constructive process demands, in fact, the acquaintance of the formative properties of the matter. Property from which the reflections gush on the relations that can be born around the architecture through the resonance produced from the forms. The void, in fact, through the form is in a position to constructing the site establishing a reciprocity relation. A reciprocity that is determined in the game between empty and full and of the forms between each other, regarding around, but also with regard to the subjective experience. The construction of a background used to amplify what is arranged on it and to clearly show the relations between the parts and at the same time able to tie itself with around opening the space of the vision, is a system that in the architecture of Moneo has one of its more effective applications in the use of the platform used like architectonic element. The spiritual force of this architectonic gesture is in the ability to define a place whose projecting intention is perceived and shared with who experience and has lived like some instrument to contact the cosmic forces, in a delicate process that lead to the equilibrium with them, but in completely physical way. The principles subtended to the construction of the form taken from the study of the void and the relations that it concurs, lead to express human values in the construction of the site. The validity of these principles however is tested from the time. The time is what Moneo considers as filter that every architecture is subordinate to and the survival of architecture, or any of its formal characters, reveals them the validity of the principles that have determined it. It manifests thus, in the tie between the spatial and spiritual dimension, between the material and the worldly dimension, the state of necessity that leads, in the construction of the architecture, to establish a contact with the forces of the universe and the intimate world, through a process that translate that necessity in elaboration of a formal system.

Study of the thermal effects induced by magnetic resonance on endocardial leads: numerical models and experimental validation

Magnetic resonance imaging (MRI) is today precluded to patients bearing active implantable medical devices AIMDs). The great advantages related to this diagnostic modality, together with the increasing number of people benefiting from implantable devices, in particular pacemakers(PM)and carioverter/defibrillators (ICD), is prompting the scientific community the study the possibility to extend MRI also to implanted patients. The MRI induced specific absorption rate (SAR) and the consequent heating of biological tissues is one of the major concerns that makes patients bearing metallic structures contraindicated for MRI scans. To date, both in-vivo and in-vitro studies have demonstrated the potentially dangerous temperature increase caused by the radiofrequency (RF) field generated during MRI procedures in the tissues surrounding thin metallic implants. On the other side, the technical evolution of MRI scanners and of AIMDs together with published data on the lack of adverse events have reopened the interest in this field and suggest that, under given conditions, MRI can be safely performed also in implanted patients. With a better understanding of the hazards of performing MRI scans on implanted patients as well as the development of MRI safe devices, we may soon enter an era where the ability of this imaging modality may be more widely used to assist in the appropriate diagnosis of patients with devices. In this study both experimental measures and numerical analysis were performed. Aim of the study is to systematically investigate the effects of the MRI RF filed on implantable devices and to identify the elements that play a major role in the induced heating. Furthermore, we aimed at developing a realistic numerical model able to simulate the interactions between an RF coil for MRI and biological tissues implanted with a PM, and to predict the induced SAR as a function of the particular path of the PM lead. The methods developed and validated during the PhD program led to the design of an experimental framework for the accurate measure of PM lead heating induced by MRI systems. In addition, numerical models based on Finite-Differences Time-Domain (FDTD) simulations were validated to obtain a general tool for investigating the large number of parameters and factors involved in this complex phenomenon. The results obtained demonstrated that the MRI induced heating on metallic implants is a real risk that represents a contraindication in extending MRI scans also to patient bearing a PM, an ICD, or other thin metallic objects. On the other side, both experimental data and numerical results show that, under particular conditions, MRI procedures might be consider reasonably safe also for an implanted patient. The complexity and the large number of variables involved, make difficult to define a unique set of such conditions: when the benefits of a MRI investigation cannot be obtained using other imaging techniques, the possibility to perform the scan should not be immediately excluded, but some considerations are always needed.

Stability, viscous dissipation and local thermal non-equilibrium in fluid saturated porous media

In this thesis, the field of study related to the stability analysis of fluid saturated porous media is investigated. In particular the contribution of the viscous heating to the onset of convective instability in the flow through ducts is analysed. In order to evaluate the contribution of the viscous dissipation, different geometries, different models describing the balance equations and different boundary conditions are used. Moreover, the local thermal non-equilibrium model is used to study the evolution of the temperature differences between the fluid and the solid matrix in a thermal boundary layer problem. On studying the onset of instability, different techniques for eigenvalue problems has been used. Analytical solutions, asymptotic analyses and numerical solutions by means of original and commercial codes are carried out.

Hybrid Methods for Resource Allocation and Scheduling Problems in Deterministic and Stochastic Environments

This work presents exact, hybrid algorithms for mixed resource Allocation and Scheduling problems; in general terms, those consist into assigning over time finite capacity resources to a set of precedence connected activities. The proposed methods have broad applicability, but are mainly motivated by applications in the field of Embedded System Design. In particular, high-performance embedded computing recently witnessed the shift from single CPU platforms with application-specific accelerators to programmable Multi Processor Systems-on-Chip (MPSoCs). Those allow higher flexibility, real time performance and low energy consumption, but the programmer must be able to effectively exploit the platform parallelism. This raises interest in the development of algorithmic techniques to be embedded in CAD tools; in particular, given a specific application and platform, the objective if to perform optimal allocation of hardware resources and to compute an execution schedule. On this regard, since embedded systems tend to run the same set of applications for their entire lifetime, off-line, exact optimization approaches are particularly appealing. Quite surprisingly, the use of exact algorithms has not been well investigated so far; this is in part motivated by the complexity of integrated allocation and scheduling, setting tough challenges for ``pure'' combinatorial methods. The use of hybrid CP/OR approaches presents the opportunity to exploit mutual advantages of different methods, while compensating for their weaknesses. In this work, we consider in first instance an Allocation and Scheduling problem over the Cell BE processor by Sony, IBM and Toshiba; we propose three different solution methods, leveraging decomposition, cut generation and heuristic guided search. Next, we face Allocation and Scheduling of so-called Conditional Task Graphs, explicitly accounting for branches with outcome not known at design time; we extend the CP scheduling framework to effectively deal with the introduced stochastic elements. Finally, we address Allocation and Scheduling with uncertain, bounded execution times, via conflict based tree search; we introduce a simple and flexible time model to take into account duration variability and provide an efficient conflict detection method. The proposed approaches achieve good results on practical size problem, thus demonstrating the use of exact approaches for system design is feasible. Furthermore, the developed techniques bring significant contributions to combinatorial optimization methods.

Microscopic Modeling on complex networks

The field of complex systems is a growing body of knowledge, It can be applied to countless different topics, from physics to computer science, biology, information theory and sociology. The main focus of this work is the use of microscopic models to study the behavior of urban mobility, which characteristics make it a paradigmatic example of complexity. In particular, simulations are used to investigate phase changes in a finite size open Manhattan-like urban road network under different traffic conditions, in search for the parameters to identify phase transitions, equilibrium and non-equilibrium conditions . It is shown how the flow-density macroscopic fundamental diagram of the simulation shows,like real traffic, hysteresis behavior in the transition from the congested phase to the free flow phase, and how the different regimes can be identified studying the statistics of road occupancy.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

Quantum Correlations in Field Theory and Integrable Systems

In this thesis we will investigate some properties of one-dimensional quantum systems. From a theoretical point of view quantum models in one dimension are particularly interesting because they are strongly interacting, since particles cannot avoid each other in their motion, and you we can never ignore collisions. Yet, integrable models often generate new and non-trivial solutions, which could not be found perturbatively. In this dissertation we shall focus on two important aspects of integrable one- dimensional models: Their entanglement properties at equilibrium and their dynamical correlators after a quantum quench. The first part of the thesis will be therefore devoted to the study of the entanglement entropy in one- dimensional integrable systems, with a special focus on the XYZ spin-1/2 chain, which, in addition to being integrable, is also an interacting model. We will derive its Renyi entropies in the thermodynamic limit and its behaviour in different phases and for different values of the mass-gap will be analysed. In the second part of the thesis we will instead study the dynamics of correlators after a quantum quench , which represent a powerful tool to measure how perturbations and signals propagate through a quantum chain. The emphasis will be on the Transverse Field Ising Chain and the O(3) non-linear sigma model, which will be both studied by means of a semi-classical approach. Moreover in the last chapter we will demonstrate a general result about the dynamics of correlation functions of local observables after a quantum quench in integrable systems. In particular we will show that if there are not long-range interactions in the final Hamiltonian, then the dynamics of the model (non equal- time correlations) is described by the same statistical ensemble that describes its statical properties (equal-time correlations).