LES of autoignition process in a turbulent droplet-laden mixing layer – impact of an evaporation model and discretisation method

Abstract

This paper focuses on the numerical modelling of the autoignition process in an ethanol-air mixture within a turbulent droplet-laden temporally evolving mixing layer, utilising the Large Eddy Simulation (LES) method. We compare results obtained by implementing two notably distinct discretisation methods (second-order Total Variation Diminishing – TVD, and fifth-order Weighted Essentially Non-Oscillatory – WENO) in the species and energy transport equations. Additionally, we examine three widely recognised evaporation models (the ‘$D^2$’ model, the Abramzon–Sirignano model, and a non-equilibrium model based on the Langmuir-Knudsen law). For modelling the combustion process, we employ the Eulerian Stochastic Field method and a laminar chemistry model. The simulations encompass two droplet diameter distributions, one with low and the other with significant diameter non-uniformity. Furthermore, nine distinct initial conditions are considered. Our findings indicate that adjusting the discretisation scheme or the evaporation model can yield simulation outcomes differing to a substantial degree, similar to the impact of modifying the actual properties of the spray droplets or their initial distribution. This is particularly pronounced in the case of sprays featuring a narrowed droplet diameter distribution. We found that the Abramzon–Sirignano model results in the most rapid autoignition. The TVD scheme contributes to a more even spatial distribution of the evaporated fuel, consequently facilitating quicker evaporation and faster autoignition. Using the WENO scheme generates higher fuel vapor gradients near the droplets, leading to less intense evaporation and slower fuel transport in space.

Keywords:

• LES
• spray combustion
• mixing layer
• numerical models

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was supported by the statutory funds of the Department of Thermal Machinery (Czestochowa University of Technology). Jakub Stempka received financial support for this research within the framework of the doctoral grant (no. 2019/32/T/ST8/00407) from the National Science Centre, Poland. The computations were carried out using PL-Grid (Cracow) and PCSS (Poznan) infrastructure.