An Eulerian–Lagrangian method for the simulation of the oxygen concentration dissolved by a two-phase turbulent jet system