Asymmetric drift and halo flattening in disk galaxies

The main result in this work is the solution of the Jeans equations for an axisymmetric galaxy model containing a baryonic component (distributed according to a Miyamoto-Nagai profile) and a dark matter halo (described by the Binney logarithmic potential). The velocity dispersion, azimuthal velocity and some other interesting quantities such as the asymmetric drift are studied, along with the influence of the model parameters on these (observable) quantities. We also give an estimate for the velocity of the radial flow, caused by the asymmetric drift. Other than the mathematical beauty that lies in solving a model analytically, the interest of this kind of results can be mainly found in numerical simulations that study the evolution of gas flows. For example, it is important to know how certain parameters such as the shape (oblate, prolate, spherical) of a dark matter halo, or the flattening of the baryonic matter, or the mass ratio between dark and baryonic matter, have an influence on observable quantities such as the velocity dispersion. In the introductory chapter, we discuss the Jeans equations, which provide information about the velocity dispersion of a system. Next we will consider some dynamical quantities that will be useful in the rest of the work, e.g. the asymmetric drift. In Chapter 2 we discuss in some more detail the family of galaxy models we studied. In Chapter 3 we give the solution of the Jeans equations. Chapter 4 describes and illustrates the behaviour of the velocity dispersion, as a function of the several parameters, along with asymptotic expansions. In Chapter 5 we will investigate the behaviour of certain dynamical quantities for this model. We conclude with a discussion in Chapter 6.

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.

Regularization methods for the solution of a nonlinear least-squares problem in tomography

In this work we study a polyenergetic and multimaterial model for the breast image reconstruction in Digital Tomosynthesis, taking into consideration the variety of the materials forming the object and the polyenergetic nature of the X-rays beam. The modelling of the problem leads to the resolution of a high-dimensional nonlinear least-squares problem that, due to its nature of inverse ill-posed problem, needs some kind of regularization. We test two main classes of methods: the Levenberg-Marquardt method (together with the Conjugate Gradient method for the computation of the descent direction) and two limited-memory BFGS-like methods (L-BFGS). We perform some experiments for different values of the regularization parameter (constant or varying at each iteration), tolerances and stop conditions. Finally, we analyse the performance of the several methods comparing relative errors, iterations number, times and the qualities of the reconstructed images.