Reaksiyon içeren türbülanslı akışların modellenmesi için birleşik olasılık yoğunluk fonksiyonu (boyf) model denklemlerinin çözümünde kullanılan yeni bir gürbüz hibrit sonlu-hacimler/parçacık-tabanlı monte carlo metodu

A new robust hybrid finite-volume (FV)/particle method is developed for solving joint probability density function (JPDF) model equations of statistically stationary turbulent reacting flows. The method is designed to remedy the deficiencies of the hybrid algorithm developed by Muradoglu et al. [J. Comp. Phys.,154, 342(1999); J. Comp. Phys.,172, 842(2001)]. The density-based FV solver in the original hybrid algorithm has been found to be excessively dissipative and yet not very robust. To remedy these deficiencies, a pressure-based PISO algorithm in the open source FV package, OpenFOAM, is used to solve the Favre-averaged mean mass and momentum equations while a particle-based Monte Carlo algorithm is employed to solve the fluctuating velocity-turbulence frequency-composition JPDF transport equation. The mean density is computed as a particle field and passed to the FV method. Thus the redundancy of the density fields in the original hybrid method is removed making the new hybrid algorithm more consistent at the numerical solu- tion level. The new hybrid algorithm is first applied to simulate non-swirling cold and reacting bluff-body flows. The convergence of the method is demonstrated. In contrast with the original hybrid method, the new hybrid algorithm is very robust with respect to grid refinement and achieves grid convergence without any unphysical vortex shedding in the cold bluff-body flow case. In addition, the results are found to be in good agreement with the earlier PDF calculations and also with the available experimental data. Then the new hybrid algorithm is successfully applied to simulate the more complicated Sydney swirling bluff-body flame 'SM1'. The method is also very robust for this difficult test case and the results are in good agreement with the available experimental data. In all the cases, the PISO-FV solver is found to be highly resilient to the noise in the mean density field extracted from the particles. Fi- nally the method is coupled with the ISAT algorithm [53] for efficient treatment of detailed chemistry. Extensive calculations have been performed for the swirling and non-swirling bluff-body flames using ISAT.