“…Whatever the scalar is chemically-reactive (e.g., progress variable) or not (e.g., mixture fraction), its variance as well as the shape of its probability density function (PDF) are indeed settled, to a large extent, by the scalar dissipation rate (SDR), see for instance [15,16,17,18]. The SDR is the product of the scalar diffusivity with the square of the scalar gradient norm [19,20,21] and, the studies reported in references [21,22,23,24,25,26] excepted, it is noteworthy that there were quite few efforts directed towards its detailed analysis in turbulent reactive two-phase flows. Moreover, it should be emphasized that most of the seldom studies that have been conducted in this direction were focused on the direct effect of the additional source terms associated with the evaporation process, see for instance [23,24,25].…”