Simulating neutron propagations with FLUKA, GEANT4 and MCNP