  We consider the Navier-Stokes-Fourier system governing the motion of a
general compressible, heat conducting, Newtonian fluid driven by random
initial/boundary data. Convergence of the stochastic collocation and Monte
Carlo numerical methods is shown under the hypothesis that approximate
solutions are bounded in probability. Abstract results are illustrated by
numerical experiments for the Rayleigh-Benard convection problem.
