skip to main content


Title: Reynolds number scaling of burning rates in spherical turbulent premixed flames
In the flamelet regime of turbulent premixed combustion the enhancement in the burning rates originates primarily from surface wrinkling. In this work we investigate the Reynolds number dependence of burning rates of spherical turbulent premixed methane/air flames in decaying isotropic turbulence with direct numerical simulations. Several simulations are performed by varying the Reynolds number, while keeping the Karlovitz number the same, and the temporal evolution of the flame surface is compared across cases by combining the probability density function of the radial distance of the flame surface from the origin with the surface density function formalism. Because the mean area of the wrinkled flame surface normalized by the area of a sphere with radius equal to the mean flame radius is proportional to the product of the turbulent flame brush thickness and peak surface density within the brush, the temporal evolution of the brush and peak surface density are investigated separately. The brush thickness is shown to scale with the integral scale of the flow, evolving due to decaying velocity fluctuations and stretch. When normalized by the integral scale, the wrinkling scale defined as the inverse of the peak surface density is shown to scale with Reynolds number across simulations and as turbulence decays. As a result, the area ratio and the burning rate are found to increase as ${Re}_{\lambda }^{1.13}$ , in agreement with recent experiments on spherical turbulent premixed flames. We observe that the area ratio does not vary with turbulent intensity when holding the Reynolds number constant.  more » « less
Award ID(s):
1805921
NSF-PAR ID:
10206898
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
906
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The burning rate in a spherically turbulent premixed flame is explored using direct numerical simulations, and a model of ordinary differential equations is proposed. The numerical dataset, from a previous work, is obtained from direct numerical simulations of confined spherical flames in isotropic turbulence over a range of Reynolds numbers. We begin the derivation of the model with an equation for the burning rate for the domain under consideration, and using a thin flame assumption and a two-fluid approach, we find the normalized turbulent burning rate to be controlled by the increase in flame surface area due to turbulent wrinkling, and correction factor which is observed to be consistently less than unity. A Reynolds scaling hypothesis for the flame turbulent wrinkling from a previous work using the same numerical dataset is used to model the term controlling the increase in flame surface area. The correction factor is hypothesized to reflect flame stretch effects, and hence this factor is modeled using Markstein theory applied to global averaged quantities. The ordinary differential equations are rewritten to reflect easily observable quantities such as the chamber pressure and mean flame radius, and then expressed in dimensionless form to assess dependence on various dimensionless parameters. The model predictions are found to be in good agreement with the numerical data within expected variances. Additionally, Markstein theory is found to be sufficient in describing the effects of flame stretch in the turbulent premixed flames under consideration. 
    more » « less
  2. The objective of this work is to provide physical insight into the mechanisms governing flame–turbulence interactions and explore the impact of the ubiquitous Darrieus–Landau instability on the propagation. It is based on the hydrodynamic theory of premixed flames that considers the flame thickness much smaller than all other fluidynamical length scales. In this asymptotic limit, the flame is thus confined to a surface whilst the diffusion and reaction processes occurring inside the flame zone are accounted for by two parameters: the unburned-to-burned density ratio and the Markstein length. The robust model, which is free of phenomenology and turbulence modelling assumptions, makes transparent the mutual interactions between the flame and the fluid flow, and permits examining trends in flame and flow characteristics while varying the turbulence intensity and mixture properties. It is used in this study to examine the morphological changes of the flame surface that result from the intertwined effects of the turbulence and instability, as demonstrated by the local displacement and curvature of the flame front, the extent of wrinkling and folding of the flame surface, and the overall flame brush thickness. It also provides a direct evaluation of the turbulent flame speed and its dependence on the mean flame curvature and on the hydrodynamic strain that it experiences. Also discussed are the effects of the flame on the flow by examining the various mechanisms of enstrophy and scalar gradient production/destruction, the degree of anisotropy created in the burned gas, and the restructuring of the vortical motion beyond the flame. 
    more » « less
  3. Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( ${>}300\delta _i$ , with $\delta _i$ the inflow boundary-layer thickness) at three different Mach numbers, $2.5$ , $4.9$ and $10.9$ , with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number $Re_\tau \approx 1200$ ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region. 
    more » « less
  4. Soot formation is quantified in detail (volume fraction, particle size, number concentration, and light emissivity dispersion exponent) in a series of partially premixed counterflow flames of ethylene at equivalence ratios equal to 6.5, 5.0, and 4.0, and with maximum temperature spanning approximately 200 K. The focus is to investigate the effect of peak temperature and equivalence ratio on soot formation while maintaining constant global strain and stoichiometric mixture fraction. Oxygen is progressively displaced from the oxidizer to the fuel stream of a diffusion flame to stabilize partially premixed flames of decreasing, showing a double-flame structure consisting of a rich premixed flame component stabilized on the fuel side of the stagnation plane and a diffusion flame component stabilized on the oxidizer side. Soot is detected in the region sandwiched between the two flame components, is formed in both of them, and is convected away radially at the Particle Stagnation Plane (PSP). At fixed , raising the peak temperature invariably raises the soot volume fraction throughout the probed region. Vice versa, at fixed peak temperature, lowering the equivalence ratio causes the premixed flame component to shift away from the diffusion flame component, with the consequent broadening of the soot forming region and an increase in both soot volume fraction as well as soot particle sizes through an enhancement of surface growth. Detailed probing of the region in the vicinity of the PSP offers evidence of soot oxidation from molecular oxygen. Furthermore, when the maximum temperature is sufficiently low, the net soot production rate turns negative because surface oxidation overwhelms surface growth. Comparing the soot number production rate inferred from experiments to the dimerization rate of benzene, naphthalene, and pyrene reveals that only the smallest aromatics are present in flames at sufficiently large concentrations to account for soot nucleation. This observation applies to both the diffusion flame and the premixed flame components and confirms previous findings in strictly diffusion flames. 
    more » « less
  5. We investigate the modes of deformation of an initially spherical bubble immersed in a homogeneous and isotropic turbulent background flow. We perform direct numerical simulations of the two-phase incompressible Navier–Stokes equations, considering a low-density bubble in the high-density turbulent flow at various Weber numbers (the ratio of turbulent and surface tension forces) using the air–water density ratio. We discuss a theoretical framework for the bubble deformation in a turbulent flow using a spherical harmonic decomposition. We propose, for each mode of bubble deformation, a forcing term given by the statistics of velocity and pressure fluctuations, evaluated on a sphere of the same radius. This approach formally relates the bubble deformation and the background turbulent velocity fluctuations, in the limit of small deformations. The growth of the total surface deformation and of each individual mode is computed from the direct numerical simulations using an appropriate Voronoi decomposition of the bubble surface. We show that two successive temporal regimes occur: the first regime corresponds to deformations driven only by inertial forces, with the interface deformation growing linearly in time, in agreement with the model predictions, whereas the second regime results from a balance between inertial forces and surface tension. The transition time between the two regimes is given by the period of the first Rayleigh mode of bubble oscillation. We discuss how our approach can be used to relate the bubble lifetime to the turbulence statistics and eventually show that at high Weber numbers, bubble lifetime can be deduced from the statistics of turbulent fluctuations at the bubble scale. 
    more » « less