  In the present paper we consider the initial data, external force, viscosity
coefficients, and heat conductivity coefficient as random data for the
compressible Navier--Stokes--Fourier system. The Monte Carlo method, which is
frequently used for the approximation of statistical moments, is combined with
a suitable deterministic discretisation method in physical space and time.
Under the assumption that numerical densities and temperatures are bounded in
probability, we prove the convergence of random finite volume solutions to a
statistical strong solution by applying genuine stochastic compactness
arguments. Further, we show the convergence and error estimates for the Monte
Carlo estimators of the expectation and deviation. We present several numerical
results to illustrate the theoretical results.
