Exploring Vortex–Flame Interactions and Combustion Dynamics in Bluff Body-Stabilized Diffusion Flames: Effects of Incoming Flow Velocity and Oxygen Content

: An afterburner encounters two primary features: high incoming flow velocity and low oxygen concentration in the incoming airflow, which pose substantial challenges and contribute significantly to the deterioration of combustion performance. In order to research the influence of oxygen content on the dynamic combustion characteristics of the afterburner under various inlet velocities, the effect of oxygen content (14–23%) on the field structure of reacting bluff body flow, flame morphology, temperature pulsation, and pressure pulsation of the afterburner at different incoming flow velocities (0.1–0.2 Ma) was investigated in this study by using a large eddy simulation method. The results show that two different instability features, BVK instability and KH instability, are observed in the separated shear layer and wake, and are influenced by changes in the O 2 mass fraction and Mach number. The oxygen content and velocity affected the oscillation amplitude of the downstream flow. As the O 2 mass fraction decreases, the flame oscillation amplitude increases, the OH concentration in the combustion chamber decreases, and the flame temperature decreases. Additionally, the amplitude of the temperature pulsation in the bluff body flame was primarily influenced by the temperature intensity of the flame and BVK instability. Moreover, the pressure pulsation is predominantly affected by the dynamic characteristics of the flow field behind the bluff body. When the BVK instability dominated, the primary frequency of the pressure pulsation aligned with that of the temperature pulsation. Conversely, under the dominance of the KH instability, the temperature pulsation did not exhibit a distinct main frequency. At present, the influence of oxygen content and incoming flow rate on the combustion performance of the combustion chamber is not clear. The study of the effect of oxygen content on the combustion characteristics of the combustion chamber at different incoming flow rates provides a reference for improving the performance of the combustion chamber and enhancing the combustion stability.


Introduction
The afterburner is a critical component for thrust augmentation in military aircraft engines [1], enabling a substantial increase in the engine thrust from 60% to 70%.The airflow entering the afterburner is the exhaust gas from the main combustor in the aero-engine, which contains one-quarter less oxygen than pure air.Low oxygen content is detrimental to stable combustion [2][3][4][5], significantly reducing the combustion reaction rate and flame propagation speed in the afterburner [6][7][8].Moreover, the incoming flow velocity in the afterburner is exceptionally high, and minor perturbations can potentially lead to combustion instability [9][10][11][12][13].Combustion instability manifests as periodic pressure and temperature pulsations [14][15][16], posing a risk of severe mechanical vibration and localized overheating of the aero-engine.This can damage key components of the afterburner, thereby affecting the overall engine combustion performance [17][18][19][20].Therefore, it is of great significance to explore the influence of the O 2 mass fraction on the dynamic characteristics of afterburner flames under various incoming flow velocities.
Given the substantial impact of the incoming flow conditions on bluff body flames, numerous studies have investigated the phenomenon of combustion instability in bluff body flames.Hosseini et al. [21] analysed the influence of the O 2 volume fraction on the combustion performances of micro-flameless combustion through numerical simulations.Their findings revealed a decrease in both the maximum and exit temperatures of the bluff body micro-flameless combustion as the O 2 volume fraction decreased from 7% to 5%.Noor et al. [22] examined the influence of the O 2 mole fraction on the flameless combustion performances of a bluff body.They found that decreased oxygen content resulted in a more uniform temperature distribution within the combustion chamber.Liu et al. [23] studied bluff body flame characteristics under low-oxygen-content conditions.The results indicated that the temperature and OH concentration in the reaction zone decreased at lower oxygen contents, suggesting weaker reactivity in the reaction zone.Roy et al. [24] analysed the influence of diluent addition on bluff flames using numerical simulation methods.The temperature of the flame reduced as the oxygen content decreased, and an increase in the diluent concentration and incoming flow velocity led to the flame approaching extinction.Mishra et al. [25] introduced N 2 into a bluff flame and observed that a decrease in the oxygen content resulted in a decrease in the flame temperature and an increase in the length of the bluff flame.Shanbhogue et al. [26] conducted experimental research on the impact of the incoming flow velocity on premixed bluff flames.Their findings revealed that the preferential diffusion effect of the premixed bluff flames intensified with increasing velocity.Wan et al. [27] studied the influence of the incoming flow velocity on CH 4 /air bluff flames using numerical simulations and found that the area of the recirculation zone behind the bluff body expands with increasing incoming velocity.A high incoming flow velocity induced a high strain rate, ultimately resulting in flame extinction.Kalathoor et al. [28] investigated the impact of the incoming flow Mach number on combustion instability in an afterburner using a numerical simulation method.Their findings revealed that with an increase in the Reynolds number, the dominant frequency shifted from the acoustic mode to the hydrodynamic mode, and the length of the heat release rate zone became shorter.
Lapenna et al. [29] investigated the combustion instability of propane-air flames by experimental and numerical simulation methods, and showed that the flame is mainly affected by Darrieus-Landau (DL) instability at low Re numbers, and is gradually dominated by turbulence as the Re number increases.Emerson et al. [30] investigated the influence of the density ratio on reacting bluff body flow field characteristics.Their results indicated that as the density ratio changes, the flame morphology behind the bluff body transitions from Bénard-von Kármán (BVK) to Kelvin-Helmholtz (KH) instability.Additionally, Erickson et al. [31] studied the influence of the temperature ratio on the dynamics of bluff body-stabilized flames.Their research also demonstrated that as the temperature ratio changes, the flame transitions between BVK and KH instabilities.Deng et al. [32] used numerical simulation to study the effects of incoming temperature and oxygen content on the flame instability of a bluff, and found that the incoming temperature and oxygen content significantly affect the position of the interconversion between the KH instability and the BVK instability.
The influence of the oxygen content on the dynamic combustion characteristics of bluff body flames at different inlet velocities is not yet fully understood.To investigate how the oxygen content affects the dynamic combustion characteristics of the bluff body flame under various inlet velocities, this study utilizes the large eddy simulation (LES) method to investigate the impact of oxygen mass fraction (14~23%) on the dynamic combustion characteristics of the non-premixed bluff body flame under the Mach numbers ranging from 0.1 to 0.2, encompassing the vortex structure, flame pattern, temperature pulsation, and pressure pulsation.

Large Eddy Simulation (LES) Model
LES directly computes large eddy and expresses small-scale eddy using a subgrid-scale (SGS) turbulence model [33,34].The filtered LES control equations are as follows: Continuous equation: ∂ρ ∂t Momentum equation: Energy equation: Component equation: Here, τ ij is the viscous flux, q j is the heat flux, h is the enthalpy flux, Y k is the convection component flux, τ ij is the stress tensor, and u ij is the viscosity coefficient.
The effect of small-scale eddies on the flow field can be calculated by adding an SGS stress term [35] to the LES.The SGS stress term can be characterised by assuming the eddy viscosity as follows: where µ sgs is the SGS viscosity coefficient, and k is the SGS turbulent kinetic energy, which is defined as In order to avoid the emergence of finely scaled turbulent structures in the near-wall region, the RANS solution model is used for the near-wall boundary layer [36,37].The formula used was as follows: (7) where κ = 0.41 and d w are the wall distances, S is the strain rate, y + is the dimensionless distance near the wall, C Smag = 0.2, and ∆ is the filter size.

Discrete Phase Model
The liquid fuel was used in this study.The numerical simulation must incorporate the mass, momentum, energy, and component transfers between phases.A discrete phase model simulated the force conditions, motion, and heat transfer of fuel.
The differential equation for the force on a droplet is where is the droplet drag force per unit mass, Re p = is the droplet Reynolds number, µ is the gas phase dynamic viscosity, ρ is the gas phase density, d p is the droplet diameter, → F represents other forms of force exerting on the droplets, and C D is the drag force coefficient.
The differential equation for the droplet trajectory is This equation can be simultaneously integrated along several coordinate directions during any time step ∆t to obtain the trajectory and velocity of the droplets.
The differential equation for droplets undergoing heating or cooling (without phase transition) is where T p and T vap represent the droplet temperature and evaporation temperature, respectively, and m p and m p,0 represent the droplet current and initial masses, respectively.During evaporation, the mass transfer equation is The mass consumption by the droplet is denoted as The evaporation heat transfer equation is In the equation above, dm p /dt represents the evaporation rate.The boiling evaporation rate is expressed as [38] where T ∞ and T p represent the continuous phase and liquid temperatures, respectively.

Combustion Model
A partially premixed combustion model was used to simulate the turbulent combustion in a bluff body.The mixing fraction f and process variable C were introduced as independent variables for the FGM combustion model [39].These are expressed as follows: Processes 2024, 12, 622 5 of 14 Here, Y indicates the mass fraction of the substance, M represents the relative atomic mass of the substance, and the subscripts 1 and 2 denote the fuel and oxidant inlets, respectively.Carbon and hydrogen are represented by the letters C and H, respectively.Only diffusion and convection affect the mixing fraction f.The transport equation is written as follows: The control equation for the process variable C is In the two equations above, D represents the laminar diffusion coefficient and • ω represents the reaction rate; they are expressed as follows: .
Here, Y i represents the component i, • ω Yi is the reaction rate of component i, H Yi is the enthalpy of component i at temperature T, c p is the specific heat at constant pressure of component i at temperature T, and n represents the total number of components.
The C 12 H 23 reaction mechanism described by Kundu et al. [40] was used to simulate the combustion reaction, replacing fuel in the model's combustion chamber to simulate combustion.

Numerical Configuration and Boundary Conditions 2.4.1. Model Structure
The model used for the calculations is the combustion chamber model developed by Roach et al. [41] The length of the computational domain is 1150 mm, the height is 100 mm, and the width is 50 mm.The V-shaped bluff body structure has a 30 • top angle and a 40 mm wide trailing edge groove (see Figure 1).

Meshing
The combustor computational domain is meshed to generate a structured mesh with a mesh size of 3.4 million, and the mesh shape is polyhedral.Periodic boundary conditions are applied to the front and rear faces of the combustor.The mesh size is densified at the position of the flame stabilizer and recirculation zone.The minimum mesh size of the flame stabilizer surface in the normal direction was 0.2 mm.The number of encrypted layers is set to 10 in order to fulfil the requirements of the wall function calculation (see Figure 2).

Meshing
The combustor computational domain is meshed to generate a structured mesh with a mesh size of 3.4 million, and the mesh shape is polyhedral.Periodic boundary conditions are applied to the front and rear faces of the combustor.The mesh size is densified at the position of the flame stabilizer and recirculation zone.The minimum mesh size of the flame stabilizer surface in the normal direction was 0.2 mm.The number of encrypted layers is set to 10 in order to fulfil the requirements of the wall function calculation (see Figure 2).a mesh size of 3.4 million, and the mesh shape is polyhedral.Periodic boundary conditions are applied to the front and rear faces of the combustor.The mesh size is densified at the position of the flame stabilizer and recirculation zone.The minimum mesh size of the flame stabilizer surface in the normal direction was 0.2 mm.The number of encrypted layers is set to 10 in order to fulfil the requirements of the wall function calculation (see Figure 2).To verify the mesh independence, the number of meshes was tested.The velocity time series in the core region of the recirculation zone behind the bluff body and in the shear layer were analysed spectrally.The spectral curves in Figure 3 reflect Kolmogorov's −5/3 power law, indicating that the setup mesh size meets the requirements of the LES simulation [42].Therefore, a mesh of 3.4 million nodes was chosen in this study for stable and unstable calculations.To verify the mesh independence, the number of meshes was tested.The velocity time series in the core region of the recirculation zone behind the bluff body and in the shear layer were analysed spectrally.The spectral curves in Figure 3 reflect Kolmogorov's −5/3 power law, indicating that the setup mesh size meets the requirements of the LES simulation [42].Therefore, a mesh of 3.4 million nodes was chosen in this study for stable and unstable calculations.
a mesh size of 3.4 million, and the mesh shape is polyhedral.Periodic boundary conditions are applied to the front and rear faces of the combustor.The mesh size is densified at the position of the flame stabilizer and recirculation zone.The minimum mesh size of the flame stabilizer surface in the normal direction was 0.2 mm.The number of encrypted layers is set to 10 in order to fulfil the requirements of the wall function calculation (see Figure 2).To verify the mesh independence, the number of meshes was tested.The velocity time series in the core region of the recirculation zone behind the bluff body and in the shear layer were analysed spectrally.The spectral curves in Figure 3 reflect Kolmogorov's −5/3 power law, indicating that the setup mesh size meets the requirements of the LES simulation [42].Therefore, a mesh of 3.4 million nodes was chosen in this study for stable and unstable calculations.

Model Verification
In order to confirm the accuracy of the model, Deng et al. [32] conducted model validation under the experimental conditions in reference [43], using the same numerical model used in this study.Figure 4 compares the calculated results and experimental data for different radial positions of reactive and nonreactive flows, and the results are consistent.

Model Verification
In order to confirm the accuracy of the model, Deng et al. [32] conducted model validation under the experimental conditions in reference [43], using the same numerical model used in this study.Figure 4   Comparison between numerical simulation and experimental results [32].Black dots/black square: experimental data; red line: numerical simulation results.

Vortex Structure
BVK instability is an absolute instability caused by asymmetric vortex shedding in the wake, and KH instability is a convective instability caused by symmetric vortex Comparison between numerical simulation and experimental results [32].Black dots/black square: experimental data; red line: numerical simulation results.

Vortex Structure
BVK instability is an absolute instability caused by asymmetric vortex shedding in the wake, and KH instability is a convective instability caused by symmetric vortex separation within the shear layer.The von Karman street is suppressed at high combustion product temperatures compared to reactant temperatures, resulting in a wake dominated by high-frequency KH instabilities.The mutual transformation between BVK instability and KH instability causes violent exothermic and large pressure gradient changes [44].When the incoming Mach number was 0.15, under the conditions of 23% and 20% O 2 mass fraction, the flow in the afterburner was characterised by BVK instability, an uneven velocity distribution in the flow field, and a large oscillation of the bluff wake.When the O 2 mass fraction was decreased from 20% to 17%, the number of vortex structures in the combustion chamber was reduced, the vortex tubes on the shear layer were symmetrically distributed and approximately linear, the velocity distribution was more uniform, the BVK instability was suppressed, and the overall flow was characterised by KH instability.When the oxygen content was reduced from 17% to 14%, the shear layer still showed KH instability.However, BVK instability started to appear downstream of the combustion chamber approximately 150 mm from the trailing edge of the bluff body.It was found that during the reduction in the oxygen content from 17% to 14%, the temperature strength was weakened, the expansion vortex and oblique pressure vortex decreased more significantly relative to the viscous vortex, the vortex in the shear layer of the bluff body increased, and the BVK instability was enhanced.At an incoming Mach number of 0.2, the flow exhibited BVK instability at 23% oxygen content, and the number of vortex structures in the stabilizer decreased when the O 2 mass fraction was reduced from 23% to 20%; the BVK instability was suppressed, and the vortex structures were distributed in a linear pattern.During the process of reduction in the oxygen content from 20% to 14%, the number of vortex structures inside the bluff increased.BVK instability occurred in the wake, and the area where the BVK instability occurred gradually moved upstream (see Figure 5).
icantly relative to the viscous vortex, the vortex in the shear layer of the bluff body increased, and the BVK instability was enhanced.At an incoming Mach number of 0.2, the flow exhibited BVK instability at 23% oxygen content, and the number of vortex structures in the stabilizer decreased when the O2 mass fraction was reduced from 23% to 20%; the BVK instability was suppressed, and the vortex structures were distributed in a linear pattern.During the process of reduction in the oxygen content from 20% to 14%, the number of vortex structures inside increased.BVK instability occurred in the wake, and the area where the BVK instability occurred gradually moved upstream (see Figure 5).

Flame Structure
The variations in OH concentration and temperature distribution with oxygen content in the combustion chamber at different incoming Mach numbers are shown in Figure 6 and Figure 7, respectively.The figures show that the flame pattern was essentially the same as the vortex structure pattern.At a constant incoming flow Mach number, the OH

Flame Structure
The variations in OH concentration and temperature distribution with oxygen content in the combustion chamber at different incoming Mach numbers are shown in Figures 6 and 7

Spectral Analysis
Temperature and pressure pulsations are commonly used to monitor combustion kinetics [45][46][47] and are analysed in this section for XO2 = 14-23%, Ma = 0.1-0.2, and p = 0.1 MPa.The time domain curves of temperature and pressure pulsations were analysed for pressure and temperature pulsation data at the monitoring points using the Fast Fourier Transform (FFT) shown in Equation ( 20) to obtain the frequency domain information.
Due to the significant impact of the separated free shear layer on dynamic combustion characteristics, the selected monitoring point is situated at 40 mm downstream of the When the incoming Mach number is 0.1, under the condition of 23% oxygen content, the flame of the shear layer of the flame stabilizer is dominated by KH instability, and the flame shows a straight line.At approximately 100 mm from the downstream of the flame stabilizer, the flame is severely stretched and almost fractured, and the flame stability is poor.As the oxygen content decreased, the combustion reaction rate and flame temperature decreased significantly.When the oxygen content was 20%, there were two dynamic characteristics of KH instability and BVK instability in the flow field, and the stabilizer near-field flame still showed obvious convective instability characteristics: a symmetrical distribution of the flame in the shear layer, and a wavy shape at the interface between the flame and the external gas mixture.However, with an increase in the axial distance, the flow field is transformed into absolute instability, and the flame in the wake stream oscillates up and down, showing an S-shaped flame.When the O 2 mass fraction was decreased from 20% to 17%, the flame was more zigzagged and almost broken at a position approximately 100 mm from the downstream of the bluff body, and the flame stability was further decreased.When the O 2 mass fraction was decreased to 14%, the area of BVK instability moved upstream, but the flame continuity was simultaneously enhanced.When the incoming Mach number was 0.15, the flame after the stabilizer was mainly composed of swirling flames under 23% oxygen content, showing an evident vortex structure.When the O 2 mass fraction was decreased from 23% to 17%, the BVK instability start position moved downstream, and the blunt body near-field flame was linear.When the O 2 mass fraction further decreased to 14%, the BVK instability was enhanced, and the starting position of the flame oscillation moved upstream.When the incoming Mach number was 0.2, the flame after the stabilizer was mainly dominated by BVK instability under 23% oxygen content.When the O 2 mass fraction was decreased from 23% to 20%, the BVK instability was suppressed, the flame was approximately linear and smoother, and the flame stability was improved.As the O 2 mass fraction was further decreased, the BVK instability began to appear downstream of the combustion chamber, and the flame approached the extinguishing state when the O 2 mass fraction was reduced to 14%.
At an oxygen content of 23%, the OH concentration and temperature in the bluff body decreased, and the main reaction zone moved downstream as the incoming Mach number increased.When the incoming Mach number was increased from 0.1 to 0.15, the near-field flame of the bluff body exhibited a large oscillation.In contrast, the BVK instability start position moved downstream as the incoming Mach number was increased to 0.2.It was found that with an increase in airflow velocity, owing to the relatively obvious increase in the viscous vortex volume, the vortex volume of the shear layer of the bluff body increases, which increases the BVK instability.In contrast, in the process of increasing the incoming Mach number from 0.15 to 0.2, owing to the weakening of the temperature strength in the bluff body, the flame temperature decreases substantially, and the BVK instability is relatively weakened.Under the conditions of 20% and 17% oxygen content, with an increase in the incoming Mach number, the OH concentration and temperature in the bluff body showed an overall increasing trend, the BVK instability was suppressed, and the flame oscillation was reduced.Under the 14% oxygen content condition, when the incoming Mach number was increased from 0.1 to 0.15, the flame temperature inside the passivate increased, and the BVK instability was suppressed.However, when the incoming Mach number was increased to 0.2, the BVK instability was enhanced.It was found that under the condition of 14% oxygen content, as the incoming Mach number was increased to 0.2, the flame was close to extinction, the expansion vortex and oblique pressure vortex decreased significantly relative to the viscous vortex, the vortex of the shear layer of the passivate body increased, and the BVK instability was enhanced.

Spectral Analysis
Temperature and pressure pulsations are commonly used to monitor combustion kinetics [45][46][47] and are analysed in this section for X O2 = 14-23%, Ma = 0.1-0.2, and p = 0.1 MPa.The time domain curves of temperature and pressure pulsations were analysed for pressure and temperature pulsation data at the monitoring points using the Fast Fourier Transform (FFT) shown in Equation (20) to obtain the frequency domain information.
Due to the significant impact of the separated free shear layer on dynamic combustion characteristics, the selected monitoring point is situated at 40 mm downstream of the trailing edge of the bluff body, positioned within the recirculation zone and the mainstream area of the interface, where the pressure pulsations reach their peak intensity and exhibit the greatest amplitude.
The temperature and pressure pulsation spectra are shown in Figures 8 and 9, respectively.Upon varying the incoming flow velocity and oxygen content, the reactive flow formed by the fluid as it passed through the bluff body switched between KH and BVK instabilities.When the BVK instability dominates the flow field, the pressure and temperature pulsations exhibit consistent dominant frequencies in the 750 to 1600 Hz range.This is due to the high-temperature combustion of the bluff body flame coupled with fluid dynamics.The temperature pulsation amplitude increases with increasing temperature intensity and BVK instability, and the pressure pulsation is mainly affected by the dynamic properties of the flow field behind the bluff body.The pressure pulsation amplitude increases with increasing BVK instability.
In the case of the KH instability-dominated flow field, the main frequency and amplitude of the pressure pulsation decreased, and the main frequency was in the range of 100-300 Hz.As the fluctuations in the flow field dominated by the KH were significantly reduced and only smaller fluctuations were generated in the downstream region, combustion was more stable and combustion exotherm was more uniform within the shear layer than in the flow field, which was mainly affected by the BVK instability.Consequently, the temperature pulsations had no obvious dominant frequency, and the frequency distribution in the spectrum was highly uniform.
In the case of the KH instability-dominated flow field, the main frequency and amplitude of the pressure pulsation decreased, and the main frequency was in the range of 100-300 Hz.As the fluctuations in the flow field dominated by the KH instability were significantly reduced and only smaller fluctuations were generated in the downstream region, combustion was more stable and combustion exotherm was more uniform within the shear layer than in the flow field, which was mainly affected by the BVK instability.Consequently, the temperature pulsations had no obvious dominant frequency, and the frequency distribution in the spectrum was highly uniform.

Conclusions
This paper investigates the vortex structure, flame morphology, temperature pulsation, and pressure pulsation of a bluff body flame in the range of 14-23% oxygen content and 0.1-0.2Ma using LES numerical simulation.The results revealed that with changes in the O2 mass fraction and Mach number, two distinct instability characteristics were observed behind the bluff body: BVK instability and KH instability.As the O2 mass fraction reduced, the flame temperature intensity decreased significantly, with a decrease in the OH concentration and temperature.As the inlet velocity increases, the viscous vortex vol-

Conclusions
This paper investigates the vortex structure, flame morphology, temperature pulsation, and pressure pulsation of a bluff body flame in the range of 14-23% oxygen content and 0.1-0.2Ma using LES numerical simulation.The results revealed that with changes in the O 2 mass fraction and Mach number, two distinct instability characteristics were observed behind the bluff body: BVK instability and KH instability.As the O 2 mass fraction reduced, the flame temperature intensity decreased significantly, with a decrease in the OH concentration and temperature.As the inlet velocity increases, the viscous vortex volume increases and the vortex volume increases in the bluff body shear layer, which may lead to enhanced BVK instability.When the oxygen content decreases, the number of vortex structures decreases, the vortex tubes in the shear layer are symmetrically distributed, the velocity distribution is more uniform, and the KH instability is enhanced.However, as the oxygen content decreases, the temperature will decrease, and the expansion vortex and oblique pressure vortex decrease significantly, leading to an increase in the number of vortices in the shear layer and an enhancement of the BVK instability.The variation between BVK instability and KH instability is nonlinear with changes in oxygen content and inlet velocity under the coupled influence of inlet velocity and oxygen content.At a low Mach number (Ma = 0.1), the separated shear layer transitions from KH to BVK instability as the oxygen content decreases.As Ma ≥ 0.2, the flame and flow characteristics downstream from the flame stabilizer gradually shift from BVK instability to KH instability.In particular, the bluff body flame approached extinction at Ma = 0.2, and the oxygen content was 14%, indicating a significant challenge in achieving efficient and stable combustion in real aero-engines under high-speed and oxygen-lean conditions.Additionally, under the domination of the BVK instability, the pressure and temperature pulsation exhibited a consistent dominant frequency in the 750 to 1600 Hz range.Under the KH instability dominance, the temperature pulsation had no significant dominant frequency, and the pressure pulsation had a lower frequency of about 100-300 Hz.

Figure 2 .
Figure 2. Schematic diagram of the grid used.

Figure 2 .
Figure 2. Schematic diagram of the grid used.

Figure 2 .
Figure 2. Schematic diagram of the grid used.

Figure 3 .
Figure 3. Verification of mesh independence.(a) Location of the shear layer.(b) Core location of the recirculation zone.
compares the calculated results and experimental data for different radial positions of reactive and nonreactive flows, and the results are consistent.

Figure 4 .
Figure 4. Comparison between numerical simulation and experimental results[32].Black dots/black square: experimental data; red line: numerical simulation results.

(Figure 5 .
Figure 5. Influence of O2 mass fraction on vortex structure shedding behind bluff body under different velocities.

Figure 5 .
Figure 5. Influence of O 2 mass fraction on vortex structure shedding behind bluff body under different velocities.

Figure 6 .Figure 6 .
Figure 6.Influence of O2 mass fraction on OH concentration field under different velocities.

Figure 6 .Figure 7 .
Figure 6.Influence of O2 mass fraction on OH concentration field under different velocities.

Figure 7 .
Figure 7. Influence of O 2 mass fraction on temperature field under different velocities.

Figure 8 .
Figure 8. Changes in temperature pulsation with O2 mass fraction under various inflow velocities.

Figure 8 .
Figure 8. Changes in temperature pulsation with O 2 mass fraction under various inflow velocities.

Figure 9 .
Figure 9. Changes in pressure pulsation with O2 mass fraction under various inflow velocities.

Figure 9 .
Figure 9. Changes in pressure pulsation with O 2 mass fraction under various inflow velocities.

Table 1 .
2.4.3.Boundary ConditionsANSYS FLUENT 14.0 was used for numerical simulations of the model combustor.Control variables were used to investigate the influence of O 2 mass fraction on the performance of the bluff body at different inlet Mach numbers.The operating conditions are listed in Table1.Operating conditions to study the effect of O 2 mass fraction on the combustion performance of the bluff body.