The dimensionless kinetic energy dissipation rate C_epsilon is estimated from numerical simulations of statistically stationary isotropic box turbulence that is slightly compressible. The Taylor microscale Reynolds number Re_lambda range is 20 < Re_lambda < 220 and the statistical stationarity is achieved with a random phase forcing method. The strong Re_lambda dependence of C_epsilon abates when Re_lambda approx. 100 after which C_epsilon slowly approaches approx 0.5 a value slightly different to previously reported simulations but in good agreement with experimental results. If C_epsilon is estimated at a specific time step from the time series of the quantities involved it is necessary to account for the time lag between energy injection and energy dissipation. Also, the resulting value can differ from the ensemble averaged value by up to +-30%. This may explain the spread in results from previously published estimates of C_epsilon.