Gas–liquid bubble column reactors are often used in industry because of their favorable mass transfer characteristics. The bubble mass boundary layer in these systems is generally one order of magnitude thinner than the momentum boundary. To resolve it in simulations, a subgrid scale model will account for the sharp concentration variation in the vicinity of the interface. In this work, the subgrid scale model of Aboulhasanzadeh et al., Chem Eng Sci, 2012, 75:456–467 embedded in our in‐house front tracking framework, has been improved to prevent numerical mass transfer due to remeshing operations. Furthermore, two different approximations of the mass distribution in the boundary layer have been tested. The local and global predicted Sherwood number has been verified for mass transfer from bubbles in the creeping and potential flow regimes. In addition, the correct Sherwood number has been predicted for free rising bubbles at several Eötvös and Morton numbers with industrial relevant Schmidt numbers (103–105).