A mathematical model that describes the heterogeneous reaction-diffusion process involved in penicillin G hydrolysis in a batch reactor with immobilized penicillin G acylase is presented. The reaction system includes the bulk liquid phase containing the dissolved substrate (and products) and the solid biocatalyst phase represented by glyoxyl-agarose spherical porous particles carrying the enzyme. The equations consider reaction and diffusion components that are presented in dimensionless form. This is a complex reaction system in which both products of reaction and the substrate itself are inhibitors. The simulation of a batch reactor performance with immobilized penicillin G acylase is presented and discussed for the internal diffusional restrictions impact on effectiveness and productivity. Increasing internal diffusional restrictions, through increasing catalyst particle size and enzyme loading, causes impaired catalyst efficiency expressed in a reduction of effectiveness factor and specific productivity. High penicillin G initial concentrations decrease the impact of internal diffusional restrictions by increasing the mass transfer towards porous catalyst until product inhibition becomes significant over approximately 50 mM of initial penicillin G, where a drop in conversion rate and a maximum in specific productivity are then obtained. Results highlight the relevance of considering internal diffusional restrictions, reactor performance, and productivity analysis for proper catalyst and reactor design.