Approximation of joint PDFs by discrete distributions generated with Monte Carlo methods

Abstract

A new method for the approximation of multivariate scalar probability density functions (PDFs) in turbulent reacting flow by means of a joint presumed discrete distribution (jPDD) is presented. The jPDDs can be generated with specified mean values and variances as well as covariances. Correlations between variables -e.g. fluctuating mixture fractions and/or reaction progress -can thereby be taken into account. In this way the new approach overcomes an important limitation of ordinary presumed PDF methods, where statistical independence between the variables is often assumed. Different methods are presented to generate discrete distributions, based either on biased random number generators or on mixing models familiar from PDF transport models.The new approach is extensively validated on a turbulent flow configuration with simultaneous mixing and reaction. Large eddy simulation data as well as results from a transported PDF model are used for the validation of the jPDD approach. The comparison shows that in particular distributions generated with mixing models are able to predict mean reaction rates accurately. For the configuration considered, the neglect of correlations results in significant underestimation of reaction rates. Moreover it is found that higher statistical moments (e.g. the skewness) can influence reaction rates. The consequences for the generation of jPDDs are discussed.In summary, the new jPDD model has the potential to be significantly more accurate than established presumed PDF methods, because correlations between fluctuating variables can be taken into account. At the same time, the new approach is nearly as efficient as standard presumed PDF formulations, since mean rates are computed in a pre-processing step and stored in look-up tables as a function of the first and second moments of the relevant variables.

Trends in Multiphase Modeling and Simulation of Sprays

“…It is now widely acknowledged that this method is the most accurate EE representation of sprays [14]. However, these methods are computationally not tractable on real systems.…”

“…The QBMM approach can also be coupled with LES (instead of Lagrangian tracking of the dispersed phase elements) and this approach could transform spray computations in the near future [14].…”

Large Eddy Simulation of Autoignition in a Turbulent Hydrogen Jet Flame Using a Progress Variable Approach

“…The flamelet model with lookup table was generalized to more than one progress variable by De Goey and Ten Thije Boon Kkamp [19], who coined the popular term "flamelet generated manifolds." An adaptation to autoignition in a turbulent jet-in-cross-flow configuration, as it is found in reheat combustors, was proposed by Brandt et al [3,20,21], with a progress variable based on intermediate species formed during the induction period. The lookup table was built from 0D reactors with detailed chemistry, rather than reactivediffusive structures.…”

“…In the present paper, the model of Brandt et al [3,20,21] is developed further. Firstly, a progress variable based on a combination of intermediate and final product species is introduced, such that it describes both the buildup of radical species during the induction period and the subsequent heat release.…”

“…In the case of binary mixing, a beta PDF shape is most commonly used for the mixture fraction. For multivariate PDFs, the definition of a suitable functional form remains a challenge [21].…”