In this thesis we present the results of an extensive campaign of direct numerical simulations of Rayleigh-Bénard convection at high Prandtl numbers (10 −1 ≤ P r ≤ 10 4 ) and moderate Rayleigh numbers (10 5 ≤ P r ≤ 10 9 ). The computational domain is a cylindrical cell of aspect-ratio (diameter over cell height) Γ = 1/2, with the no-slip condition imposed to the boundaries.By scaling the results, we find a 1/ √ P r correction to apply to the free-fall velocity, obtaining a more appropriate representation of the large scale velocity at high P r. We investigate the Nusselt and the Reynolds number dependence on Ra and P r, comparing the results to previous numerical and experimental work. At high P r the scaling behavior of the Nusselt number with respect to Ra is generally consistent with the power-law exponent 0.309.The Nusselt number is independent of P r, even at the highest Ra simulated. The Reynolds number scales as Re ∼ √ Ra/P r, neglecting logarithmic corrections. We analyze the global and local features of viscous and thermal boundary layers and their scaling behavior with respect to Rayleigh and Prandtl numbers, and with respect to Reynolds and Peclet numbers.We find that the flow approaches a saturation regime when Reynolds number decreases below the critical value Re s ≃ 40. The thermal boundary layer thickness turns out to increase slightly even when the Peclet number increases. We explain this behavior as a combined effect of the Peclet number and the viscous boundary layer influences.The range of Ra and P r simulated contains steady, periodic and turbulent solutions.A rough estimate of the transition from steady to unsteady flow is obtained by monitoring the time-evolution of the system until it reaches stationary solutions (Ra U ≃ 7.5 × 10 6 at P r = 10 3 ). We find multiple solutions as long-term phenomena at Ra = 10 8 and P r = 10 3 which, however, do not result in significantly different Nusselt number. One of these multiple solutions, even if stable for a long time interval, shows a break in the mid-plane symmetry of the temperature profile. The result is similar to that of some non-Boussinesq effects. We analyze the flow structures through the transitional phases by direct visualizations of the temperature and velocity fields. We also describe how the behavior of the flow structures changes for increasing P r. A wide variety of large-scale circulations and plumes structures are found. The single-roll circulation is characteristic only of the steady and periodic solutions. For other solutions, at lower P r, the mean flow generally consists of two opposite toroidal structures; at higher P r, the flow is organized in multi-cell structures extending mostly in the vertical direction. At high P r, plumes detach from sheet-like structures. The different large-scale-structure signatures are generally reflected in the data trends with rei spect to Ra, but not in those with respect to P r. In particular, the Nusselt number is independent of P r, even when the flow structures appear strongly diffe...