Monte Carlo calculations of neutron thermalization in a heterogeneous system