Numerical model for cooling a gravitational wave detector below 1-K