Soot particle size distribution reconstruction in a turbulent sooting flame with the split-based extended quadrature method of moments

Abstract

The Method of Moments (MOM) has largely been applied to investigate sooting laminar and turbulent flames. However, the classical MOM is not able to characterize a continuous particle size distribution (PSD). Without access to information on the PSD, it is difficult to accurately take into account particle oxidation, which is crucial for shrinking and eliminating soot particles. Recently, the Split-based Extended Quadrature Method of Moments (S-EQMOM) has been proposed as a numerically robust alternative to overcome this issue [Salenbauch et al., “A numerically robust method of moments with number density function reconstruction and its application to soot formation, growth, and oxidation,” J. Aerosol Sci. 128, 34–49 (2019)]. The main advantage is that a continuous particle number density function can be reconstructed by superimposing kernel density functions (KDFs). Moreover, the S-EQMOM primary nodes are determined individually for each KDF, improving the moment realizability. In this work, the S-EQMOM is combined with a large eddy simulation/presumed-probability density function flamelet/progress variable approach for predicting soot formation in the Delft Adelaide Flame III. The target flame features low/high sooting propensity/intermittency and comprehensive flow/scalar/soot data are available for model validation. Simulation results are compared with the experimental data for both the gas phase and the particulate phase. Good quantitative agreement has been obtained especially in terms of the soot volume fraction. The reconstructed PSD reveals predominantly unimodal/bimodal distributions in the first/downstream portion of this flame with particle diameters smaller than 100 nm. By investigating the instantaneous and statistical sooting behavior at the flame tip, it has been found that the experimentally observed soot intermittency is linked to mixture fraction fluctuations around its stoichiometric value that exhibits a bimodal probability density function.