Characteristics of two-phase flow and heat transfer were numerically investigated in a submerged gas injection system when temperature of the injected gas was different from that of the liquid. The Eulerian approach was used for both the continuous and dispersed phases. The turbulence in the liquid phase was modeled using the standard $k-varepsilon$$varepsilon$ turbulence model. The interphase friction and heat transfer coefficient were calculated from the correlations available in the literature. The turbulent dispersion of the phases was modeled by a "dispersion Prandtl number". In the case with heat transfer where the temperature of the injected gas is higher than the mean liquid temperature, the axial and the radial velocities are lower in comparison with the case of homogeneous temperatures. The results in the present research are of interest in the design and operation of a wide variety of material and chemical processes.