The subcooled boiling is an effective heat transfer method by using latent coolant heat in modern engine cooling passages. Thus, the cylinder head temperature of internal combustion (IC) engine could be maintained below the damaged temperature through controlled nuclear boiling. Hence，in this paper, The Eulerian multiphase and RPI subcooled boiling models coupled with Hibiki, Kocamustafaogullari-Ishiii and Kurul-Podowski sub-models were employed to study the heat and mass transfer in cylinder head cooling passages. Firstly, the two-phase RPI boiling model was validated by Robinson’s experimental data with (50% ethylene glycol aqueous) solution of 50% volume fraction as working coolant. The benchmark simulated results were in good agreement with coincide well with experimental data. Then, the two-phase RPI boiling model could be applied to simulate subcooled flow boiling heat transfer phenomenon in practical engine cooling passages. Furthermore, in order to investigate the cylinder head temperature field was measured with 20 points. Compared with the experiment results, the temperature error between the numerical simulation results and measure is within 5 %.