Computing reactive flows with a field Monte Carlo formulation and multi-scale methods Original Research Article

Two novel techniques are merged together to compute efficiently flows involving chemical reaction or combustion. In order to predict the evolution of the reacting flow, a probability density function approach is applied to the chemical species, which is described by a set of pure Eulerian transport equations with stochastic source terms. The proposed model, which has been extended in this case to a multi-scalar situation, has the advantage that the chemistry terms are closed, requiring no modeling. The formulation is solved by a stabilized finite element method based on the variational multi-scale theory, which enables an accurate and stable treatment of the difficult and stiff source terms employing iterative techniques. The combination of these ingredients results in a robust and accurate tool for the computation of reactive flows, as will be shown in the numerical examples.