The present investigation considers the interactions among active sites or bubbles and the nonlinear dynamic characteristics of the boiling processes to analyze the stochastic and nonlinear features of boiling systems. A bubble size distribution function is derived by using bubble number balance. A new method is then proposed for predicting pool nucleate boiling heat transfer flux. The present model is in good agreement with previous experimental results

which indicates that the theory presented is more reasonable than traditional theories.