A numerical study of turbulent condensing vapor jet submerged in subcooled liquids has been conducted. A physical model of the process is presented employing the locally homogeneous flow approximation of two phase flow in conjunction with a $kappa$-$varepsilon$-g model of turbulence properties. In this model the turbulence is represented by differential equations for its kinetic energy and dissipation. A differential equation for the concentration fluctuations is solved and a clipped normal probability distribution function is proposed for the mixture fraction. Effects of steam mass flux, pool temperature and nozzle internal diameter on the condensing vapor jet are also analyzed. The model is evaluated using existing data for turbulent condensing vapor jets. The agreement between the predictions and the available experimental data is good.