Experimental Flame Front Characterisation in a Lean Premix Burner Operating with Syngas Simpliﬁed Model Fuel

: The recent growing attention to energy saving and environmental protection issues has brought attention to the possibility of exploiting syngas from gasiﬁcation of biomass and coal for the ﬁring of industrial plants included in the, so called, Integrated Gasiﬁcation Combined Cycle power plants. In order to improve knowledge on the employ of syngas in lean premixed turbulent ﬂames, a large scale swirl stabilized gas-turbine burner has been operated with a simpliﬁed model of H 2 enriched syngas from coal gasiﬁcation. The experimental campaign has been performed at atmospheric pressure, with operating conditions derived from scaling the real gas turbines. The results are reported here and consist of OH-PLIF (OH Planar Laser Induced Fluorescence) measurements, carried out at decreasing equivalence of air / fuel ratio conditions and analysed together with the mean aerodynamic characterisation of the burner ﬂow ﬁeld in isothermal conditions obtained through LDV (Laser Doppler Velocimetry) and PIV (Particle Image Velocimetry) measurements. The OH concentration distributions have been analysed statistically in order to obtain information about the location of the most reactive zones, and an algorithm has been applied to the data in order to identify the ﬂame fronts. In addition, the ﬂame front locations have been successively interpreted statistically to obtain information about their main features and their dependence on the air to fuel ratio


Introduction
The attention to environmental protection and energy saving has increased in the last few decades, leading to more and more stringent regulations for power generation and transport emissions. Synthesis gas or "syngas" has shown to be a good opportunity to substitute, at least partially, fossil fuels in the near future. The syngas may be obtained through a large number of industrial activities, such as the gasification of coal, biomass and organic waste, or through steam reforming of natural gas [1]. Syngas is mainly composed of H 2 and CO, along with variable amounts of CH 4 and diluents, such as N 2 , CO 2 , and water vapour. The absence of particulate matter emissions and corrosive elements in the ashes constitutes a strong advantage if compared to conventional fuels, such as oils and coal [2]. The differences in the adopted gasification technology or in the coal quality may lead to large variations in the syngas composition [3]. Such a variation may affect several parameters of the syngas combustion, such as the flame speed [4,5], the ignition delay [6,7], and the extinction limit [8]. Consequently, the combustor design and operation may be largely complicated, especially in the case of gas turbine employ, such as in the Integrated Gasification Combined Cycle (IGCC) power plants. In fact, most of the stationary gas turbine combustors are employing lean premixed combustion technology (LPC) which has been shown to result in lower NO X emissions, if compared The test section was characterized by a variable geometry with the possibility of modifying its length by the insertion or removal of a 317 mm long cylindrical section, depending on the test operating conditions. Furthermore, the combustion chamber has, downstream of the test section strictly considered, a convergent segment driving the flow through a constricted cross-section, realized by means of the insertion of a steel septum, aimed at generating a pressure loss simulating the presence of the first stage of a gas turbine.
As previously mentioned, in order to allow in depth investigations of the combustive processes taking place within the combustion chamber, by means of laser-based non-intrusive measurement techniques, the test section was provided with properly designed optical accesses. Three 200 mm × 80 mm quartz glass windows were placed on three sides (right side, left side, and bottom side) of the chamber just downstream of the burner outlet.
In order to support the possibility of feeding multi-fuel burners equipped with multiple injection systems, the rig was provided with four independent fuel lines, fully instrumented with all the measurement and monitoring instrumentations and retroacted control devices. For the experimental activities presented here, the fuel lines were equipped with mass flow meters/controllers (IN-FLOW CTA model, maximum mass flow rate 140 kg/h) by Bronkhorst (High Tech B.V., Ruurlo-Netherlands), with an accuracy equal to ±0.5% of reading plus ±0.1% full scale.
The fuel lines can both be connected to the natural gas network and to alternative fuel sources (i.e., syngas or special gas mixtures tanks).

The burner prototype
Within the frame of the present study, a swirl stabilized burner prototype (a qualitative sketch is reported in Figure 2) was designed and manufactured by the Savona Combustion Laboratory, following the same philosophy underlying the 64.3 family of Ansaldo burners; however, it results As previously mentioned, in order to allow in depth investigations of the combustive processes taking place within the combustion chamber, by means of laser-based non-intrusive measurement techniques, the test section was provided with properly designed optical accesses. Three 200 mm × 80 mm quartz glass windows were placed on three sides (right side, left side, and bottom side) of the chamber just downstream of the burner outlet.
In order to support the possibility of feeding multi-fuel burners equipped with multiple injection systems, the rig was provided with four independent fuel lines, fully instrumented with all the measurement and monitoring instrumentations and retroacted control devices. For the experimental activities presented here, the fuel lines were equipped with mass flow meters/controllers (IN-FLOW CTA model, maximum mass flow rate 140 kg/h) by Bronkhorst (High Tech B.V., Ruurlo-Netherlands), with an accuracy equal to ±0.5% of reading plus ±0.1% full scale.
The fuel lines can both be connected to the natural gas network and to alternative fuel sources (i.e., syngas or special gas mixtures tanks).

The burner prototype
Within the frame of the present study, a swirl stabilized burner prototype (a qualitative sketch is reported in Figure 2) was designed and manufactured by the Savona Combustion Laboratory, following the same philosophy underlying the 64.3 family of Ansaldo burners; however, it results from a specifically-conceived aerodynamics and fuel injectors design. It is characterized by a nominal thermal power output (atmospheric pressure operations) of 500 kW and is equipped with two swirlers: An axial inner swirler, whose air mass flow rate is about 15% of the total swirler, and a centripetal outer swirler elaborating the remaining part (85%) of the air mass flow rate. from a specifically-conceived aerodynamics and fuel injectors design. It is characterized by a nominal thermal power output (atmospheric pressure operations) of 500 kW and is equipped with two swirlers: An axial inner swirler, whose air mass flow rate is about 15% of the total swirler, and a centripetal outer swirler elaborating the remaining part (85%) of the air mass flow rate. The burner can be fed with both natural gas ((NG) not exploited for the experimental activity presented here) and H2-enriched syngas (SG) by means of dedicated injection systems. The fuel nozzles are located slightly upstream of the inner axial swirler, which provides a diffusive flame aimed at the overall combustion stabilization. In contrast, in the centripetal swirler, the fuel nozzles are placed along the vanes (NG operations) or immediately downstream (SG operations), and such a flow generates the main premixed flame. In the present experimental campaign, only the SG injection system was used.
A cylindrical duct (whose radius is Rout = 59 mm) was placed downstream of the swirler assembly and collected the incoming flows before the combustion chamber inlet. A sketch of the burner is depicted in Figure 2, showing the airflow paths and the fuel nozzles locations.
The experimental investigations reported in the present work have been carried out in the combustor near field (close to the cylindrical duct outlet section) within the combustion chamber.

LDV, PIV, and LIF Diagnostic Systems
In order to characterize the flow field in the burner outlet region, the isothermal (350 °C) nonreactive velocity distribution was investigated by means of Laser Doppler Velocimetry and Particle Image Velocimetry measuring campaigns. The LDV instrumentation was composed of a 300 mW argon ion laser light source, operating at 514.5 nm (green) and 488 nm (blue). The transmitting probe had a front lens of 400 mm focal length and a beam separation of 38 mm, resulting in a geometrical control volume of 90 μm diameter and 1.9 mm length.
The probe was stiffly mounted on a three-axis computer-controlled traversing mechanism. The motion was transmitted to the carriages by stepping motors through a preloaded ball-screw assembly with a minimum linear translation step of 8 μm. A Bragg cell was used to apply a frequency shift (40 MHz) to one of each pair of beams and allows to solve directional ambiguity and to reduce angle bias. The flow was seeded with a 0.5-2.5 μm atomized spray of silicon oil injected in the settling chamber. The Doppler signals were transmitted by optical fibres to the receiving unit (Dantec The burner can be fed with both natural gas ((NG) not exploited for the experimental activity presented here) and H 2 -enriched syngas (SG) by means of dedicated injection systems. The fuel nozzles are located slightly upstream of the inner axial swirler, which provides a diffusive flame aimed at the overall combustion stabilization. In contrast, in the centripetal swirler, the fuel nozzles are placed along the vanes (NG operations) or immediately downstream (SG operations), and such a flow generates the main premixed flame. In the present experimental campaign, only the SG injection system was used.
A cylindrical duct (whose radius is R out = 59 mm) was placed downstream of the swirler assembly and collected the incoming flows before the combustion chamber inlet. A sketch of the burner is depicted in Figure 2, showing the airflow paths and the fuel nozzles locations.
The experimental investigations reported in the present work have been carried out in the combustor near field (close to the cylindrical duct outlet section) within the combustion chamber.

LDV, PIV, and LIF Diagnostic Systems
In order to characterize the flow field in the burner outlet region, the isothermal (350 • C) non-reactive velocity distribution was investigated by means of Laser Doppler Velocimetry and Particle Image Velocimetry measuring campaigns. The LDV instrumentation was composed of a 300 mW argon ion laser light source, operating at 514.5 nm (green) and 488 nm (blue). The transmitting probe had a front lens of 400 mm focal length and a beam separation of 38 mm, resulting in a geometrical control volume of 90 µm diameter and 1.9 mm length.
The probe was stiffly mounted on a three-axis computer-controlled traversing mechanism. The motion was transmitted to the carriages by stepping motors through a preloaded ball-screw assembly with a minimum linear translation step of 8 µm. A Bragg cell was used to apply a frequency shift (40 MHz) to one of each pair of beams and allows to solve directional ambiguity and to reduce angle bias. The flow was seeded with a 0.5-2.5 µm atomized spray of silicon oil injected in the settling chamber. The Doppler signals were transmitted by optical fibres to the receiving unit (Dantec Dynamics, Skovlunde-Denmark, Dual PDA 58N81 detector unit), where they were filtered to minimize the background noise. The signal transduced by the photomultipliers was collected by a Burst Spectrum Analyzer (Dantec Dynamics BSA P70 processor). The frequency range extended from 122 Hz up to 120 MHz, with accuracy better than 0.1% of the bandwidth. Dedicated software was used for the data post-processing. The measurements were performed on a radial traverse, located 40 mm downstream of the burner discharge section (Figure 2), with a radial extension equal to nearly twice the burner outlet radius R out . A 5 mm radial step was adopted. For each measuring point, 100,000 velocity samples were collected. Since a statistical bias could occur because the arrival times of the measurable particles were not statistically independent on the flow velocity, which brought them into the probe volume to correctly compute the average values, the data was weighted with the residence time of the particles in measuring volume. For the present experiment, the uncertainty of the instantaneous velocity was evaluated to be less than 1% [24,25].
The PIV instrumentation was composed by a double-cavity Nd:Yag pulsed laser (maximum energy 2 × 400 mJ per pulse at 532 nm, pulse duration 8 ns, repetition rate 10 Hz). The two beams, employed at nearly 220 mJ per pulse, were combined and formed a light sheet with a thickness of 2 mm, which entered the combustion chamber from the bottom window and was focused in a vertical meridional plane (Figure 1a). The light scattered by the seeding particles was recorded on a Dantec Dynamics HiSense MkII digital camera with a Charge-Coupled Device (CCD) matrix of 1344 × 1024 pixels. The camera's maximum frame rate in the double frame mode was 6 Hz, and the minimum interframe interval was 200 ns. The investigated area (not indicated in Figure 2 for sake of brevity) had an extension of 1.9 R out in the axial direction and of 2.6 R out in the radial direction, and was centred with respect of the burner axis. A total amount of 3000 couples of images were acquired. Once the PIV images were captured and digitized the velocity field was obtained, using a spatial cross-correlation function calculated over 64 × 64 pixels, 50% overlapped, interrogation areas. The experimental uncertainty for the instantaneous velocity was estimated to be 4%, according to the PIV error analyses of Grant [26] and Prasad et al. [27]. After a cross-correlation operation, data were post-processed to remove outliers, applying range scale validation and peak validation.
The laser-based techniques (LIF) measurement system exploited a laser source to excite the OH Q1(6) line in the A 2 Σ-X 2 Π (1,0) band, consisting of a diode pumped Nd:YAG laser (Quantel SA, Les Ulis-France, Brilliant B) capable of generating a pulsed (repetition rate up to 10 Hz) beam at 1064 nm wavelength with 850 mJ energy and equipped with a temperature-stabilized second-harmonic generator and a second beam at 532 nm wavelength with 400 mJ energy. The laser pump worked in combination with Tuneable Dye Laser utilizing, for the experimental activity presented here, a solution of Rhodamine 590 in ethanol. The output UV beam had energy equal to 15 mJ per pulse at 283 nm wavelength and was driven throughout an optical path to a standard UV-coated sheet forming optics, which generated a 50 mm long light sheet used to illuminate a vertical meridional plane in the combustor outlet region.
The light coming from the excited OH radicals was collected by an image intensified (Hamamatsu Photonics, Shizuoka-Japan, Image Intensifier Unit C-9546) camera (Dantec Dynamics HiSense Mk II), with a full resolution of 1344 × 1024 pixels and a quantum efficiency above 70% within the wavelength range 450-550 nm, adopting a narrowband UV-filter with a centre frequency of 309 nm. The camera allowed the pixel binning in 2 × 2 and 4 × 4 modality, in order to increase both sensitiveness and frame rate. The signal digitalisation was performed on a 12-bit basis in both single frame and double frame mode.
The synchronization of laser, camera, and image intensifier of the LIF system was controlled from the DynamicStudio (Dantec Dynamics) software platform via an 8-channel high-resolution (250 ps) synchronizer unit. The synchronizer unit was also equipped with two additional input channels, supporting the synchronization of measurements with external triggering signals.
The measurements were carried out in two axial locations (Figure 2), where the first was located as close as possible to the burner outlet section and the second was shifted axially downstream of 25 mm. The overall investigated area, after the sampled image post processing, had an extension of 1.09 R out in the axial direction and of 1.3 R out in the radial direction.

Air Flow Field Characterisation
The flow field generated by the burner was investigated in isothermal non-reactive conditions, representative of the nominal power output working point. Hence, the adopted mass flow rate was the same of such an operating condition (0.426 kg/s @ 500 kW). The air was preheated at 350 • C, which was the same preheating temperature used in the combustive tests. The flow field referred to unreactive conditions, even if characterised by the same values of air mass flow rate and similar preheating temperature, could not be uncritically assumed as stringently representative of the reactive aerodynamics in the combustion chamber. Nevertheless, the main flow features and the order of magnitude of the aerodynamic phenomena could be reasonably preserved when moving to the reactive conditions [28]. For such a reason, the unreactive flow characteristics were deemed useful to interpret the OH-PLIF results.
The radial profiles of the three mean velocity components u a , u r , and u θ (axial, radial, and tangential velocity components, respectively) made non-dimensional with the bulk velocity u bulk , measured at x/R out = 0.67, where x represents the axial direction, are reported in Figure 3. The u bulk velocity was obtained according to the following equation: where m represents the air mass flow rate.

Air Flow Field Characterisation
The flow field generated by the burner was investigated in isothermal non-reactive conditions, representative of the nominal power output working point. Hence, the adopted mass flow rate was the same of such an operating condition (0.426 kg/s @ 500 kW). The air was preheated at 350 °C, which was the same preheating temperature used in the combustive tests. The flow field referred to unreactive conditions, even if characterised by the same values of air mass flow rate and similar preheating temperature, could not be uncritically assumed as stringently representative of the reactive aerodynamics in the combustion chamber. Nevertheless, the main flow features and the order of magnitude of the aerodynamic phenomena could be reasonably preserved when moving to the reactive conditions [28]. For such a reason, the unreactive flow characteristics were deemed useful to interpret the OH-PLIF results.
The radial profiles of the three mean velocity components ua, ur, and uθ (axial, radial, and tangential velocity components, respectively) made non-dimensional with the bulk velocity ubulk, measured at x/Rout = 0.67, where x represents the axial direction, are reported in Figure 3. The ubulk velocity was obtained according to the following equation: where m represents the air mass flow rate. A zone characterized by negative axial velocity could be observed, ranging from the duct centreline to r/Rout = 0.29. For larger radii, such a component increases to a maximum equal to 1.1 ubulk located at r/Rout = 0.92, from which it decreases with a rapid decay to an outer zone, characterized by values close to 0.1 ubulk. Hence, the axial component profile indicated the presence of an inner recirculation zone (IRZ), typical in swirl-stabilized injectors [29][30][31][32][33][34], having, in the selected axial position, a radial extension corresponding to 29% of the burner outlet diameter. Similarly, the flow leaving the burner was concentrated in an annular jet, which showed a centrifugal behaviour, as suggested by the ur profile characterized by continuously positive values (the maximum value is equal to 0.36 ubulk at r/Rout = 0.92). The uθ profile indicated that the tangential velocity component behaved like a Rankine vortex with free vortex distribution in the external region and solid body rotation near the jet centreline. The flow turbulence was evaluated according to the following equation: A zone characterized by negative axial velocity could be observed, ranging from the duct centreline to r/R out = 0.29. For larger radii, such a component increases to a maximum equal to 1.1 u bulk located at r/R out = 0.92, from which it decreases with a rapid decay to an outer zone, characterized by values close to 0.1 u bulk . Hence, the axial component profile indicated the presence of an inner recirculation zone (IRZ), typical in swirl-stabilized injectors [29][30][31][32][33][34], having, in the selected axial position, a radial extension corresponding to 29% of the burner outlet diameter. Similarly, the flow leaving the burner was concentrated in an annular jet, which showed a centrifugal behaviour, as suggested by the u r profile characterized by continuously positive values (the maximum value is equal to 0.36 u bulk at r/R out = 0.92). The u θ profile indicated that the tangential velocity component behaved like a Rankine vortex with free vortex distribution in the external region and solid body rotation near the jet centreline. The flow turbulence was evaluated according to the following equation: where random fluctuation u' is obtained, considering the instantaneous velocity u in a generic position as the sum of the mean velocity plus its random fluctuation: u = u + u . For sake of simplicity, in the legend of Figure 3 and in the rest of the paper, the overbar indicating the mean value has been omitted. The Tu profile indicated that in the IRZ and the annular, the turbulence level was particularly high, with values in the interval 0.30 u bulk ÷ 0.4 u bulk and with a local maximum (0.32 u bulk ) placed at r/R out = 1.09, which was probably the shear layer of the outer part of the jet.
In order to provide a more immediate representation of the flow field in a meridional plane, a PIV measuring campaign was performed, and the results are reported in Figure 4 as contour plots of the axial velocity component and of Tu, with the vector field superposed. The turbulence intensity, which is represented, was obtained by neglecting the part related to the tangential component, since it cannot be measured by such a technique.
where random fluctuation u' is obtained, considering the instantaneous velocity u in a generic position as the sum of the mean velocity plus its random fluctuation: For sake of simplicity, in the legend of Figure 3 and in the rest of the paper, the overbar indicating the mean value has been omitted.
The Tu profile indicated that in the IRZ and the annular, the turbulence level was particularly high, with values in the interval 0.30 ubulk ÷ 0.4 ubulk and with a local maximum (0.32 ubulk) placed at r/Rout = 1.09, which was probably the shear layer of the outer part of the jet.
In order to provide a more immediate representation of the flow field in a meridional plane, a PIV measuring campaign was performed, and the results are reported in Figure 4 as contour plots of the axial velocity component and of Tu, with the vector field superposed. The turbulence intensity, which is represented, was obtained by neglecting the part related to the tangential component, since it cannot be measured by such a technique. The distribution of ua confirms that the flow was concentrated in an annular jet confined by two recirculation zones. Close to the burner outlet section, the jet showed large positive values which gradually decreased with axial coordinate increase. Moving downstream, the jet diverged and the cross section became larger. At the same time, due to turbulent mixing, the distribution tended to attenuate non-uniformities. Close to the burner outlet at large radii, a trace of the corner recirculation could have been seen, while across the duct axis, the IRZ recirculation zone, associated to the large tangential momentum and the sudden jet expansion in the combustion chamber, appeared evident. The maximum negative values of ua was around 0.4 ubulk located in the region 1 < r/Rout < 1.5. The superimposed vector plot provides an immediate representation of the meridional mean flow field, suggesting that the ring like jet may have been separated into two regions: An upstream region characterized by the reverse flow leaving the IRZ, and a downstream region where the flow entered the separation bubble (i.e., the IRZ) moving from the outer diameter towards the duct axis. The Tu distribution indicated that the largest random fluctuations of the axial and radial components were located at the interface between the annular jet and the internal separated flow.
Even if such a flow field analysis simply based on time-averaged values could not reveal the actual dynamic interaction between the IRZ and the surrounding high velocity region, it was able to provide the basic characteristics, which may be useful for the following analysis of the experimental characterisation of the combustion process, the object of the present work. The distribution of u a confirms that the flow was concentrated in an annular jet confined by two recirculation zones. Close to the burner outlet section, the jet showed large positive values which gradually decreased with axial coordinate increase. Moving downstream, the jet diverged and the cross section became larger. At the same time, due to turbulent mixing, the distribution tended to attenuate non-uniformities. Close to the burner outlet at large radii, a trace of the corner recirculation could have been seen, while across the duct axis, the IRZ recirculation zone, associated to the large tangential momentum and the sudden jet expansion in the combustion chamber, appeared evident. The maximum negative values of u a was around 0.4 u bulk located in the region 1 < r/R out < 1.5. The superimposed vector plot provides an immediate representation of the meridional mean flow field, suggesting that the ring like jet may have been separated into two regions: An upstream region characterized by the reverse flow leaving the IRZ, and a downstream region where the flow entered the separation bubble (i.e., the IRZ) moving from the outer diameter towards the duct axis. The Tu distribution indicated that the largest random fluctuations of the axial and radial components were located at the interface between the annular jet and the internal separated flow.
Even if such a flow field analysis simply based on time-averaged values could not reveal the actual dynamic interaction between the IRZ and the surrounding high velocity region, it was able to provide the basic characteristics, which may be useful for the following analysis of the experimental characterisation of the combustion process, the object of the present work.

OH-PLIF Measurement
Combustion tests were carried out in operating conditions corresponding to the working point chosen as the reference condition for the flow field characterisation. This is so that the same air mass flow rate and preheating temperature could be utilized (0.426 kg/s at 350 • C).
The tests were carried out in ultra-lean, partially premixed flame conditions, with the aim of exploring the effects on flame positioning, shape, and spatial stability due to air to fuel ratio modification in those that could reasonably be possible gas turbine conditions.
In particular, attention was dedicated to operations of the experimental burner with a hydrogen rich syngas model. The syngas composition was chosen to be 85% H 2 and 15% N 2 by volume (29% H 2 and 71% N 2 by mass, respectively). The chosen fuel model was provided already mixed and ready to be used in cylinders and its composition was certified by the supplier, having a preparation accuracy equal to ±0.10 ∆C/C (being C the nominal molar concentration of hydrogen in the mixture).
The results presented here refer to three different values of the equivalent air/fuel ratio (λ = α α st , ratio of actual air to fuel ratio to the stoichiometric one) at constant air mass flow rate (i.e., at a variable thermal power output).
Testing conditions are summarized in the following Table 1. The three operating conditions were obtained by scaling down the thermal power output of a real turbogas to atmospheric pressure. All tests were carried out addressing the 16% of total fuel mass flow rate to the nozzles feeding the central stabilizing, pilot, and diffusive flame.
In order to characterise the flame and its position, shape, and spatial stability, the OH-PLIF diagnostic technique was exploited.
According to many classical sources in scientific literature [35][36][37][38], theoretical calculations on laminar premixed flames provide information about the determination of molecules that can be possibly exploited as a tracer for the flame front position determination. Generally, three molecules, OH, CH, and CH 2 O, are the most frequently used in premixed turbulent flame experimental characterisations.
Independently of how good tracers of flame fronts could be recognised to be CH and formaldehyde (both showing limiting shortcomings, i.e., quite low concentrations of CH and possible significant formation of CH 2 O quite before the actual flame front), they are not exploitable in the present case due to that fact that no carbon kinetic was involved in the combustion of the utilized fuel.
On the other hand, the concentration of the OH radical is relatively high in premixed flames (both hydrocarbon/air flames and hydrogen/air flames) and, additionally, its concentration increases rapidly in the reaction zone. This rapid increase can be exploited in experimental approaches to mark the beginning of the flame front. The disadvantage of OH is that it is widely spread in the "post-flame" region, so that mere concentration maps of OH cannot be directly adopted for locating flame fronts. However, this limitation can be handled by adopting a proper flame front detection method, as is described in the following. Moreover, the utilization of laser-based techniques (LIF) for the OH in flame detection are well studied and accounted [16,[39][40][41][42][43][44][45][46][47][48]. Hence, all the results that are presented in the following will refer to OH-PLIF measurements.
It is well known that this technique provides data that is sensitive not only to the OH radical concentration, but also to the temperature field.
In light of this well-known equation: where I out is the fluorescence emission, I in the incident excitation, c dye the fluorescence dye concentration, ε(λ) the absorption coefficient, Φ(λ, T) the fluorescence quantum yield, and λ the radiation wavelength, OH-LIF measures either concentration distribution or temperature fields, assuming the other of the two parameters constant. Nevertheless, in the case of a flame, both parameters vary in the whole physical domain: The OH radicals (the dye, in this case) are produced in the flame front, heated up, then cooled down and finally destroyed. Despite no attempt to distinguish the effects due to the OH radical concentration variation from the ones due to temperature variation, this has not been considered a shortcoming in the flame front identification process. In fact, the two effects give concordant contributions to the detected fluorescence intensity. The raw image processing, which has been utilized in order to obtain information about shape, position, and spatial fluctuation, of the flame fronts is similar to the image processing used in [42] and [45].
A total of 600 OH-PLIF raw images were collected for each operating condition, in order to reduce the statistical uncertainty in the data. Such a number of samples was found to be a good trade-off between a sufficiently high number of samples and a reasonable data acquisition and post-processing time. The adopted procedure to post process the acquired images is indicated in Figure 5 and will be hereby explained.
Sets of 40 background images were collected without a flame, but with the laser running, in order to take into account both the radiation coming from the combustion chamber walls and the scattered laser light as background noise. One set of 20 images, corresponding to each of the operating conditions, was acquired before the measurement ran and 20 further images were acquired at the conclusion of the measurement run. Averaging the two sets gave the average background for each acquisition session. This average background image was used to correct the raw images.
In order to account for the axial unevenness of the laser sheet (due to possible misalignments or imperfect neatness of optics within the tuneable dye laser or along the optical path to the test section), a properly designed target was adopted to collect a set of images which allowed us to identify the mean laser sheet footprint, which resulted in the real axial distribution of light. Therefore, it was possible to determine a correction function, which was applied to each raw image to take into account such an unevenness. Since the axial extension of the collected images was set larger than the laser sheet, such a function was also used to black out the not illuminated (or not-enough illuminated) areas. The effects of this two-step post processing is shown in Figure 5a-c.
Subsequent filtering was required to reduce the noise content of the original OH-PLIF image, which in turn improved the accuracy of the subsequent flame-front detection. Smoothing of images by reducing the high-frequency pixilation noise was obtained by applying a 9 × 9 median filter (Figure 5d). Furthermore, to remove potentially spurious intensity fluctuations of the laser beam, the collected images have been filtered by analysing the intensity histogram and equalizing the max-min range to allow a correct and uniform use of the greyscale dynamic range (0÷255 basis for an 8 bit signal).
The effect of such a post processing is depicted in Figure 5e. The conditioned and balanced images were used to perform a statistical analysis of the data, obtaining the average distribution of OH and the RMS one (Figure 5j,k). Successively, the averaged and RMS distributions obtained from the two measurement locations were joined, as shown in Figure 5l,m. In order to account for the axial unevenness of the laser sheet (due to possible misalignments or imperfect neatness of optics within the tuneable dye laser or along the optical path to the test section), a properly designed target was adopted to collect a set of images which allowed us to identify the mean laser sheet footprint, which resulted in the real axial distribution of light. Therefore, it was possible to determine a correction function, which was applied to each raw image to take into account such an unevenness. Since the axial extension of the collected images was set larger than the laser sheet, such a function was also used to black out the not illuminated (or not-enough illuminated) areas. The effects of this two-step post processing is shown in Figure 5a-c.
Subsequent filtering was required to reduce the noise content of the original OH-PLIF image, which in turn improved the accuracy of the subsequent flame-front detection. Smoothing of images by reducing the high-frequency pixilation noise was obtained by applying a 9 × 9 median filter ( Figure 5d). Furthermore, to remove potentially spurious intensity fluctuations of the laser beam, the collected images have been filtered by analysing the intensity histogram and equalizing the max-min range to allow a correct and uniform use of the greyscale dynamic range (0÷255 basis for an 8 bit signal). The effect of such a post processing is depicted in Figure 5e. The conditioned and balanced images were used to perform a statistical analysis of the data, obtaining the average distribution of OH and the RMS one (Figure 5j,k). Successively, the averaged and RMS distributions obtained from the two measurement locations were joined, as shown in Figure 5l,m.
In order to extract the flame front from the conditioned and balanced images, a further postprocessing algorithm, based on the detection of the OH signal intensity gradient in the reaction zone, In order to extract the flame front from the conditioned and balanced images, a further post-processing algorithm, based on the detection of the OH signal intensity gradient in the reaction zone, was applied. In fact, a combination of median (already used in the step shown in Figure 5d) and non-linear diffusion (NLD) filters was used to condition the OH-PLIF images. The NLD filter required a threshold value, which was derived by the statistical treatment of the intensity histograms. According to [42] and [49], thresholds are set in a neighbourhood of the 15% of the maximum intensity of the mean distribution of each image set.
Hence, the application of the NLD filter [50,51] lead to a "remove outliers" effect (Figure 5f), where regions of the selected image, characterized by values larger than the threshold, were smoothed to a uniform value. In contrast, no smoothing occurred between or across the boundary of such regions. This kind of application of the NLD filter is well accounted as a powerful tool in the flame front detection process, as reported in [52] (concerning diffusive flames), [53] (in a two stroke internal combustion engine case), and [45] (for lean partially-premixed flames).
By using median and NLD filters, the flame-front edges were preserved, and their position was not artificially shifted. Furthermore, the pre-process, by means of a median filter, improved the quality of the NLD filter, reducing the noise content and lowering the number of erroneous flame fronts identified by the later processing.
Finally, the flame front detection was performed by applying a Sobel gradient detector (Figure 5g), allowing the separation of the image processing procedure from the edge detection operation. Among other gradient operators, the Sobel gradient detector represented a good choice as it had a simple implementation. Furthermore, by adopting a central differencing scheme on a 3 × 3 pixel matrix, it only considered pixel intensity in a close neighbourhood [54].
Successively, the single image flame fronts were processed with a simple statistical analysis obtaining the average distribution of the flame fronts (Figure 5h). Similar to the OH distributions, the average flame front distributions obtained from the two measurement locations were joined, as shown in Figure 5i.

Discussion
The contour plots of the mean distributions of the OH radical, as well as the RMS, are reported in Figures 6 and 7, where the three selected operating conditions are depicted from left to right for λ = 3.60, 3.25 and 2.90, respectively (thus, from left to right side, λ decreases). Moreover, to improve the clarity of the figures, a sketch of the burner was added (the burner cylindrical outlet was shortened in the picture for editing reasons). The spatial scale of the sketch and of the contour plots was kept equal, and the duct axis location was preserved. Consequently, this layout should be helpful in improving the correct location of the general features of the reacting areas (it is useful to remember that areas characterised by large values of OH concentration are not necessarily the same areas of the flame front location). Such large concentrations indicate that the OH radical kinetic was still far from being quenched. quality of the NLD filter, reducing the noise content and lowering the number of erroneous flame fronts identified by the later processing.
Finally, the flame front detection was performed by applying a Sobel gradient detector (Figure 5g), allowing the separation of the image processing procedure from the edge detection operation. Among other gradient operators, the Sobel gradient detector represented a good choice as it had a simple implementation. Furthermore, by adopting a central differencing scheme on a 3 × 3 pixel matrix, it only considered pixel intensity in a close neighbourhood [54].
Successively, the single image flame fronts were processed with a simple statistical analysis obtaining the average distribution of the flame fronts (Figure 5h). Similar to the OH distributions, the average flame front distributions obtained from the two measurement locations were joined, as shown in Figure 5i.

Discussion
The contour plots of the mean distributions of the OH radical, as well as the RMS, are reported in Figures 6,7, where the three selected operating conditions are depicted from left to right for λ = 3.60, 3.25 and 2.90, respectively (thus, from left to right side, λ decreases). Moreover, to improve the clarity of the figures, a sketch of the burner was added (the burner cylindrical outlet was shortened in the picture for editing reasons). The spatial scale of the sketch and of the contour plots was kept equal, and the duct axis location was preserved. Consequently, this layout should be helpful in improving the correct location of the general features of the reacting areas (it is useful to remember that areas characterised by large values of OH concentration are not necessarily the same areas of the flame front location). Such large concentrations indicate that the OH radical kinetic was still far from being quenched.  For all the working points, the reaction zones show some common features. An intense OH reacting zone was confined within the portion of the meridional plane close to the duct axis where the IRZ was present. Additionally, a second intense reaction region was located across the shear layer between the IRZ and the surrounding annular jet, although it is characterized by OH concentration values lower than the ones previously highlighted in the IRZ. The two main OH reaction zones showed opposite behaviours as λ decreased. In fact, the reaction region placed within the IRZ appeared to be progressively thinned in radial direction for enriched combustion conditions, whilst the reaction region located on the inner shear layer strengthened (i.e., the OH concentration increased) and widened for reduced equivalence air/fuel ratios. Such an effect was stronger at the large value of x/R out . Furthermore, at λ = 3.25 and 2.90, a third reaction area may have been detected and was located on the external shear layer (i.e., between the annular jet and the outer surrounding flow). The resulting scenario indicated that, for present experiments, the reduction of the equivalence ratio determined a broadening in the radial direction of the reaction area placed on the inner shear layer and the emerging of a new area located on the outer shear layer.
The analysis of the contour plots of the RMS of the OH concentrations, depicted in Figure 7, indicates that the largest fluctuations were located in the OH reaction zone close to the duct axis and particularly in its downstream part. In fact, for λ = 3.60, the largest fluctuations extended in an area whose axial range was 0.9 < x/R out < 1.1. For λ = 3.25. such an area results axially increased in the range of 0.75 < x/R out < 1.1. Finally, for λ = 2.90, large value fluctuations may have been detected already at x/R out = 0.4. A similar effect, even if less evident, was also present in the zone located across the inner shear layer. The general trend of the RMS contour plots suggested that the OH concentration fluctuations became higher as the equivalence air/fuel ratios were reduced. ratio determined a broadening in the radial direction of the reaction area placed on the inner shear layer and the emerging of a new area located on the outer shear layer.
The analysis of the contour plots of the RMS of the OH concentrations, depicted in Figure 7, indicates that the largest fluctuations were located in the OH reaction zone close to the duct axis and particularly in its downstream part. In fact, for λ = 3.60, the largest fluctuations extended in an area whose axial range was 0.9 < x/Rout < 1.1. For λ = 3.25. such an area results axially increased in the range of 0.75 < x/Rout < 1.1. Finally, for λ = 2.90, large value fluctuations may have been detected already at x/Rout = 0.4. A similar effect, even if less evident, was also present in the zone located across the inner shear layer. The general trend of the RMS contour plots suggested that the OH concentration fluctuations became higher as the equivalence air/fuel ratios were reduced. It is well known that the OH concentration contours can provide information on the location of the chemical activity as well as on the temperature field, but they cannot be exploited in their standard form for the flame front identification.
For this reason, the previously described algorithm for the flame front identification was applied to the OH-PLIF images. This resulted in a set of 600 instantaneous flame fronts for each of the investigated locations and operating conditions. Three instantaneous flame fronts referred to each of the operating conditions discussed here (from left to right, λ = 3.60, 3.25, and 2.90, respectively) are depicted as an example of flame front detection in Figure 8. It is well known that the OH concentration contours can provide information on the location of the chemical activity as well as on the temperature field, but they cannot be exploited in their standard form for the flame front identification.
For this reason, the previously described algorithm for the flame front identification was applied to the OH-PLIF images. This resulted in a set of 600 instantaneous flame fronts for each of the investigated locations and operating conditions. Three instantaneous flame fronts referred to each of the operating conditions discussed here (from left to right, λ = 3.60, 3.25, and 2.90, respectively) are depicted as an example of flame front detection in Figure 8. It is useful to remember that all the images were balanced and masked in order to make the investigated region representative of an even light distribution in the laser sheet. Some spurious effects may have been present close to the illuminated region limits.
In all the three operating conditions, the individuated flame fronts appeared uniformly corrugated along the entire frontiers, being assignable to the corrugated flamelets region in the wellknown Borghi/Peters diagram (not reported for sake of brevity). Moreover, the previously identified ( Figure 6) reaction areas located on the external shear layer for λ = 3.25 and 2.90 may have been ascribed to small sized flame fronts located in the same regions. The shape of the flame fronts (corrugated flamelets along all frontiers) was compatible with the high reactivity and diffusivity of hydrogen, even in significantly fuel-lean conditions.
Looking at the flame structures for the different burner loads, corresponding to the different equivalence air/fuel ratios, some observation could be done about the mutual interaction between It is useful to remember that all the images were balanced and masked in order to make the investigated region representative of an even light distribution in the laser sheet. Some spurious effects may have been present close to the illuminated region limits.
In all the three operating conditions, the individuated flame fronts appeared uniformly corrugated along the entire frontiers, being assignable to the corrugated flamelets region in the well-known Borghi/Peters diagram (not reported for sake of brevity). Moreover, the previously identified ( Figure 6) reaction areas located on the external shear layer for λ = 3.25 and 2.90 may have been ascribed to small sized flame fronts located in the same regions. The shape of the flame fronts (corrugated flamelets along all frontiers) was compatible with the high reactivity and diffusivity of hydrogen, even in significantly fuel-lean conditions.
Looking at the flame structures for the different burner loads, corresponding to the different equivalence air/fuel ratios, some observation could be done about the mutual interaction between turbulence and flame fronts. Considering the blow-up (Figure 9) of a portion, as an example, of the λ = 3.60 single shot flame front previously shown (consider that at a local and small scale level, the other flame fronts, even if at different equivalence air/fuel ratios, appear to behave similarly), it can be noticed that, at local and small-scale level, the flame front showed a typical feature of high velocity-lean hydrogen-enriched flame, whose convolutions broke up into islands of reactants and products [23].
The high local propagation velocity, as implied by hydrogen combustion, induced an increase of the turbulence characteristic scale, able to produce the flame front circumvolution. In accordance with the Gibson criterion, the turbulence length scales, effective in corrugating the flame front, had characteristic dimensional scales increased due to increased reactivity and propagation velocity [55]. Furthermore, a higher flame temperature produced a more intense action of flow "laminarization" [56], inhibiting the effect of the smallest turbulence scales in thickening the flame front. These two actions described above, lead to the not-so-corrugated or -thickened flame front observed, even in the case of a self-evident corrugated flamelet, such as the one considered here.
Moreover, the analysis of Figure 8 indicates the presence of small scale closed contour flame fronts surrounding the main flame front, which are also clearly visible in the enlarged area of Figure 9. Such a closed flame front became more evident for decreasing λ. In order to evaluate if the presence of the small scale closed fronts could be related to possible residual noise still affecting the processed images, the flame front detection procedure was repeated with different median filter sizes, but no significant changes in the flame fronts appearance was observed.
In order to provide a possible explanation for this characteristic, two consideration can be done. First of all, it has to be considered that, in the case of a swirl stabilized burner, as in the present case, both the motion flow field and the thermochemical phenomena taking place within it were characterised by significant tangential components, so that the motion of fresh gas towards the flame fronts must also be considered in the direction normal to the illuminated (measurement) plane. Furthermore, the flame fronts themselves were three-dimensional corrugated surfaces, presenting lobes that could be cut by the laser sheet and, as a consequence, projected on the measurement plane's closed contours.
In addition, the inception and progressive development of the so-called flamelets in the eddies phenomenon, as accounted by some authors [57], may be useful to provide an interpretation to the peculiar features of the instantaneous flame fronts shown here, making possible the interpretation of the small scale closed fronts, not only as the sectioning of a structure with a predominant development in the tangential direction, but also as actual parcelled out portion of flame front. Furthermore, a higher flame temperature produced a more intense action of flow "laminarization" [56], inhibiting the effect of the smallest turbulence scales in thickening the flame front. These two actions described above, lead to the not-so-corrugated or -thickened flame front observed, even in the case of a self-evident corrugated flamelet, such as the one considered here.
Moreover, the analysis of Figure 8 indicates the presence of small scale closed contour flame fronts surrounding the main flame front, which are also clearly visible in the enlarged area of Figure 9. Such a closed flame front became more evident for decreasing λ. In order to evaluate if the presence of the small scale closed fronts could be related to possible residual noise still affecting the processed images, the flame front detection procedure was repeated with different median filter sizes, but no significant changes in the flame fronts appearance was observed.
In order to provide a possible explanation for this characteristic, two consideration can be done. First of all, it has to be considered that, in the case of a swirl stabilized burner, as in the present case, both the motion flow field and the thermochemical phenomena taking place within it were characterised by significant tangential components, so that the motion of fresh gas towards the flame fronts must also be considered in the direction normal to the illuminated (measurement) plane. Furthermore, the flame fronts themselves were three-dimensional corrugated surfaces, presenting lobes that could be cut by the laser sheet and, as a consequence, projected on the measurement plane's closed contours.
In addition, the inception and progressive development of the so-called flamelets in the eddies phenomenon, as accounted by some authors [57], may be useful to provide an interpretation to the peculiar features of the instantaneous flame fronts shown here, making possible the interpretation of the small scale closed fronts, not only as the sectioning of a structure with a predominant development in the tangential direction, but also as actual parcelled out portion of flame front.
As reported in Figure 5j, for each operating condition, the flame front locations of each measuring area were averaged and joined. The results are depicted in Figure 10 and represented as contour plots, where the three selected operating conditions are depicted from left to right for λ = 3.60, 3.25, and 2.90, respectively (thus, from left to right side, λ decreases). It can be observed that the average siting of the flame front follows the averaged OH concentration maps. Thus, the flame fronts may be individuated mainly into the aforementioned two regions: One within the IRZ and the second one on the inner shear layer. Furthermore, at lower values of λ, traces of the presence of flame fronts could be individuated in correspondence with the external shear layer, at the interface between the swirling jet and the outer flow field. Nonetheless, such a presence of flame fronts appears weak and almost negligible, where the OH concentration in the same location appears more relevant. This allows the consideration that the outer shear layer (for lower global air to fuel ratios) was more interested in a protracting OH kinetic, elsewhere incepted, than in an independent combustion inception (though weakly witnessed).
The average flame front distributions were successively post-processed in order to individuate the local maxima as a function of the axial coordinate (independently within the two main flame front locations: Within the IRZ and on the inner shear layer). This allowed the finding of two lines (depicted in red in Figure 11), representing the most probable location (MPL) of the flame fronts (one for each of the two main reaction zones) for each operating condition. In fact, the flame fronts in a single snapshot were represented by a binarization of the greyscale image. As a consequence, the mean distribution of the flame fronts corresponds with a spatial probability density function, where higher pixel values corresponded to more flame fronts in the specific pixel [49,58]. allows the consideration that the outer shear layer (for lower global air to fuel ratios) was more interested in a protracting OH kinetic, elsewhere incepted, than in an independent combustion inception (though weakly witnessed). The average flame front distributions were successively post-processed in order to individuate the local maxima as a function of the axial coordinate (independently within the two main flame front locations: Within the IRZ and on the inner shear layer). This allowed the finding of two lines (depicted in red in Figure 11), representing the most probable location (MPL) of the flame fronts (one for each of the two main reaction zones) for each operating condition. In fact, the flame fronts in a single snapshot were represented by a binarization of the greyscale image. As a consequence, the mean distribution of the flame fronts corresponds with a spatial probability density function, where higher pixel values corresponded to more flame fronts in the specific pixel [49,58]. In Figure 10, the average distributions of flame fronts for the three operating conditions are depicted, positioned downstream of the burner (sketched) outlet. It can be noticed that the average siting of the flame front follows the averaged OH concentration maps. Thus, two main regions where fronts mainly can be found are individuated: One within the IRZ, the second on the inner shear layer. Furthermore, at lower values of λ, traces of the presence of flame fronts could be individuated in correspondence with the external shear layer, at the interface between the swirling jet and the outer motionless flow field. Nonetheless, this external presence of flame fronts appears weak and almost negligible, where the OH concentration in the same location appears more relevant. This allows considerating that the outer shear layer (for lower global air to fuel ratios) was more affected by the protracting of the OH kinetic, elsewhere incepted, than in an independent combustion inception (though weakly witnessed).
The average flame front distribution maps have been processed in order to individuate the signal local maxima (independently within the two main flame front locations: Within the IRZ and on the inner shear layer) as a function of the axial coordinate. This allowed the finding of two lines (depicted In Figure 10, the average distributions of flame fronts for the three operating conditions are depicted, positioned downstream of the burner (sketched) outlet. It can be noticed that the average siting of the flame front follows the averaged OH concentration maps. Thus, two main regions where fronts mainly can be found are individuated: One within the IRZ, the second on the inner shear layer. Furthermore, at lower values of λ, traces of the presence of flame fronts could be individuated in correspondence with the external shear layer, at the interface between the swirling jet and the outer motionless flow field. Nonetheless, this external presence of flame fronts appears weak and almost negligible, where the OH concentration in the same location appears more relevant. This allows considerating that the outer shear layer (for lower global air to fuel ratios) was more affected by the protracting of the OH kinetic, elsewhere incepted, than in an independent combustion inception (though weakly witnessed).
The average flame front distribution maps have been processed in order to individuate the signal local maxima (independently within the two main flame front locations: Within the IRZ and on the inner shear layer) as a function of the axial coordinate. This allowed the finding of two lines (depicted in red in Figure 11), representing the most probable location (MPL) of flame fronts for each operating condition.
Reporting the MPL lines on Cartesian diagrams ( Figure 11) allows the observation of different behaviours of the two reaction zones. Whilst the location of flame fronts on the shear layer appeared to be independent on the global air to fuel ratio, the MPL of flame fronts within the IRZ progressively moved towards the inner radii due to decreasing λ. In addition, for the operating conditions #1 and #2, it can be observed ( Figure 10) that, in the region within the IRZ and for x/Rout ranging from 0.1 to 0.5, an area which may have been identified as a dark nucleus, where the flame fronts were poorly probable, may be identified. As λ decreased, the probability of having flame fronts also within that nucleus grew to a uniform probability density of the presence of flame fronts in the volume enclosed within the MPL line, for λ = 2.90. This behaviour could suggest that premixed air and fuel, together with partially oxidised products brought back towards the burner outlet by the reversed flow within the IRZ more promptly, are ignited because of the more intense thermal release and higher temperatures in the presence of globally richer combustion.

Conclusions
The present investigation has been aimed at highlighting the main features and behaviours of the flame field obtained, employing a simplified syngas fuel model in a combustor, representative of turbogas real scale combustors, and operating in lean conditions at atmospheric pressure.
The flow field generated by the combustor was highly unsteady and characterized by a strong inner recirculation zone, which extended towards r/Rout = 0.52 at x/Rout = 0.95 and, hence. occupied a large portion of the downstream flow section. The turbulence level, based on the axial and radial velocity component fluctuations and on the bulk velocity, was characterized by large values across the shear layer, which developed between the ring-like jet and the IRZ, where Tu reached values close to 0.4. The OH-PLIF measurements were operated in three working points, characterized by equivalence air/fuel ratio values equal to λ = 3.60, 3.25. and 2.90, respectively. Overall, the investigated conditions of the reaction zones were characterized by some common features: A reacting zone confined within the inner recirculation zone and a second reaction zone located across the shear layer between the IRZ and the surrounding annular jet. Moreover, at λ = 3.25 and 2.90, a third reaction area was present and was located in the outer shear layer between the annular jet and the outer flow. The lowering of the equivalence air/fuel ratio leads to a reduction in the radial dimension of the reaction zone located in the IRZ, while it induces a broadening of the radial dimensions of the reaction zone placed across the inner shear layer and the emerging of new areas In addition, for the operating conditions #1 and #2, it can be observed ( Figure 10) that, in the region within the IRZ and for x/R out ranging from 0.1 to 0.5, an area which may have been identified as a dark nucleus, where the flame fronts were poorly probable, may be identified. As λ decreased, the probability of having flame fronts also within that nucleus grew to a uniform probability density of the presence of flame fronts in the volume enclosed within the MPL line, for λ = 2.90. This behaviour could suggest that premixed air and fuel, together with partially oxidised products brought back towards the burner outlet by the reversed flow within the IRZ more promptly, are ignited because of the more intense thermal release and higher temperatures in the presence of globally richer combustion.

Conclusions
The present investigation has been aimed at highlighting the main features and behaviours of the flame field obtained, employing a simplified syngas fuel model in a combustor, representative of turbogas real scale combustors, and operating in lean conditions at atmospheric pressure.
The flow field generated by the combustor was highly unsteady and characterized by a strong inner recirculation zone, which extended towards r/R out = 0.52 at x/R out = 0.95 and, hence. occupied a large portion of the downstream flow section. The turbulence level, based on the axial and radial velocity component fluctuations and on the bulk velocity, was characterized by large values across the shear layer, which developed between the ring-like jet and the IRZ, where Tu reached values close to 0.4. The OH-PLIF measurements were operated in three working points, characterized by equivalence air/fuel ratio values equal to λ = 3.60, 3.25. and 2.90, respectively. Overall, the investigated conditions of the reaction zones were characterized by some common features: A reacting zone confined within the inner recirculation zone and a second reaction zone located across the shear layer between the IRZ and the surrounding annular jet. Moreover, at λ = 3.25 and 2.90, a third reaction area was present and was located in the outer shear layer between the annular jet and the outer flow. The lowering of the equivalence air/fuel ratio leads to a reduction in the radial dimension of the reaction zone located in the IRZ, while it induces a broadening of the radial dimensions of the reaction zone placed across the inner shear layer and the emerging of new areas located on the outer regions.
The OH-PLIF distributions were post-processed with a suitable algorithm in order evaluate the instantaneous flame fronts, which appeared to be largely corrugated along their borders and, for decreasing λ, resulted in the presence of small-scale closed flame fronts.
A statistical analysis of the flame front data has allowed us to identify the most probable location of flame fronts in each operating condition and for each reaction area identified in the mean OH-PLIF distribution analysis. As λ was lowered, the MPL of flame fronts located in the IRZ progressively became thinner in the radial direction, while the flame front located along the shear layer between the IRZ end the ring-like jet remained unvaried.
As the data acquisitions were not triggered with the signals generated by the periodical phenomena potentially taking place within the test section, as typically is observed with a downstream a swirl stabilized burner, no period-dependant analysis of the combustive phenomena was carried out within the frame of the present work. This is the object of an ongoing understanding of the experimental activity presented here that will be the object of further works in the future.