In this paper, a novel simulation methodology based on the reproducing kernels is proposed for solving the fractional order integro‐differential transport model for a nuclear reactor. The analysis carried out in this paper thus forms a crucial step in the process of development of fractional calculus as well as nuclear science models. The fractional derivative is described in the Captuo Riemann–Liouville sense. Results are presented graphically and in tabulated forms to study the efficiency and accuracy of method. The present scheme is very simple, effective, and appropriate for obtaining numerical simulation of nuclear science models. Copyright © 2017 John Wiley & Sons, Ltd.