This paper presents a VOF-based approach for the numerical simulation of non-reactive or reactive mass transfer from a rising bubble into the ambient liquid. The mathematical model and the numerical method allow for different solubilities of the species in the respective fluid phases, while volume changes due to mass transfer are neglected. The discontinuity of the transfer species concentration at the interface is modeled by means of Henry’s law. Simulations are carried out with an extended version of the highly parallelized code FS3D, which employs an advanced Volume-Of-Fluid (VOF) method. To analyze complex reaction schemes a local selectivity is introduced. First numerical results demonstrate the potential of the method for contributing to process intensification of gas-liquid reactors.