In the present work, a new simulation of nanofluid/vapor two-phase flow inside the 2-D rectangular boiling chamber was numerically investigated. The Eulerian–Eulerian approach used to predict the boiling curve and the interaction between two phases. The surface modification during pool boiling of silica nanofluid represented by surface roughness and wettability is put into the account in this simulation. New closure correlations regarding the nucleation sites density and bubble departure diameter during boiling of silica nanofluid were inserted to extend the boiling model in this work. Besides, the bubble waiting time coefficient which involved in quenching heat flux under heat flux partitioning HFP model was corrected to improve the results of this study. The numerical results validated with experimental works in the literature, and they revealed good agreements for both pure water and nanofluids. The results found that when improving the heat flux partitioning model HFP by considering the surface modification of nucleate pool boiling parameters, it will give more mechanistic sights compared to the classical model, which is used for predicting of boiling heat transfer of pure liquid.