A theoretical simulation was carried out for predicting the critical heat flux (CHF) of convective boiling for a round saturated water jet impingement on the stagnation zone of a hot surface. The study was focused on the effect of the solid–liquid contact angle on the CHF. A theoretical model based on the Long wave instability was applied to calculate the maximum liquid sub-film thickness under boiling bubbles and finally a semi-empirical and semi-theoretical correlation was proposed by combining the simulated calculation and the experimental data from the common metal heating surface. The correlation revealed the comprehensive effects of solid–liquid contact angles, jet velocity and jet diameter on the CHF and agreed well with the experimental data proposed by authors in the previous study.