Although the renewable energy are attracting attention in the world, combustion still holds its special place in the energy supply of the power plant industries as well as in traortation, especially in aerospace equipments.Much research into combustion and enhancing the efficiency of this procees has done, but subjects related to envirionment are still interesting for many researchers. Therefore, in order to reduce combustion pollutants, it is necessary to predict the concentration of pollutants due to combustion. Since most of industrial burners operate in turbulent regimes and flow turbulence and chemical reactions of combustion has interaction on each other, hence, their interaction should be taken into account. In this research in order to simulate flow turbulence, chemical reactions in combustion and interaction between combustion and turbulence standard , laminar flamelet and probability density functions are used as the flamelet model, respectively. The flamelet model is a combustion model with a finite chemical reaction rate that has relatively high simulation precision in problems in which the chemical time scale is much smaller than the flow time scale and it needs less computional effort than models that use detailed chemical mechanism. In order to implement adiabatic and non-adiabatic flamelet models are formulated as UDF codes in C programming language. One of the necessities to utilize flamelet model is to generate a data base after solving a one or two dimensional laminar flame problem as pre-processing. In the present work, the flamelet used for adiabatic flamelet model survey, is generated using Cantera and Chemkin softwares in physical coordinates and Fluent software in in mixture fraction space. In adiabatic flamelet model, after solving traort equations for the laminar flame and flamelet generation, considering the effects of fluctuations is done with integrating using the probability density functions for the mixture fraction and the scalar dissipation rate and a library is generated in which thermochemical quantities are introduced as a function of mixture fraction mean and its variance and mean scalar dissipation rate. After generating the library, in order to simulate the turbulent flame, in each iteration, in addition to the traort equations commonly used in fluid dynamics problems, equations for the mean mixture fraction, mixture fraction variance, and mean scalar dissipation rate are also solved in each iteration and temperature and mass fraction of species in each computational cell in the domain with interpolation is calculated using these three quantities and this cycle will continue until convergence. In the simulation of some turbulent flames, due to the importance of heat transfer effects from flame or combustion chamber walls, the use of non-adiabatic flamelet model has been suggested by the researchers. In this study, the flamelet used in the non-adiabatic flamelet model is generated in the Cantera software and with considering the effect of thermal losses on the temperature distribution by help of a variable called enthalpy defect.This table has one more independent variable, enthalpy defect that is calculated by subtracting actual enthalpy obtained from energy equation from adiabatic enthalpy. Considering that non-adiabatic flamelet model can reduce the computational error just in and its needed time for simulation is nearly 20 times than adiabatic flamelet model; it is better to use adiabatic flamelet model in simulation of flames that do not have much heat loss. Keywords: Combustion, Steady laminar flamelet, Non-premixed, Bluff-body flame, Non-adiabatic