A fractal model for nucleate pool boiling of nanofluids at high heat flux and critical heat flux (CHF) is developed based on the fractal distribution of nanoparticles and nucleation sites on boiling surfaces in this paper. The formula of calculating high heat flux and CHF for nanofluids in nucleate pool boiling is given by taking into account heat convection between nanoparticles and liquids due to the Brownian motion of nanoparticles in fluids. The proposed model is expressed as a function of temperature of nanofluids, the effective thermal conductivity of nanofluids, the average size of nanoparticles, the fractal dimension of nanoparticles and nucleation sites, the nanoparticles volume fraction of suspension, and physical properties of fluids. No additional/new empirical constant is introduced in this fractal model. An agreement between the proposed model predictions and experimental data is found. The validity of the fractal model for nucleate pool boiling of nanofluids at high heat flux and CHF is thus verified.