Fractional Brownian motion is a widely used stochastic process that is particularly suited to model anomalous diffusion. We focus on capturing the mean and variance of fractional Brownian motion reflected at level 0. As explicit expressions or numerical techniques are not available, we base our analysis on Monte Carlo simulation. Our main findings concern closed-form approximations of the mean and variance, with a near-perfect fit.