A new consistent hybrid algorithm for solution of the PDF equations of turbulent reactive flow

Abstract

This paper presents a newly developed consistent hybrid finite-volume (FV)/particle algorithm for solution of joint PDF (JPDF) model equation of turbulent reacting flows. In this approach, the open source FV package of OpenFOAM is employed to solve the Favre-averaged mean mass and momentum equations using pressure-based PISO algorithm while a particle-based Monte Carlo algorithm is used to solve the fluctuating velocity-turbulence frequency-compositions JPDF transport equation. In the earlier hybrid method [2, 3], a density-based FV algorithm was used to solve the mean flow equations but it has been found to be too dissipative and yet not very robust for incompressible or nearly incompressible flows mainly due to stiffness of the compressible flow equations in the low Mach number limit. In the this work, the density-based FV algorithm is first replaced with a pressure-based PISO algorithm to tackle this problem and then applied for simulation of the Sydney swirl stabilized bluff-body flame SM1. All the equations solved by the FV and particle algorithms are directly derived from the modeled JPDF transport equation so the present method is completely consistent at the level of governing equations. The position and velocity correction algorithms [3] are used to enforce full constancy at the numerical solution level. The results are found to be in a good agreement with the available experimental data and the recent computational results of De Meester et al. [1].