Two-fluid model is a common method to simulate the subcooled flow boiling heat transfer, in which the wall boiling model is mainly used for the partition of wall heat flux and the mass transfer between two phases on the wall. The model determines the amount of vapor phase and predicts the cross-sectional void fraction in the channel, nucleate site density and bubble departure diameter play an important role in the accurate prediction of wall boiling model. Eulerian two-fluid model coupled with Rensselaer Polytechnic Institute (RPI) wall boiling model is employed to simulate the heat transfer characteristics and boiling phenomena in vertical narrow rectangular channels by using FLUENT code. Based on the experimental data of subcooled boiling in vertical narrow rectangular channel, different combinations of nucleate site density and bubble departure diameter correlations are used to calculate under different conditions of heat flux and inlet subcooling. Comparing the calculated heat transfer coefficients along the vertical height with experimental results, it can be found that these two parameters have a significant effect on the subcooled boiling heat transfer in narrow rectangular channels. Different parameter combinations lead to differences in wall heat flux distribution, different heat flux and inlet subcooling also have different effects on these models, which eventually lead to different evaporative heat flux, thus affecting the prediction of void fraction.