Boiling flows are widely encountered in several engineering and industrial processes. They have a special interest in nuclear industry, where a Computational Fluid Dynamic (CFD) thermohydraulic investigation becomes very popular for design and safety. Many attempts to model numerically subcooled nucleate boiling flows can be found in the literature, where several interfacial forces acting on bubbles which are interacting on the bulk fluid were neglected, due to the hard convergence of the calculations, or to the bad accuracy of the obtained results. In this paper, a sensitivity analysis is carried out for the interfacial forces acting on bubbles during subcooled nucleate boiling flows. For this purpose, 2D CFD axisymmetric simulations based on an Eulerian approach are performed. The developed models aim to mimic the subcooled nucleate boiling flows in concentric pipes, operating at high pressure. The predicted spatial fields of boiling quantities of interest are presented and commented. The numerical results are compared against the available experimental data, where it is shown that neglecting some interfacial forces like the lift or the wall lubrication forces will yield to good predictions for some quantities but will fail the prediction for others.  The models leading to the best predictions are highlighted and proposed as recommendations for future CFD simulations of subcooled nucleate boiling flows.