The assumed β-pdf is widely used in turbulent flame simulation with the moment closure method. However, implementation of the β-pdf in turbulent flame simulation may sometimes encounter singularity difficulties in numerical calculation. The study proposes a robust β-pdf treatment to overcome these numerical difficulties. The proposed β-pdf treatment associated with the partial equilibrium chemistry model and the k-ϵ turbulence model is firstly applied to the case of H 2/N 2-air nonpremixed flame to demonstrate its robustness. Then, the present method is applied to the study of thermal NO formation in H 2-air nonpremixed flames. Comparison of results shows the present model to be capable of yielding fair agreement with the measurements. It is competitive with other more sophisticated methods, such as the pdf evolution method and the conditional moment closure method.