A mathematical model was written for simulating the removal of phenol from wastewater in enzyme-loaded membrane reactor (EMR). The numerical simulation program was developed so as to predict the degradation of phenol through an EMR. Numerical model proves to be effective in searching for optimal operating conditions and creating an optimal microenvironment for the biocatalyst in order to optimize productivity. In this study, several dimensionless parameters such as Thiele Modulus (Φ2, dimensionless Michaelis-Menten constant (ξ), Peclet number (Pe) were introduced to simplify their effects on system efficiency. In particular, the study of phenol conversion at different feed compositions shows that low phenol concentrations and high Thiele Modulus values lead to higher reactant degradation.