Monte Carlo simulation of thermal neutron flux of americium–beryllium source used in neutron activation analysis