The purpose of this study was to develop a numerical model to estimate the oxygen-transfer rate for a laboratory-scale bottom aeration system at a 1.28 L reactor volume and to contribute to fundamental knowledge regarding the oxygenation of surfactant solutions. The primary goal of the study has been to develop a computational fluid dynamics (CFD) model using Euler–Euler (EE) mixture model coupled with the advection-diffusion equation to predict the oxygen-transfer rate in bubble columns containing clean water. The secondary goal has been to apply the model to water-based solutions containing the surfactant lauric acid (DDA) to identify options for further development of the model to make it applicable to surfactant solution systems. The Sauter mean diameter (SMD) was calculated to represent the average bubble diameter, based on available experimental data for different combinations of superficial velocities rate and DDA concentration. The oxygen-transfer rate in clean water fit well with experimental data at lower superficial velocities, and the differences in volumetric mass-transfer coefficients were 0.7% and 3.3% for superficial velocities of 0.24 cm/s and 0.48 cm/s, respectively. For surfactant solutions, the model overestimates the oxygen-transfer rate due to surfactant adsorption at the bubble/water interface and the consequent decrease in the mass-transfer coefficient not being modeled sufficiently. A correction factor for the mass-transfer coefficient based on a larger sample size of experimental data may need to be calculated and applied to improve model predictability.