Particle laden flows with reactive particles are common in industrial applications. Chemical reactions inside the particle can generate a Stefan flow that affects heat, mass and momentum transfer between the particle and the bulk flow. This study aims at investigating the effect of Stefan flow on the drag coefficient of a spherical particle immersed in a uniform flow under isothermal conditions. Fully resolved simulations were carried out for particle Reynolds numbers ranging from 0.2 to 14 and Stefan flow Reynolds numbers from (-1) to 3, using the immersed boundary method for treating fluid-solid interactions. Results showed that the drag coefficient decreased with an increase of the outward Stefan flow. The main reason was the change in viscous force by the expansion of the boundary layer surrounding the particle. A simple model was developed based on this physical interpretation. With only one fitting parameter, the performance of the model to describe the simulation data were comparable to previous empirical models.