Hydrogen combustion attracted much attention recently because of the need for a clean alternative energy. For the theoretical/numerical study of hydrogen combustion, there is a need for modeling capabilities for turbulent hydrogen flames. The present work examines the applicability of probability density function (pdf) turbulence models. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm is validated using experimental data for a heated turbulent plane jet. A study of H-2-F-2 diffusion flames has been carried out using this algorithm. Numerical results show that the pdf method is capable of correctly simulating turbulence effects on hydrogen combustion.