A numerical model for predicting diffusivity in a biocatalyst particle with heterogeneous reactions