The present paper proposes an application for an Euler–Lagrangian one-way coupling method for CO2 in water bubble tower reactor design. As a first step to model a large-scale bubble reactor a single bubble while rising and dissolving can be modeled as a spherical body submitted to punctual forces. Combined with a description for the decreasing mass and size of a dissolving bubble in an aqueous solution this will be basis for bubble swarm modeling. A link between mass transfer boundary layer and bubble force is established through the mass time derivative of the dissolving bubble as function of the local Reynolds, Schmidt and Sherwood numbers. Furthermore, the turbulent mixer fluid environment of a stirred reactor is presented as the basis of this research work for a slice of fluid and a rotating mesh model is presented as the basis for the simulation. A turbulent flow interface (k–ɛ turbulence) from a commercial software is used. Finally, bubble swarm displacement is introduced and bubbles are tracked individually using a Lagrangian spherical particle approach. Aqueous CO2 concentration as function of time and position is studied. Results are of relevance for CO2 mineralization processes, specifically for the bubble reactor used in the so-called Slag2PCC process being scaled up in Finland.