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 (10 3 -10 5 ). K E Y W O R D Sboundary layer, bubble columns, front tracking, mass transfer, subgrid scale modeling

The mass transfer between a rising bubble and the surrounding liquid is mainly determined by an extremely thin layer of dissolved gas near the bubble interface. Resolving this concentration boundary layer in numerical simulations is computationally expensive and limited to low Péclet numbers. Subgrid-scale (SGS) models mitigate the resolution

The mass transfer between a rising bubble and the surrounding liquid is mainly determined by an extremely thin layer of dissolved gas forming at the liquid side of the gas-liquid interface. Resolving this concentration boundary layer in numerical simulations is computationally expensive. Subgrid-scale models mitigate the resolution requirements enormously and allow approximating the mass transfer in industrially relevant flow conditions with high accuracy. However, the development and validation of such models is difficult as only integral mass transfer data for steady-state conditions are available. Therefore, it is difficult to assess the validity of the sub-grid models in transient conditions. In this contribution, we compare the local and global mass transfer of an improved subgrid-scale model for rising bubbles (Re = 72-569 and Sc = 10^2-10^4) to a single-phase simulation approach, which maps the two-phase flow field to a highly-resolved mesh comprising only the liquid phase.