The ever increasing power rates of modern IC engines pose a great challenge for maintaining the optimal temperature conditions. The current trend in the thermal management of IC engines is to utilize the limited and controlled nucleate boiling to achieve higher heat transfer rates. The modern CFD instrumentation offers a wide range of methods for predicting nucleate and film boiling. Of the available models the mechanistic RPI wall boiling model is selected in present paper in the cojuntion with the Eulerian multiphase model available in ANSYS Fluent solver to simulate subcooled nucleate boiling in engine cooling jacket.Results obtained from the simulation of boiling flow in cooling passages is then compared with the experimental data from author’s previous work and show a decent agreement.