Direct numerical simulation and analysis of a hydrogen/air swirling premixed flame in a micro combustor

A three dimensional spatially developing hydrogen/air premixed flame in a micro combustor with a moderate Reynolds number and a high swirl number is studied using direct numerical simulation. The inflow mixture is composed of hydrogen and air at an equivalent ratio of 1.0 in the jet core region, and pure air elsewhere. The maximum axial velocity at the inlet is 100 m/s. A fourth-order explicit Runge–Kutta method for time integration and an eighth-order central differencing scheme for spatial discretization are used to solve the full Navier–Stokes (N–S) equation system. A 9 species 19-step reduced mechanism for hydrogen/air combustion is adopted. Vortex and turbulence characteristics are examined. Two instabilities, namely Kalvin–Helmholtz instability and centrifugal instability, are responsible for the transition from laminar flow to turbulence. A cone-like vortex breakdown is observed both in the isothermal swirling flow and in the swirling flame. One dimensional premixed laminar flame is studied, the structure of which is compared with that of the multi-dimensional one. Probability density functions of the curvature and tangential strain rate are presented. It is shown that the flame curvature has a near zero mean, and the flame aligns preferentially with extensive strain. Finally, the turbulent premixed flame regime diagram is used to characterize the flame. It is found that most of the flame elements lie in the laminar flame regime and the thin reaction zones regime.     

Measurement of the instantaneous flame front structure of syngas turbulent premixed flames at high pressure

Instantaneous flame front structure of syngas turbulent premixed flames including the local radius of curvature, the characteristic radius of curvature, the fractal inner cutoff scale and the local flame angle were derived from the experimental OH-PLIF images. The CO/H₂/CO₂/air flames as a model of syngas/air combustion were investigated at pressure of 0.5 MPa and compared to that of CH₄/air flames. The convex and concave structures of the flame front were detected and statistical analysis including the PDF and ADF of the local radius of curvature and local flame angle were conducted. Results show that the flame front of turbulent premixed flames at high pressure is a wrinkled flame front with small scale convex and concave structures superimposed with large scale flame branches. The convex structures are much more frequent than the concave ones on flame front which reflects a general characteristic of the turbulent premixed flames at high pressure. The syngas flames possess much wrinkled flame front with much smaller fine cusps structure compared to that of CH₄/air flames and the main difference is on the convex structure. The effect of turbulence on the general wrinkled scale of flame front is much weaker than that of the smallest wrinkled scale. The general wrinkled scale is mainly dominated by the turbulence vortex scale, while, the smallest wrinkled scale is strongly affected by the flame intrinsic instability. The effect of flame intrinsic instability on flame front of turbulent premixed flame is mainly on the formation of a large number of convex structure propagating to the unburned reactants and enlarge the effective contact surface between flame front and unburned reactants.     

Flame speed and tangential strain measurements in widely stratified partially premixed flames interacting with grid turbulence

Methane-air partially premixed flames subjected to grid-generated turbulence are stabilized in a two-slot burner with initial fuel concentration differences leading to stratification across the stoichiometric concentration. The fuel concentration gradient at the location corresponding to the flame base is measured using planar laser induced fluorescence (PLIF) of acetone in the non-reacting mixing field. Simultaneous PLIF of the OH radical and particle image velocimetry (PIV) measurements are performed to deduce the flow velocity and the flame front. These flames exhibit a convex premixed flame front and a trailing diffusion flame, with flow divergence upstream of the flame, as indicated by the instantaneous OH–PLIF, Mie scattering images, and PIV data. The mean streamwise velocity profile attains a global minimum just upstream of the flame front due to expansion of a gases caused by heat release. The flame speed measured just upstream of the flame leading edge is normalized with respect to the turbulent stoichiometric flame speed that takes into account variations in turbulent intensity and integral length scale. The turbulent edge flame speed exceeds the corresponding stoichiometric premixed flame speed and reaches a peak at a certain concentration gradient. The mean tangential strain at the flame leading edge locally peaks at the concentration gradient corresponding to the peak flame speed. The strain varies non-monotonically with the flame curvature unlike in a non-stratified curved premixed flame. The mechanism of peak flame speed is explained as the competition between availability of hot excess reactants from the premixed flame branches to the flame stretch induced due to flame curvature. The results suggest that the stabilization of lifted turbulent partially premixed flames occurs through an edge flame even at a relatively gentle concentration gradient. The strain is also evaluated along the flame front; it peaks at the flame leading edge and decreases gradually on either side of the leading edge. The present results also show qualitatively similar trends as those of laminar triple flames.     

Experimental investigation of three-dimensional flame-front structure in premixed turbulent combustion—I: hydrocarbon/air bunsen flames

OH concentrations and three-dimensional gradients of the reaction progress variable have been measured in turbulent liquefied petroleum gas/air and compressed natural gas/air premixed flames stabilized on a Bunsen-type burner with a combined two-sheet Rayleigh scattering and planar LIF-OH imaging technique. The progress variable is observed to undergo a transition from lamella-like to non-flamelet front structure with increasing turbulence. This is consistent with the recently proposed change of the combustion regimes from complex-strain to turbulent flame front regime on a recently proposed premixed combustion diagram. The anisotropy of local flame-front orientation in three-dimensional space is explained by the forward propagation ability of the planar turbulent flame brush. Weighting functions have thus been derived for the isotropic pdf distributions of the in-plane and out-of-plane orientation angles to agree better with the experimental data. A linear scaling is found between the overall flame surface area and the turbulence intensity normalized by the laminar burning velocity. However, flames with excess backward-facing flame fronts do not comply with this linear relationship, showing enhanced flame surface folding. The thin-flame assumption breaks down when non-flamelet broadening effects become important, although the pdf’s of the progress variable are still bimodal-like. Non-unity Lewis-number combined curvature effects are evident for LPG/air flames of weak turbulence, in that the conditional mean scalar dissipation increases steadily from the unburnt to burnt side across the flame brush. A consistent correlation exists between the Favre-averaged scalar dissipation and progress variable variance. This implies that small-scale scalar dissipation of local flame-fronts is linked to large-scale scalar fluctuations. Sub- or super-flamelet OH concentration is found in lean LPG/air or CNG/air premixed flames, respectively, and occurs in line with a positive or negative correlation between OH concentrations and magnitudes of the progress variable gradient.     

The effects of hydrogen substitution on turbulent premixed methane–air kernels using direct numerical simulation

Direct numerical simulations (DNS) are used to assess the effects of hydrogen substitution on lean premixed methane–air kernels during the early stages of growth in freely decaying turbulence. Two-dimensional simulations with a detailed 68-step reaction mechanism are carried out at equivalence ratios of ? = 0.53 and ? = 0.625, both with and without the substitution of a 30% fuel mole fraction of hydrogen. A comparative analysis is made into the changes in turbulent flame speeds, global stretch rates, and flame wrinkling at different turbulence intensities. The underlying causes of these changes are investigated through the distributions of the surface-conditioned displacement speed, strain rate and curvature. Direct comparison is made with the planar flame results of Hawkes and Chen [1] to assess the qualitative effects of kernel geometry in combination with a hydrogen-enriched fuel. It was found that the reduced effective fuel Lewis number and preferential diffusion of hydrogen, combined with the higher stretch rates and mean positive curvature of the kernel make the effects of hydrogen enrichment much more pronounced in kernels compared to planar flames.     

Investigations of swirl flames in a gas turbine model combustor: II. Turbulence–chemistry interactions

The thermochemical states of three swirling CH₄/air diffusion flames, stabilized in a gas turbine model combustor, were investigated using laser Raman scattering. The flames were operated at different thermal powers and air/fuel ratios and exhibited different flame behavior with respect to flame instabilities. They had previously been characterized with respect to their flame structures, velocity fields, and mean values of temperature, major species concentrations, and mixture fraction. The single-pulse multispecies measurements presented in this article revealed very rapid mixing of fuel and air, accompanied by strong effects of turbulence–chemistry interactions in the form of local flame extinction and ignition delay. Flame stabilization is accomplished mainly by hot and relatively fuel-rich combustion products, which are transported back to the flame root within an inner recirculation zone. The flames are not attached to the fuel nozzle, and are stabilized approximately 10 mm above the fuel nozzle, where fuel and air are partially premixed before ignition. The mixing and reaction progress in this area are discussed in detail. The flames are short (<50 mm), especially that exhibiting thermoacoustic oscillations, and reach a thermochemical state close to adiabatic equilibrium at the flame tip. The main goals of this article are to outline results that yield deeper insight into the combustion of gas turbine flames and to establish an experimental database for the validation of numerical models.     

From diffusion to premixed flames in an H2/air opposed-jet burner: the role of edge flames

Edge flames obtained on a hydrogen/air non-premixed opposed-jet burner after the local extinction of the disk-shaped diffusion flame are investigated with 2-D direct numerical simulations using detailed chemical kinetics and transport. Over a large range of flowrates, edge flames were found to coexist with the well-known strongly burning diffusion flames corresponding to the upper branch of the S-shaped curve. The critical flowrates of the strong hysteresis associated with the transitions between the two solution branches were identified: re-establishment of the diffusion flame is controlled by the propagation of the edge flame and cannot be represented simply by the extinction scalar dissipation rate. It was also observed that in all the flow conditions simulated, the edge flame was able to consume all the supplied fuel by re-orienting itself, varying its flame surface area, or changing its structure. The latter was found to depend on the flow conditions (which strongly affects the degree of mixing ahead of the edge flame) and can take on different configurations ranging from a triple flame to an essentially premixed flame. Because of flame curvature and the preferential diffusion of hydrogen, the propagation speed of the edge flames was found to be higher than that of the corresponding planar premixed flames.     

Pressure influence on the flame front curvature of turbulent premixed flames: comparison between experiment and theory

In gas turbines, lean premixed combustion is executed in strongly turbulent flow fields and under high-pressure to allow large thermal loads within small-size combustors. Previous research on turbulent premixed flames has revealed the vital importance of flame-vortex interactions, but most of these investigations have been performed only at atmospheric pressure disregarding the large pressure dependency of the flame front dynamics. We report about spatially high-resolved laser-induced predissociation fluorescence imaging of OH (OH-LIPF) in premixed, high-pressure bluff-body stabilized methane/air flames. For each of the two measurement series with different equivalence ratio (φ = 0.7 and φ = 1.0), the planar flame topology at different pressures (0.1 to 1.1 MPa) but constant exit velocity was detected and stored for analysis. As the pressure was increased, the flame front contour of both equivalence ratios became strongly wrinkled with formation of highly curved flame front elements. For quantification of this phenomenon, the probability density function of flame curvature was evaluated with definition of the mean curvature radius as representative folding scale. To discuss different mechanisms of flame front disturbances according to their relevance, the flame curvature is compared with characteristic turbulence scales of the flow field and with the expected folding scale derived with Sivashinsky‘s formulation of linear flame instability theory. Significant changes become obvious especially if the pressure is increased up to 0.5 MPa. The mean curvature radius decreases distinctly and can be linked to the decreasing size of the Taylor length. Additionally, the formation of highly convoluted flame front elements is enforced by the increasing flame instability behavior. As the results show, the flame stoichiometry has a strong impact on the flame front topology at increasing pressures due to the differences of their flame dynamics.     

Preferential diffusion effects on the burning rate of interacting turbulent premixed hydrogen-air flames

The upstream interaction of twin premixed hydrogen-air flames in 2-D turbulence is studied using direct numerical simulations with detailed chemistry. The primary objective is to determine the effect of flame stretch on the overall burning rate during various stages of the interaction. Preferential diffusion effects are accounted for by varying the equivalence ratio from symmetric rich-rich to lean-lean interactions. The results show that the local flame front response to turbulence is consistent with previous understanding of laminar premixed flames, in that rich premixed flames become intensified in regions of negative strain or curvature, while the opposite response is found for lean premixed flames. The overall burning rate history with respect to the surface density variation is found to depend on the mixture condition; the consumption rate enhancement advances (follows) the surface enhancement for the rich-rich (lean-lean) case. For the lean-lean case, a self-turbulization mechanism results in a large positive skewness in the area-weighted mean tangential strain statistics. Because of the statistical dominance of positive stretch on the flame surface, the lean-lean case results in a significantly larger burning enhancement (over a twofold increase) in addition to the surface density production. For the case of rich-rich interaction, the abundance in hydrogen species results in an instantaneous overshoot of the radical pool in the post-flame region, resulting in an additional “burst” in the reactant consumption rate history, suggesting its potential impact on the pollutant formation process.     

Tubular premixed and diffusion flames: Effect of stretch and curvature

Tubular flames are ideal for the study of stretch and curvature effects on flame structure, extinction, and instabilities. Tubular flames have uniform stretch and curvature and each parameter can be varied independently. Curvature strengthens or weakens preferential diffusion effects on the tubular flame and the strengthening or weakening is proportional to the ratio of the flame thickness to the flame radius. Premixed flames can be studied in the standard tubular burner where a single premixed gas stream flows radially inward to the cylindrical flame surface and products exit as opposed jets. Premixed, diffusion and partially premixed flames can be studied in the opposed tubular flame where opposed radial flows meet at a cylindrical stagnation surface and products exit as opposed jets. The tubular flame flow configurations can be mathematically reduced to a two-point boundary value solution along the single radial coordinate. Non-intrusive measurements of temperature and major species concentrations have been made with laser-induced Raman scattering in an optically accessible tubular burner for both premixed and diffusion flames. The laser measurements of the flame structure are in good agreement with numerical simulations of the tubular flame. Due to the strong enhancement of preferential diffusion effects in tubular flames, the theory-data comparison can be very sensitive to the molecular transport model and the chemical kinetic mechanism. The strengthening or weakening of the tubular flame with curvature can increase or decrease the extinction strain rate of tubular flames. For lean H₂-air mixtures, the tubular flame can have an extinction strain rate many times higher than the corresponding opposed jet flame. More complex cellular tubular flames with highly curved flame cells surrounded by local extinction can be formed under both premixed and non-premixed conditions. In the hydrogen fueled premixed tubular flames, thermal-diffusive flame instabilities result in the formation of a uniform symmetric petal flames far from extinction. In opposed-flow tubular diffusion flames, thermal-diffusive flame instabilities result in cellular flames very close to extinction. Both of these flames are candidates for further study of flame curvature and extinction.     

Dynamics of premixed hydrogen/air flames in microchannels

The stabilization and dynamics of lean (φ=0.5) premixed hydrogen/air atmospheric-pressure flames in planar microchannels of prescribed wall temperature are investigated with respect to the inflow velocity and channel height (0.3 to 1.0 mm) using direct numerical simulation with detailed chemistry and transport. Rich dynamics starting from periodic ignition and extinction of the flame and further transitioning to symmetric V-shaped flames, asymmetric flames, oscillating and pulsating flames, and finally again to asymmetric flames are observed as the inlet velocity is increased. The richest behavior is observed for the 1.0-mm-height channel. For narrower channels, some of the dynamics are suppressed. The asymmetric flames, in particular, vanish for channel heights roughly less than twice the laminar flame thickness. Stability maps delineating the regions of the different flame types in the inlet velocity/channel height parameter space are constructed.     

Large-eddy simulation of lean hydrogen–methane turbulent premixed flames in the methane-dominated regime

The application of large-eddy simulation (LES) to the prediction of H₂-enriched lean methane–air turbulent premixed combustion is considered. A presumed conditional moment (PCM) subfilter-scale combustion model is coupled with the flame prolongation of intrinsic low-dimensional manifold (FPI) chemistry tabulation technique. The LES and PCM-FPI modelling procedures are then applied to the prediction of laboratory-scale axisymmetric Bunsen-type turbulent premixed flames. Both premixed methane–air and H₂-enriched methane–air flames are considered and the predicted solutions are examined and compared to available experimental data. The enriched flame has 20% H₂ in terms of mole fraction and lies in the methane-dominated regime of hydrogen–methane mixtures. The LES simulations predict similar qualitative trends to those found in the experiments for flame height and curvature. The addition of H₂ decreases the flame height and broadens the curvature probability density functions, which show a Gaussian-type shape centred around zero. Moreover, the enriched flame displays a higher degree of wrinkling with sharper ridges of negative curvature and larger pockets of positive curvature. Overall, the proposed treatment for the PCM-FPI combustion model, in terms of progress variable and tabulated data, seems to perform well for the H₂-enriched methane flame in the methane-dominated regime.     

Effects of CO addition on the characteristics of laminar premixed CH4/air opposed-jet flames

The effects of CO addition on the characteristics of premixed CH₄/air opposed-jet flames are investigated experimentally and numerically. Experimental measurements and numerical simulations of the flame front position, temperature, and velocity are performed in stoichiometric CH₄/CO/air opposed-jet flames with various CO contents in the fuel. Thermocouple is used for the determination of flame temperature, velocity measurement is made using particle image velocimetry (PIV), and the flame front position is measured by direct photograph as well as with laser-induced predissociative fluorescence (LIPF) of OH imaging techniques. The laminar burning velocity is calculated using the PREMIX code of Chemkin collection 3.5. The flame structures of the premixed stoichiometric CH₄/CO/air opposed-jet flames are simulated using the OPPDIF package with GRI-Mech 3.0 chemical kinetic mechanisms and detailed transport properties. The measured flame front position, temperature, and velocity of the stoichiometric CH₄/CO/air flames are closely predicted by the numerical calculations. Detailed analysis of the calculated chemical kinetic structures reveals that as the CO content in the fuel is increased from 0% to 80%, CO oxidation (R99) increases significantly and contributes to a significant level of heat-release rate. It is also shown that the laminar burning velocity reaches a maximum value (57.5 cm/s) at the condition of 80% of CO in the fuel. Based on the results of sensitivity analysis, the chemistry of CO consumption shifts to the dry oxidation kinetics when CO content is further increased over 80%. Comparison between the results of computed laminar burning velocity, flame temperature, CO consumption rate, and sensitivity analysis reveals that the effect of CO addition on the laminar burning velocity of the stoichiometric CH₄/CO/air flames is due mostly to the transition of the dominant chemical kinetic steps.     

Effect of differential diffusion on turbulent lean premixed hydrogen enriched flames through structure analysis

This study presents the flame structure influenced by the differential diffusion effects and evaluates the structural modifications induced by the turbulence, thus to understand the coupling effects of the diffusively unstable flame fronts and the turbulence distortion. Lean premixed CH₄/H₂/air flames were conducted using a piloted Bunsen burner. Three hydrogen fractions of 0, 30% and 60% were adopted and the laminar flame speed was kept constant. The turbulence was generated with a single-layer perforated plate, which was combined with different bulk velocities to obtain varied turbulence intensities. Quasi-laminar flames without the plate were also performed. Explicit flame morphology was obtained using the OH-PLIF. The curvature, flame surface density and turbulent burning velocity were measured. Results show that the preferential transport of hydrogen produces negatively curved cusps flanked with positively curved bulges, which are featured by skewed curvature pdfs and consistent with the typical structure caused by the Darrieus-Landau instability. Prevalent bulge-cusp like wrinkles remain with relatively weak turbulence. However, stronger turbulence can break the bulges to be finer, and induce random positively curved cusps, therefore to destroy the bulge-cusp structures. Evident positive curvatures are generated in this process modifying the skewed curvature pdfs to be more symmetric, while the negative curvatures are not affected seriously. From low to high turbulence intensities, the hydrogen addition always strengthens the flame wrinkling. The augmentation of flame surface density and turbulent burning velocity with hydrogen is even more obvious at higher turbulence intensity. It is suggested that the differential diffusion can persist and even be strengthened with strong turbulence.     

DNS of spark ignition and edge flame propagation in turbulent droplet-laden mixing layers

A parametric study of forced ignition at the mixing layer between air and air carrying fine monosized fuel droplets is done through one-step chemistry direct numerical simulations to determine the influence of the size and volatility of the droplets, the spark location, the droplet-air mixing layer initial thickness and the turbulence intensity on the ignition success and the subsequent flame propagation. The propagation is analyzed in terms of edge flame displacement speed, which has not been studied before for turbulent edge spray flames. Spark ignition successfully resulted in a tribrachial flame if enough fuel vapour was available at the spark location, which occurred when the local droplet number density was high. Ignition was achieved even when the spark was offset from the spray, on the air side, due to the diffusion of heat from the spark, provided droplets evaporated rapidly. Large kernels were obtained by sparking close to the spray, since fuel was more readily available. At long times after the spark, for all flames studied, the probability density function of the displacement speed was wide, with a mean value in the range 0.55-0.75S_L, with S_L the laminar burning velocity of a stoichiometric gaseous premixed flame. This value is close to the mean displacement speed in turbulent edge flames with gaseous fuel. The displacement speed was negatively correlated with curvature. The detrimental effect of curvature was attenuated with a large initial kernel and by increasing the thickness of the mixing layer. The mixing layer was thicker when evaporation was slow and the turbulence intensity higher. However, high turbulence intensity also distorted the kernel which could lead to high values of curvature. The edge flame reaction component increased when the maximum temperature coincided with the stoichiometric contour. The results are consistent with the limited available experimental evidence and provide insights into the processes associated with ignition of practical spray flames.     

Characterization of the effects of hydrogen addition in premixed methane/air flames

The use of hydrogenated fuels shows considerable promise for applications in gas turbines and internal combustion engines. In the present work, the effects of hydrogen addition in methane/air flames are investigated using both a laminar flame propagation facility and a high-pressure turbulent flame facility. The aim of this research is to contribute to the characterization of lean methane/hydrogen/air premixed turbulent flames at high pressures, by studying the flame front geometry, the flame surface density and the instantaneous flame front thermal thickness distributions. The experiments and analyses show that a small amount of hydrogen addition in turbulent premixed methane–air flames introduces changes in both instantaneous and average flame characteristics.     

LES of a premixed jet flame DNS using a strained flamelet model

Turbulent premixed flames in the thin and broken reaction zones regimes are difficult to model with Large Eddy Simulation (LES) because turbulence strongly perturbs subfilter scale flame structures. This study addresses the difficulty by proposing a strained flamelet model for LES of high Karlovitz number flames. The proposed model extends a previously developed premixed flamelet approach to account for turbulence’s perturbation of subfilter premixed flame structures. The model describes combustion processes by solving strained premixed flamelets, tabulating the results in terms of a progress variable and a hydrogen radical, and invoking a presumed PDF framework to account for subfilter physics. The model is validated using two dimensional laminar flame studies, and is then tested by performing an LES of a premixed slot-jet direct numerical simulation (DNS). In the premixed regime diagram this slot-jet is found at the edge of the broken reaction zones regime. Comparisons of the DNS, the strained flamelet model LES, and an unstrained flamelet model LES confirm that turbulence perturbs flame structure to leading order effect, and that the use of an unstrained flamelet LES model under-predicts flame height. It is shown that the strained flamelet model captures the physics characterizing interactions of mixing and chemistry in highly turbulent regimes.     

LES–ODT study of turbulent premixed interacting flames

The LES–ODT model is implemented for the study of twin turbulent premixed flames in decaying isotropic turbulence. The approach is based on the coupling of large-eddy simulation (LES) for mass and momentum with a fixed 3D lattice of 1D fine-grained solutions based on the one-dimensional turbulence (ODT) model. The ODT solutions for momentum and reactive scalars are designed to capture subgrid scale physics that is not captured by LES. The LES–ODT formulation is capable of capturing important fine-scale processes, such as flame–flame interactions, which play an important role in flame shortening in turbulent premixed flames, and the role of preferential diffusion on curved flames’ structures.     

LES flamelet modeling of hydrogen combustion considering preferential diffusion effect

A flamelet-generated manifold (FGM) method that explicitly considers the preferential diffusion effect, referred to as FGM-PD method, is employed for large-eddy simulations (LESs) of a lean-premixed H₂/air low-swirl lifted flame, and the validity is examined by comparing with the experiment. First, the applicability of the FGM-PD method is investigated by one-dimensional numerical simulations of planar laminar premixed H₂/air flames. Next, LESs of a lean-premixed H₂/air low-swirl lifted flame are performed employing the FGM-PD and conventional FGM methods. Results of the one-dimensional numerical simulations show the importance of considering preferential diffusion to accurately predict species concentrations near the flame front. The FGM-PD method accurately predicts this, and therefore, reproduces the laminar burning velocity and spatial distributions of temperature and mixture fraction. Three-dimensional LES results confirm that the prediction accuracy of the velocities near the flame front is improved by employing this FGM-PD method. Additionally, the OH mass fraction distribution predicted by the FGM-PD method exhibits the inhomogeneous finger-like structure, which has been observed in previous experiments. This inhomogeneity of OH mass fraction distribution, which corresponds to that of the reaction rate, predicted by the FGM-PD method, strongly affects the flame front structure.     

Effects of the external turbulence on centrally-ignited spherical unstable Ch4/H2/air flames in the constant-volume combustion bomb

In the present study, we conducted experiments to investigate the effects of external turbulence on the development of spherical H₂/CH₄/air unstable flames developments at two different equivalence ratios associated with different turbulent intensities using a spherical constant-volume turbulent combustion bomb and high speed schlieren photography technology. Flame front morphology and acceleration process were recorded and different effects of weak external turbulent flow field and intrinsic flame instability on the unstable flame propagation were compared. Results showed the external turbulence has a great influence on the unstable flame propagation under rich fuel conditions. For fuel-lean premixed flames, however, the effects of external turbulence on the morphology of the cellular structure on the flame front was not that obvious. Critical radius decreased firstly and then kept almost unchanged with the augment of the turbulence intensity. This indicated the dominating inhibiting effect of flame stretch on the turbulent premixed flame at the initial stage of the flame front development. Beyond the critical radius, the acceleration exponent was found increasing with the enhancement of initial turbulence intensity for fuel-lean premixed flames. For fuel-rich conditions, however, the initial turbulence intensity had little effect on acceleration exponent. In order to evaluate the important impact of the intrinsic flame instability and external turbulent flow field for spherical propagating premixed flames, intrinsic flame instability scale and average diameter of vortex tube were calculated. Intrinsic flame instability scale decreased greatly and then stayed unchanged with the propagation of the flame front. The comparison between intrinsic flame instability scale and average diameter of vortex tube demonstrated that the external turbulent flow filed will be more important for the evolution of wrinkle structure in the final stage of the flame propagation, when the turbulence intensity was more than 0.404 m/s.