In the present work, a numerical approach is proposed to simulate two phase pool boiling on a staggered tube bundle. The simulations are carried out using the Eulerian-Eulerian framework for two phase flows, accounting for the nucleate pool boiling on the surface of the tube bundle. A favourable agreement was reported between the numerical results and previously published experimental data for pool boiling heat transfer. The numerical results are presented for four working fluids, including water, ethanol, R134a and R22, each at three different saturation temperatures. The effect of saturation temperature, the variation in the heat transfer within the tube bundle and the effect of the working fluid are fully addressed. (C) 2018 Elsevier Ltd. All rights reserved.