Recent experimental results reveal that with rising concentration of surfactants, the generated bubble‐induced turbulence (BIT) increases as well, although the bubbles rise slower. Herein, it is shown that this phenomenon has not be considered so far in the standard Reynolds‐averaged Navier‐Stokes (RANS) modeling for BIT, using the Ansatz with source terms added to the transport equation of turbulent kinetic energy and dissipation, respectively. Some problematic issues related to the common concepts for closure of the liquid‐phase turbulence kinetic energy equation are also discussed. Furthermore, this issue is demonstrated by simulating the experiments using two‐fluid approach.