A dynamic Euler-Euler model for gas-liquid flows subject to mass transfer and simultaneous chemical reaction is assembled. The bubble size distribution is computed from an equation governing the evolution of average bubble mass. High-performance finite element tools for numerical solution of the problem at hand are discussed. Computational results for the chemisorption of CO_2 into NaOH shed some light on the implications of reaction enhanced mass transfer in bubble columns and call for further research in this field.