Numerical investigation of spontaneous flame propagation under RCCI conditions

This paper presents results from one and two-dimensional direct numerical simulations under Reactivity Controlled Compression Ignition (RCCI) conditions of a primary reference fuel (PRF) mixture consisting of n-heptane and iso-octane. RCCI uses in-cylinder blending of two fuels with different autoignition characteristics to control combustion phasing and the rate of heat release. These simulations employ an improved model of compression heating through mass source/sink terms developed in a previous work by Bhagatwala et al. (2014), which incorporates feedback from the flow to follow a predetermined experimental pressure trace. Two-dimensional simulations explored parametric variations with respect to temperature stratification, pressure profiles and n-heptane concentration. Statistics derived from analysis of diffusion/reaction balances locally normal to the flame surface were used to elucidate combustion characteristics for the different cases. Both deflagration and spontaneous ignition fronts were observed to co-exist, however it was found that higher n-heptane concentration provided a greater degree of flame propagation, whereas lower n-heptane concentration (higher fraction of iso-octane) resulted in more spontaneous ignition fronts. A significant finding was that simulations initialized with a uniform initial temperature and a stratified n-heptane concentration field, resulted in a large fraction of combustion occurring through flame propagation. It was also found that the proportion of spontaneous ignition fronts increased at higher pressures due to shorter ignition delay when other factors were held constant. For the same pressure and fuel concentration, the contribution of flame propagation to the overall combustion was found to depend on the level of thermal stratification, with higher initial temperature gradients resulting in more deflagration and lower gradients generating more ignition fronts. Statistics of ignition delay are computed to assess the Zel'dovich (1980) theory for the mode of combustion propagation based on ignition delay gradients.