A presumed joint probability density function (pdf) model of turbulent combustion is proposed in this paper. The turbulent fluctuations of reactant concentrations and temperature are described using a presumed joint pdf of three-dimensional Gaussian distribution based on first and second-order moments of reactant concentration and temperature. Mean reaction rates in both premixed and diffusion combustion are obtained by mean of integration under the presumed joint pdf. This model is applied to predict turbulent premixed combustion of sudden-expansion flow and turbulent jet diffusion methane/air flame. For turbulent premixed combustion, the predicted results of temperature distribution and maximum temperature using the proposed model agree better with the experiment than that using the conventional eddy-breakup (EBU)-Arrhenius model. For the turbulent jet diffusion methane/air flame, the predicted results of velocity, temperature and species concentrations using the proposed model, the Arrhenius, EBU-Arrhenius, and laminar flamelet models are compared with experiment data. Results obtained with the presumed pdf model and that obtained by the laminar flamelet model both agree well with experiments, while results using the other models have a significant difference. The presumed joint pdf model is used to predict the NO formation process, which also agrees well with the experiment data. A unified turbulent combustion model, in which both effects of turbulent diffusion and chemical dynamics are considered, is established for both premixed and diffusion combustion, especially for the process of NO formation. (C) 2003 Elsevier Science Ltd. All rights reserved.