A mechanistic model for predicting critical heat flux (CHF) in rod bundles is developed based on the high-precision subchannel code ATHAS. To account for the non-uniform distribution of quality in subchannels caused by the presence of mixing vane grids (MVGs) and guide tubes (GTs), further subdivision of conventional subchannels and the incorporation of MVG crossflow models are implemented in ATHAS to provide a more detailed local flow field for the CHF model. The CHF mechanism model is developed based on the concepts of bubble crowding and assumption of local phenomenon. An equivalent-tube concept is proposed to convert the turbulent intensity distribution on the tube surface to that on the rod bundle surface. The present model is implemented in ATHAS and evaluated against CHF experimental data of rod bundles under the conditions of interest for pressurized water reactors (PWRs). The predicted results show good agreement with the experimental CHF data, with an average error and root mean square of 0.99% and 4.69%, respectively, for a total of 601 data points. Furthermore, owing to the use of the subchannel subdivision method, the present model does not require any cold wall effect correction and can accurately predict the radial position of CHF.