Turbulence–chemistry interaction models such as the conditional moment closure and various flamelet models require a presumed Probability Density Function (PDF) model of the conditioning variables. In turbulent stratified flames, a joint-PDF model for a reaction progress variable and the mixture fraction is required. The joint-PDF is often modelled by two beta functions using the first and second moments of the two conditioning variables. In this work, the performance of a joint-PDF model based on the flamelet PDF approach is compared to the double beta-PDF model. The conditional PDF of the reaction progress variable conditioned on mixture fraction is modelled with the flamelet PDF. Unlike the beta PDF, the flamelet functional form for the conditional PDF varies with the local mixture fraction. The two PDF models are coupled with a two-dimensional premixed flamelet-generated manifold chemistry model. Two stratified flames of the Cambridge–Sandia bluff-body burner are simulated using a Reynolds-Averaged Navier–Stokes (RANS) approach. The results indicate that the flamelet PDF model has a superior, albeit marginal, predictions of the temperature, major and minor species mass fractions.