A comparative study on partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI) in an optical engine

Partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI) are two new combustion modes in compression-ignition (CI) engines. However, the detailed in-cylinder ignition and flame development process in these two CI modes were not clearly understood. In the present study, firstly, the fuel stratification, ignition and flame development in PPC and RCCI were comparatively studied on a light-duty optical engine using multiple optical diagnostic techniques. The overall fuel reactivity (PRF number) and concentration (fuel-air equivalence ratio) were kept at 70 and 0.77 for both modes, respectively. Iso-octane and n-heptane were separately used in the port-injection (PI) and direct-injection (DI) for RCCI, while PRF70 fuel was introduced through direct-injection (DI) for PPC. The DI timing for both modes was fixed at –25°CA ATDC. Secondly, the combustion characteristics of PPC and RCCI with more premixed charge were explored by increasing the PI mass fraction for RCCI and using the split DI strategy for PPC. In the first part, results show that RCCI has shorter ignition delay than PPC due to the fuel reactivity stratification. The natural flame luminosity, formaldehyde and OH PLIF images prove that the flame front propagation in the early stage of PPC can be seen, while there is no distinct flame front propagation in RCCI. In the second part, the higher premixed ratio results in more auto-ignition sites and faster combustion rate for PPC. However, the higher premixed ratio reduces the combustion rate in RCCI mode and the flame front propagation can be clearly seen, the flame speed of which is similar to that in spark ignition engines but lower than that in PPC. It can be concluded that the ratio of flame front propagation and auto-ignition in RCCI and PPC can be modulated by the control over the fuel stratification degree through different fuel-injection strategies.     

Optical diagnostics on the reactivity controlled compression ignition (RCCI) with micro direct-injection strategy

Micro direct-injection (DI) strategy is often used to extend the operation range of the reactivity controlled compression ignition (RCCI) to high engine load, but its combustion process has not been well understood. In this study, the ignition and flame development of the micro-DI RCCI strategy were investigated on a light-duty optical engine using formaldehyde planar laser-induced fluorescence (PLIF) and high-speed natural flame luminosity imaging techniques. The premixed fuel was iso-octane and an oxygenated fuel of polyoxymethylene dimethyl ethers (PODE) was employed for DI. The fuel-air equivalence ratio of DI was kept at 0.09 and the premixed equivalence ratio was varied from 0 to 1. RCCI strategies with early and late DI timing at –25° and –5° crank angle after top dead center were studied, respectively. Results indicate that the early micro-DI RCCI features a single-stage high-temperature heat release (HTHR). The combustion in the low-reactivity region shows a combination of flame front propagation and auto-ignition. The late micro-DI RCCI presents a two-stage HTHR. The second-stage HTHR is owing to the combustion in the low-reactivity region that is dominated by flame front propagation when the premixed equivalence ratio approaches 1. For both early and late micro-DI RCCI, the intermediate-temperature heat release (ITHR) of iso-octane, indicated by formaldehyde, takes place in the low-reactivity region before the arrival of the flame front. This is quite different from the flame front propagation in spark-ignition (SI) engine that shows no ITHR in the unburned region. The DI fuel mass is a key factor that affects the combustion in the low-reactivity region. If the DI fuel mass is quite low, there is more possibility of flame front propagation; otherwise, sequential auto-ignition dominates. The emergence of the flame front propagation in micro-DI RCCI strategy reduces its combustion rate and peak pressure rise rate.     

Simultaneous 36 kHz PLIF/chemiluminescence imaging of fuel,CH2O and combustion in a PPC engine

The requirements on high efficiency and low emissions of internal combustion engines (ICEs) raise the research focus on advanced combustion concepts, e.g., premixed-charge compression ignition (PCCI), partially premixed compression ignition (PPCI), reactivity controlled compression ignition (RCCI), partially premixed combustion (PPC), gasoline compression ignition (GCI) etc. In the present study, an optically accessible engine is operated in PPC mode, featuring compression ignition of a diluted, stratified charge of gasoline-like fuel injected directly into the cylinder. A high-speed, high-power burst-mode laser system in combination with a high-speed CMOS camera is employed for diagnostics of the autoignition process which is critical for the combustion phasing and efficiency of the engine. To the authors’ best knowledge, this work demonstrates for the first time the application of the burst-system for simultaneous fuel tracer planar laser induced fluorescence (PLIF) and chemiluminescence imaging in an optical engine, at 36?kHz repetition rate. In addition, high-speed formaldehyde PLIF and chemiluminescence imaging are employed for investigation of autoignition events with a high temporal resolution (5 frames/CAD). The development of autoignition together with fuel or CH₂O distribution are simultaneously visualized using a large number of consecutive images. Prior to the onset of combustion the majority of both fuel and CH₂O are located in the recirculation zone, where the first autoignition also occurs. The ability to record, in excess of 100 PLIF images, in a single cycle brings unique possibilities to follow the in-cylinder processes without the averaging effects caused by cycle-to-cycle variations.     

Numerical simulation and validation of SI-CAI hybrid combustion in a CAI/HCCI gasoline engine

SI-CAI hybrid combustion, also known as spark-assisted compression ignition (SACI), is a promising concept to extend the operating range of CAI (Controlled Auto-Ignition) and achieve the smooth transition between spark ignition (SI) and CAI in the gasoline engine. In this study, a SI-CAI hybrid combustion model (HCM) has been constructed on the basis of the 3-Zones Extended Coherent Flame Model (ECFM3Z). An ignition model is included to initiate the ECFM3Z calculation and induce the flame propagation. In order to precisely depict the subsequent auto-ignition process of the unburned fuel and air mixture independently after the initiation of flame propagation, the tabulated chemistry concept is adopted to describe the auto-ignition chemistry. The methodology for extracting tabulated parameters from the chemical kinetics calculations is developed so that both cool flame reactions and main auto-ignition combustion can be well captured under a wider range of thermodynamic conditions. The SI-CAI hybrid combustion model (HCM) is then applied in the three-dimensional computational fluid dynamics (3-D CFD) engine simulation. The simulation results are compared with the experimental data obtained from a single cylinder VVA engine. The detailed analysis of the simulations demonstrates that the SI-CAI hybrid combustion process is characterised with the early flame propagation and subsequent multi-site auto-ignition around the main flame front, which is consistent with the optical results reported by other researchers. Besides, the systematic study of the in-cylinder condition reveals the influence mechanism of the early flame propagation on the subsequent auto-ignition.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

A generalized flame surface density modelling approach for the auto-ignition of a turbulent non-premixed system

Auto-ignition of turbulent non-premixed systems is encountered in practical devices such as diesel internal combustion engines. It remains a challenge for modellers, as it exhibits specific features such as unsteadiness, flame propagation and combustion far from stoichiometric conditions. In this paper, a two-dimensional DNS database of an igniting H₂/O₂/N₂ mixing layer, including detailed chemistry and transport, is extensively post-processed in order to gain physical insight into the flame structure and dynamics during auto-ignition. The results are used as a framework for the development of a generalized flame surface density modelling approach by integrating the equations over all possible mixture fraction values. The mean reaction rate is split into two contributions: a generalized flame surface density and a mean reaction rate per unit generalized flame surface density. The unsteadiness of the ignition phenomenon is accounted for via a generalized progress variable. Closures for the generalized surface average of the reaction rate and for the generalized progress variable are proposed, and the modelling approach is tested a priori versus the DNS data. The use of a laminar database for the chemistry coupled to the mean turbulent field via the generalized progress variable shows very promising results, capturing the correct ignition delay and the premixed peak in the turbulent mean heat release rate evolution. This allows confidence in future inclusion and validation of this approach in a RANS-CFD code.     

Numerical study of transient evolution of lifted jet flames: partially premixed flame propagation and influence of physical dimensions

Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.     

A numerical study of auto-ignition in turbulent lifted flames issuing into a vitiated co-flow

This paper presents a numerical study of auto-ignition in simple jets of a hydrogen–nitrogen mixture issuing into a vitiated co-flowing stream. The stabilization region of these flames is complex and, depending on the flow conditions, may undergo a transition from auto-ignition to premixed flame propagation. The objective of this paper is to develop numerical indicators for identifying such behavior, first in well-known simple test cases and then in the lifted turbulent flames. The calculations employ a composition probability density function (PDF) approach coupled to the commercial CFD code, FLUENT. The in-situ-adaptive tabulation (ISAT) method is used to implement detailed chemical kinetics. A simple k–ε turbulence model is used for turbulence along with a low Reynolds number model close to the solid walls of the fuel pipe.

The first indicator is based on an analysis of the species transport with respect to the budget of convection, diffusion and chemical reaction terms. This is a powerful tool for investigating aspects of turbulent combustion that would otherwise be prohibitive or impossible to examine experimentally. Reaction balanced by convection with minimal axial diffusion is taken as an indicator of auto-ignition while a diffusive–reactive balance, preceded by a convective–diffusive balanced pre-heat zone, is representative of a premixed flame. The second indicator is the relative location of the onset of creation of certain radical species such as HO₂ ahead of the flame zone. The buildup of HO₂ prior to the creation of H, O and OH is taken as another indicator of autoignition.

The paper first confirms the relevance of these indicators with respect to two simple test cases representing clear auto-ignition and premixed flame propagation. Three turbulent lifted flames are then investigated and the presence of auto-ignition is identified. These numerical tools are essential in providing valuable insights into the stabilization behaviour of these flames, and the demarcation between processes of auto-ignition and premixed flame propagation.     

Diesel flame lift-off stabilization in the presence of laser-ignition: a numerical study

Diesel flame lift-off and stabilization in the presence of laser-ignition were numerically investigated with the method of Eulerian stochastic fields. The aim was to scrutinise the interaction between the lifted diesel flame and an ignition kernel upstream of the lifted flame. The numerical simulation was carried out in a constant-volume combustion vessel with n-heptane as fuel. The process was studied previously in an experiment employing Diesel #2 as the fuel in the same combustion vessel. In the experiment a lifted flame was first established at a position downstream of the nozzle. An ignition kernel was then initiated using a high-energy pulse laser at a position upstream of the natural lift-off position of the diesel flame. The laser-ignition kernel was modelled using a high-temperature (～2000 K) hot spot. In both experiment and simulations the upstream front of the ignition kernel was shown to remain around the initial laser ignition site for a substantially long period of time, while the downstream front of the ignition kernel propagates rapidly towards the natural lift-off position downstream of the laser ignition site. The lift-off position oscillated before the final stabilization at the natural lift-off position. The structures and the propagation speed of the reaction fronts in the laser-ignition kernel and the main flame were analysed. Two different stabilization mechanisms, the auto-ignition mechanism and the flame propagation mechanism, were identified for the naturally lifted flame and the laser-induced reaction front, respectively. A mechanism was proposed to explain the oscillation of the lift-off position.     

Controlled detonation initiation in hypersonic flow

Experimental evidence of controlled detonation initiation and propagation in a hypersonic flow of premixed hydrogen-air is presented. This controlled detonation initiation is created in a hypersonic facility capable of producing a Mach 5 flow of hydrogen-air. Flow diagnostics such as high-speed schlieren and OH* chemiluminescence results show that a flame deflagration-to-detonation transition occurs as a combined result of turbulent flame acceleration and shock-focusing. The experimental results define three new distinct regimes in a Mach 5 premixed flow: deflagration-to-detonation transition (DDT), unsteady compressible turbulent flames, and shock-induced combustion. A two-dimensional implicit-LES (ILES) simulation, which solves the compressible, reactive Navier-Stokes equations on an adapting grid is conducted to provide additional insight into the local physical mechanism of detonation transition and propagation.     

Unsteady deflagration speed of an auto-ignitive dimethyl-ether (DME)/air mixture at stratified conditions

The propagation speed of an auto-ignitive dimethyl-ether (DME)/air mixture at elevated pressures and subjected to monochromatic temperature oscillations is numerically evaluated in a one-dimensional statistically stationary configuration using fully resolved numerical simulations with reduced kinetics and transport. Two sets of conditions with temperatures within and slightly above the negative temperature coefficient (NTC) regime are simulated to investigate the fundamental aspects of auto-ignition and flame propagation along with the transition from auto-ignitive deflagration to spontaneous propagation regimes under thermal stratification. Contrary to the standard laminar flame speed, the steady propagation speed of an auto-ignitive front is observed to scale proportionally to its level of upstream reactivity. It is shown that this interdependence is primarily influenced by the characteristic residence time and the homogeneous auto-ignition delay. Furthermore, the unsteady reaction front in either of the two cases responds distinctly to the imposed stratification. Specifically, the results in both cases show that the dynamic flame response depends on the mean temperature at the flame base T_b and the time-scale of thermal stratification. It is also found that, based on T_b and the propensity of the mixture to two-stage chemistry, the instantaneous peak propagation speed and the overall time taken to achieve that speed differs considerably. A displacement speed analysis is carried out to elucidate the underlying combustion modes that are responsible for such a variation in flame response.     

Dynamics of triple-flames in ignition of turbulent dual fuel mixture: A direct numerical simulation study

Pilot-ignited dual fuel combustion involves a complex transition between the pilot fuel autoignition and the premixed-like phase of combustion, which is challenging for experimental measurement and numerical modelling, and not sufficiently explored. To further understand the fundamentals of the dual fuel ignition processes, the transient ignition and subsequent flame development in a turbulent dimethyl ether (DME)/methane-air mixing layer under diesel engine-relevant conditions are studied by direct numerical simulations (DNS). Results indicate that combustion is initiated by a two-stage autoignition that involves both low-temperature and high-temperature chemistry. The first stage autoignition is initiated at the stoichiometric mixture, and then the ignition front propagates against the mixture fraction gradient into rich mixtures and eventually forms a diffusively-supported cool flame. The second stage ignition kernels are spatially distributed around the most reactive mixture fraction with a low scalar dissipation rate. Multiple triple flames are established and propagate along the stoichiometric mixture, which is proven to play an essential role in the flame developing process. The edge flames gradually get close to each other with their branches eventually connected. It is the leading lean premixed branch that initiates the steady propagating methane-air flame. The time required for the initiation of steady flame is substantially shorter than the autoignition delay time of the methane-air mixture under the same thermochemical condition. Temporal evolution of the displacement speed at the flame front is also investigated to clarify the propagation characteristics of the combustion waves. Cool flame and propagation of triple flames are also identified in this study, which are novel features of the pilot-ignited dual fuel combustion.     

Simulations and experiments on the ignition probability in turbulent premixed bluff-body flames

The ignition characteristics of a premixed bluff-body burner under lean conditions were investigated experimentally and numerically with a physical model focusing on ignition probability. Visualisation of the flame with a 5 kHz OH* chemiluminescence camera confirmed that successful ignitions were those associated with the movement of the kernel upstream, consistent with previous work on non-premixed systems. Performing many separate ignition trials at the same spark position and flow conditions resulted in a quantification of the ignition probability P_ign, which was found to decrease with increasing distance downstream of the bluff body and a decrease in equivalence ratio. Flows corresponding to flames close to the blow-off limit could not be ignited, although such flames were stable if reached from a richer already ignited condition. A detailed comparison with the local Karlovitz number and the mean velocity showed that regions of high P_ign are associated with low Ka and negative bulk velocity (i.e. towards the bluff body), although a direct correlation was not possible. A modelling effort that takes convection and localised flame quenching into account by tracking stochastic virtual flame particles, previously validated for non-premixed and spray ignition, was used to estimate the ignition probability. The applicability of this approach to premixed flows was first evaluated by investigating the model's flame propagation mechanism in a uniform turbulence field, which showed that the model reproduces the bending behaviour of the S_T-versus-u′ curve. Then ignition simulations of the bluff-body burner were carried out. The ignition probability map was computed and it was found that the model reproduces all main trends found in the experimental study.     

Auto-ignition and flame stabilization of pulsed methane jets in a hot vitiated coflow studied with high-speed laser and imaging techniques

The auto-ignition of a pulsed methane jet issuing into a laminar coflow of hot exhaust products of a lean premixed hydrogen/air flat flame was examined using high-speed laser and optical measurement techniques with frame rates of 5 kHz or more. OH* chemiluminescence was used to determine the downstream location of the first auto-ignition kernel as well as the stabilization height of the steady-state lifted jet flame. OH planar laser-induced fluorescence (PLIF) was used to determine further details of the auto-ignition with a higher spatial resolution. Simultaneous imaging of broadband luminosity from a viewing angle perpendicular to the OH* chemiluminescence was applied, to three-dimensionally reconstruct the ignition kernel location in space and to determine whether the first occurrence of the kernel was within or beyond the PLIF laser sheet. The development and expansion of the jet was characterized by high-speed Schlieren imaging. Statistics have been compiled for both the ignition time as well as the downstream location of the first auto-ignition kernel and the stabilization height of the steady-state lifted jet flame. From the PLIF images it was found that auto-ignition tended to occur at the interface between bulges of the inflowing jet and the coflow. For steady-state conditions, auto-ignition kernels were observed frequently below the flame base, emphasizing that the lifted jet flame is stabilized by auto-ignition.     

Characteristics of auto-ignition in a stratified iso-octane mixture with exhaust gases under homogeneous charge compression ignition conditions

Ignition and propagation of a reaction front in a counterflow system of an iso-octane/air stream mixing with an exhaust gas stream is computationally investigated to understand the fundamental characteristics of homogeneous charge compression ignition (HCCI) auto-ignition. Various mixing rates are imposed on the system and the effects of dissipation rates on auto-ignition are studied. Ignition delay and front propagation speed across the mixing layer are determined as a function of a local mixture fraction variable. The results show that mixture inhomogeneity and dissipation rate have a significant influence on ignition. Diffusive transport is found to either hamper or advance ignition depending on the initial reactivity of the mixture. Based on the relative importance of diffusion on ignition front propagation, two distinct ignition regimes are identified: the spontaneous ignition regime and the diffusion-controlled regime. The transition between these two regimes is identified using a criterion based on the ratio of the timescales of auto-ignition and diffusion. The results show that ignition in the spontaneous regime is more likely under typical HCCI operating conditions with iso-octane due to its high reactivity. The present analysis provides a means to develop an improved modelling strategy for large-scale engine simulations.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

The impact of N2 micro-jets on the V-to-M flame shape transition in a premixed swirl burner

Injection of N₂ through micro-jets located on the dump plane of a lean premixed swirl stabilized combustor is investigated as a new method for mitigating combustion instabilities. This study focuses on the chemical and fluid dynamic processes by which the N₂ micro-jets impact the flame dynamics. An experimental and numerical investigation is performed to characterize the combustion instability during the V-to-M flame shape transition in a swirl burner fueled with premixed CH₄/air, at an equivalence ratio of 0.62. Reasonable agreements have been found between the experimental measurements and simulation results. Both of them present that the flame changes from V-shape to M-shape periodically, and a low-frequency instability around 10 Hz is observed accordingly. It is confirmed that intermittent flame extinction in the outer recirculation zone (ORZ) is the source of the combustion instability. Furthermore, injection of N₂ through micro-jets located on the combustor dump plane, into the outer recirculation zone, results in a stable V shape flame. It is clearly seen that the ORZ dilution can eliminate the combustion instability without inhibiting the combustion efficiency. A special focus is placed on the impact of the diluent injection on the local flame-flow interaction. The nitrogen micro-jets increase the local nitrogen concentration by 7% on average, lowering the flame speed and extinction strain rates by 27% and 17% respectively. Moreover, the micro-jets increase the turbulence intensity in the ORZ, leading to a significant increase in the Karlovitz number and transferring the local combustion regime from the thin reaction zone regime to the broken reaction zone regime. Hence, the nitrogen micro-jets impact on both the turbulence and the chemical reaction rates prevents flame propagation into the ORZ and results in a stable flame.     

LES/FGM investigation of ignition and flame structure in a gasoline partially premixed combustion engine

This paper presents a joint numerical and experimental study of the ignition process and flame structures in a gasoline partially premixed combustion (PPC) engine. The numerical simulation is based on a five-dimension Flamelet-Generated Manifold (5D-FGM) tabulation approach and large eddy simulation (LES). The spray and combustion process in an optical PPC engine fueled with a primary reference fuel (70% iso-octane, 30% n-heptane by volume) are investigated using the combustion model along with laser diagnostic experiments. Different combustion modes, as well as the dominant chemical species and elementary reactions involved in the PPC engines, are identified and visualized using Chemical Explosive Mode Analysis (CEMA). The results from the LES-FGM model agree well with the experiments regarding the onset of ignition, peak heat release rate and in-cylinder pressure. The LES-FGM model performs even better than a finite-rate chemistry model that integrates the full-set of chemical kinetic mechanism in the simulation, given that the FGM model is computationally more efficient. The results show that the ignition mode plays a dominant role in the entire combustion process. The diffusion flame mode is identified in a thin layer between the ultra fuel-lean unburned mixture and the hot burned gas region that contains combustion intermediates such as CO. The diffusion flame mode contributes to a maximum of 27% of the total heat release in the later stage of combustion, and it becomes vital for the oxidation of relatively fuel-lean mixtures.     

Large eddy simulation of multi-regime combustion with a two-progress variable approach for carbon monoxide

Simulations of two cases in a novel multi-regime burner configuration are undertaken using a presumed joint probability density function (PDF) approach with tabulated chemistry. The flame conditions are varied by changing the central jet equivalence ratio, which produces different multi-regime combustion modes in the non-premixed inner flame. An outer premixed flame and recirculation zone behind a bluff body are present to supply heat and combustion products to stabilise the inner flame. A two-progress variable approach is tested to improve predictions of carbon monoxide (CO) in the post-flame regions, where CO oxidation occurs. The large eddy simulation set-up and sub-grid combustion model are assessed through comparisons with time-averaged measurements for radial profiles at different streamwise locations. The jet break-up length, the shear layers and the mixture fraction distribution are well captured in both cases. The temperature distribution is well captured for the inner flame in each case but the temperature and mixture fraction are over predicted in the downstream regions of the outer premixed flame, which is due to increased dilatation that suppresses air entrainment. Improved predictions of the CO mass fraction are obtained for the outer premixed flames with the two-progress variable approach. Over predictions are seen in the upstream regions of the inner flame when the CO mass fraction is obtained from a look-up table, suggesting that the CO mass fraction should be transported to include the convection/diffusion balance in regions where there is no flame. Furthermore, transporting the CO mass fraction with a one-progress variable approach produces over predictions in the burnt regions, suggesting a two-progress variable model is needed to capture the consumption region of CO. The multi-regime combustion characteristics are observed to be stronger in flame MRB26b, where non-premixed and rich premixed combustion is present. For flame MRB18b, the non-premixed contribution is smaller and weak stratified combustion is observed.