The results of numerical modeling of hydrodynamics and heat and mass transfer at boiling of subcooled liquids in conditions of forced flow in vertical heated pipes are presented. The mathematical model is based on the use of Euler’s description of mass, motion and energy conservation for liquid and gas phases, recorded within the framework of the theory of interacting continua. The turbulent characteristics of the fluid are calculated using a modified model of transfer of components of the Reynolds stress tensor, taking into account the presence of the gaseous phase in the medium. For an approximate calculation of the heat transfer coefficient for bubble boiling of a liquid near the heat-generating wall, generalized empirical dependences are used, taking into account the various mechanisms of heat transfer in a two-phase vapor-liquid medium. Comparison of the results of numerical modeling with experimental data has shown that the proposed approach allows to simulate bubble boiling modes in a wide range of pressure values, mass flow rates, heating modes of subcooled liquids during forced turbulent fluid flow in vertical pipes.