Mixing by gas injection is an operation used in industrial processes such as wastewater treatment, metallurgy, or methanization in which pressurized gas is injected into a fluid in order to reduce concentrations and temperatures gradients. This study demonstrates how the CFD toolbox OpenFOAM can be used to simulate such flows.Experimental measurements and observations have been performed on a pilot-scale reactor where pressurized air is injected in a yield stress fluid. The volume of fluid method and an adaptive mesh with refinement at the interface have been used to track the gas inclusions. The numerical model accuracy has been assessed by comparing experimental and numerical results related to the bubble's frequency, dimensions, and rising velocities as well as the fluid recirculation, yielded, and unyielded regions in the tank. The influence of injection parameters such as the injection flow rate and the fluid rheological parameters has been quantified.