Numerical modelling of rock masses interacting with dams: the case of Ridracoli

A three-dimensional Direct Finite Element procedure is here presented which takes into account most of the factors affecting the interaction problem of the dam-water-foundation system, whilst keeping the computational cost at a reasonable level by introducing some simplified hypotheses. A truncated domain is defined, and the dynamic behaviour of the system is treated as a wave-scattering problem where the presence of the dam perturbs an original free-field system. The rock foundation truncated boundaries are enclosed by a set of free-field one-dimensional and two-dimensional systems which transmit the effective forces to the main model and apply adsorbing viscous boundaries to ensure radiation damping. The water domain is treated as an added mass moving with the dam. A strategy is proposed to keep the viscous dampers at the boundaries unloaded during the initial phases of analysis, when the static loads are initialised, and thus avoid spurious displacements. A focus is given to the nonlinear behaviour of the rock foundation, with concentrated plasticity along the natural discontinuities of the rock mass, immersed in an otherwise linear elastic medium with Rayleigh damping. The entire procedure is implemented in the commercial software Abaqus®, whose base code is enriched with specific user subroutines when needed. All the extra coding is attached to the Thesis and tested against analytical results and simple examples. Possible rock wedge instabilities induced by intense ground motion, which are not easily investigated within a comprehensive model of the dam-water-foundation system, are treated separately with a simplified decoupled dynamic approach derived from the classical Newmark method, integrated with FE calculation of dam thrust on the wedges during the earthquake. Both the described approaches are applied to the case study of the Ridracoli arch-gravity dam (Italy) in order to investigate its seismic response to the Maximum Credible Earthquake (MCE) in a full reservoir condition.

A comparison between vertical mixing parameterizations in the Levantine Sea

The goal of this thesis was the study of an optimal vertical mixing parameterization scheme in a mesoscale dominated field characterized from a strong vorticity and the presence of a layer of colder, less saline water at about 100 m depth (Atlantic Waters); in these conditions we compared six different experiments, that differ by the turbulent closure schemes, the presence or not of an enhanced diffusion parameterization and the presence or not of a double diffusion mixing parameterization. To evaluate the performance of the experiments and the model we compared the simulations with the ARGO observations of temperature and salinity available in our domain, in our period of interest. The conclusions were the following: • the increase of the resolution gives better results in terms of temperature in all the considered cases, and in terms of salinity. • The comparisons between the Pacanovski-Philander and the TKE turbulent closure schemes don’t show significant differences when the simulations are compared to the observations. • The removing of the enhanced diffusion parameterization in presence of the TKE turbulent closure submodel doesn’t give positive results, and show limitations in the resolving of gravitational instabilities near the surface • The k-ϵ turbulent closure model utilized in all the GLS experiments, is the best performing closure model among the three considered, with positive results in all the salinity comparison with the in situ observation and in most of the temperature comparisons. • The double mixing parameterization utilized in the k-ϵ closure submodel improves the results of the experiments improving both the temperature and salinity in comparison with the ARGO data.