Spatial and Temporal Characteristics of OH in Turbulent Opposed-Jet Double Flames

Simultaneous high repetition-rate, two-point hydroxyl (OH) time-series measurements with associated PLIF/PIV measurements are employed to investigate spatio–temporal scales and flame-velocity interactions in turbulent opposed jets sustaining methane-air double flames. For a fuel-side equivalence ratio, ϕ _B  = 1.2, a rich premixed flame exists on the fuel side while a diffusion flame exists on the air side of the stagnation plane. The bulk Reynolds number (Re) and strain rate (SR) can be adjusted to generate flames at ϕ _B  = 1.2 with both well separated and completely merged flame fronts. Simultaneous PLIF/PIV measurements highlight distinct spatial OH structures of the premixed and diffusive fronts corresponding to variations in the flow field. The self-propagating tendency of the rich premixed front causes large-scale wrinkling, thereby enhancing the OH contour length by 15% as compared to the diffusive front. Two-point OH time-series measurements are implemented to quantify both spatial and temporal fluctuations via study of radial length and time scales. In general, these integral length and time scales follow similar trends and reach a minimum at the axial location of peak [OH]. In comparison to merged double flames having higher Re and SR, greater OH fluctuations are observed in the rich-premixed front as compared to the diffusive front for a well separated double flame. Because of the developing turbulence, the OH length scales exhibit reduced axial gradients across the reaction zone for higher Re in comparison to lower Re. A stochastic time-series simulation, using a state relationship based on a joint mixture fraction and progress variable, is utilized to extract estimated scalar time scales from those of measured OH. The simulations indicate that the hydroxyl fluctuations in double flames are only twice those of the underlying conserved scalar. “Turbulent Opposed-Jet Double Flames” is submitted for consideration as a full length article to Flow Turbulence and Combustion.     

A Non-Adiabatic Flamelet Progress�CVariable Approach for LES of Turbulent Premixed Flames

A progress variable/flame surface density/probability density function method has been employed for a Large Eddy Simulation of a CH₄/Air turbulent premixed bluff body flame. In particular, both mean and variance of the progress variable are transported and subgrid spatially filtered gradient contributes to model the flame surface density (that introduces the effect of the subgrid flame reaction zone) and to presume a probability density function (that introduces the effect of subgrid fluctuations on chemistry). Chemistry is preliminarly tabulated in terms of laminar premixed flames and enthalpy is included as a new coordinate in their tabulation to take into account heat losses in the flowfield. Then, the PDF is used to build a turbulent flamelet library. The filtered mass, momentum, enthalpy and scalar equations mentioned above are integrated by an explicit scheme using finite differences, 2nd?Corder accurate in space and third order in time, over a cylindrical non-uniform grid using a staggered mesh. The bluff-body geometry is modelled by using the Immersed Boundary Method. The numerical predictions are compared with the available experimental data.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.     

Combined CARS/LDA instrument for simultaneous temperature and velocity measurements

The performance of a combined CARS/LDA instrument capable of measuring temperature and two velocity components with a time coincidence of about 4 s is evaluated in a turbulent premixed propane-air Bunsen-burner flame. Measurements near the base of the flame exhibit negative axial correlations, indicative of normal gradient transport; those near the flame tip show strong positive axial correlations, indicative of transport counter to the temperature gradient. The radial correlations are positive both in the reaction zone and in the plume. An analysis of temperature data from measurements made (1) independent of and (2) coincidental with LDA measurements indicates that the CARS/LDA instrument provides a density-weighted velocity, temperature, and velocity temperature correlation due to the density variations in the flame.     

Time-resolved measurements of hydroxyl in stable and extinguishing partially premixed turbulent opposed-jet flames

Quantitative hydroxyl time-series measurements from a set of stable and extinguishing turbulent opposed-flow partially premixed CH₄/air flames are used to investigate the effect of Reynolds number and fuel-side equivalence ratio on the structure of turbulent partially premixed flames. The hydroxyl (OH) integral time scale, computed from the autocorrelation function, is used to characterize OH fluctuations and is found to reach a minimum at the axial location of peak OH. Analyses of the duration of and period between bursts in the OH time series are used to examine the dynamics of flame-front motion. In general, with increasing Reynolds number (Re), the distribution in OH burst times shifts towards smaller time scales. A hydroxyl intermittency parameter is also defined from the bursts to quantify the presence or absence of OH. For flames with the same fuel-side equivalence ratio, the hydroxyl intermittency at peak OH remains almost constant when going from stable to extinguishing flames. However, histograms portray an increase in burst separation times for flames displaying occasional extinction events. Hydroxyl time series for a partially premixed flame at a fuel-side equivalence ratio of 2.0 and Re = 6650 are synthesized by using mixture-fraction simulations based on calculated state relationships for OH versus mixture fraction (f). The laminar-flamelet model is employed to explore relations between OH and f so as to predict trends in mixture-fraction time scales.“Time-Series Measurements in Turbulent Opposed-Jet Flames" is submitted for consideration as a full length article to Flow Turbulence and Combustion.     

CH double-pulsed PLIF measurement in turbulent premixed flame

The flame displacement speeds in turbulent premixed flames have been measured directly by the CH double-pulsed planar laser-induced fluorescence (PLIF). The CH double-pulsed PLIF systems consist of two independent conventional CH PLIF measurement systems and laser beams from each laser system are led to same optical pass using the difference of polarization. The highly time-resolved measurements are conducted in relatively high Reynolds number turbulent premixed flames on a swirl-stabilized combustor. Since the time interval of the successive CH PLIF can be selected to any optimum value for the purpose intended, both of the large scale dynamics and local displacement of the flame front can be discussed. By selecting an appropriate time interval (100–200 μs), deformations of the flame front are captured clearly. Successive CH fluorescence images reveal the burning/generating process of the unburned mixtures or the handgrip structures in burnt gas, which have been predicted by three-dimensional direct numerical simulations of turbulent premixed flames. To evaluate the local flame displacement speed directly from the successive CH images, a flame front identification scheme and a displacement vector evaluation scheme are developed. Direct measurements of flame displacement speed are conducted by selecting a minute time interval (≈30 μs) for different Reynolds number (Re _λ = 63.1–115.0). Local flame displacement speeds coincide well for different Reynolds number cases. Furthermore, comparisons of the mean flame displacement speed and the mean fluid velocity show that the convection in the turbulent flames will affect the flame displacement speed for high Reynolds number flames.     

An Experimental Investigation of the Turbulence Structure of a Lifted H2/N2 Jet Flame in a Vitiated Co-Flow

Measurements of mean velocity components, turbulent intensity, and Reynolds shear stress are presented in a turbulent lifted H₂/N₂ jet flame as well as non-reacting air jet issuing into a vitiated co-flow by laser doppler velocimetry (LDV) technique. The objectives of this paper are to obtain a velocity data base missing in the previous experiment data of the Dibble burner and so provide initial and flow field data for evaluating the validity of various numerical codes describing the turbulent partially premixed flames on this burner. It is found that the potential core is shortened due to the high ratio of jet density to co-flow density in the non-reacting cases. However, the existence of flame suppressed turbulence in the upstream region of the jet dominates the length of potential core in the reacting cases. At the centreline, the normalized axial velocities in the reacting cases are higher than the non-reacting cases, and the relative turbulent intensities of the reacting flow are smaller than in the non-reacting flow, where a self-preserving behaviour for the relative turbulent intensities exists at the downstream region. The profiles of mean axial velocity in the lifted flame distribute between the non-reacting jet and non-premixed flame both in the axial and radial distributions. The radial distributions of turbulent kinetic energy in the lifted flames exhibit a change in distributions indicating the difference of stabilisation mechanisms of the two lifted flame. The experimental results presented will guide the development of an improved modelling for such flames.     

LES of Low to High Turbulent Combustion in an Elevated Pressure Environment

A subgrid scale flame surface density combustion model for the Large Eddy Simulation (LES) of premixed combustion is derived and validated. The model is based on fractal characteristics of the flame surface, assuming a self similar wrinkling of the flame between smallest and largest wrinkling length scales. Experimental and direct numerical simulation databases as well as theoretical models are used to derive a model for the fractal parameters, namely the cut-off lengths and the fractal dimension suitable in the LES context. The combustion model is designed with the intent to simulate low as well as high Reynolds number premixed turbulent flame propagation and with a focus on correct scaling with pressure. The combustion model is validated by simulations of turbulent Bunsen flames with methane and propane fuel at pressure levels between 0.1 MPa and 2 MPa and at turbulence levels of $0 < u^{\prime }/s_{L}^{0} < 11$ , conditions typical for spark ignition engines. The predicted turbulent flame speed is in a very good agreement with the experimental data and a smooth transition from resolved flame wrinkling to fully modelled, nearly subgrid-only wrinkling is realized. Evaluating the influence of mesh resolution shows a predicted mean flame surface and turbulent flame speed independent of mesh resolution for cases with 9–86 % resolved flame surface. Additional simulations of a highly turbulent jet flame at 0.1 MPa and 0.5 MPa and the comparison with experimental data in terms of flame shape, velocity field and turbulent fluctuations validates the model also at conditions typical for gas turbines.     

Conditional Moment Closure for Large Eddy Simulations

A conditional moment closure (CMC) based combustion model for large-eddy simulations (LES) of turbulent reacting flow is proposed and evaluated. Transport equations for the conditionally filtered species are derived that are consistent with the LES formulation and closures are suggested for the modelling of the conditional velocity, conditional scalar dissipation and the fluctuations around the conditional mean. A conventional β-probability density distribution of the scalar is used together with dynamic modelling for the sub-grid fluxes. The model is validated by comparison of simulations with measurements of a piloted, turbulent methane-air jet diffusion flame.     

Space-time velocity correlations in the impeller stream of a Rushton turbine

 Simultaneous velocity measurements at two points have been carried out with a two-component laser Doppler velocimeter in order to characterise the turbulence along the impeller stream of a Rushton turbine in a water tank. In addition to mean values of the velocity and its fluctuations, space–time correlations have been studied to enable the determination of the convection velocity as well as of the scales of the turbulent structures and to investigate the validity of Taylor’s hypothesis. A correction to the relationship between space and time correlations is proposed. Received: 22 March 1996/Accepted: 23 April 1997     

Topology of turbulent premixed flame fronts resolved by simultaneous planar imaging of LIPF of OH radical and rayleigh scattering

 Planar images of Rayleigh scattering and laser-induced predissociative OH-fluorescence (OH-LIPF) have been obtained simultaneously in turbulent premixed jet flames on a single-shot basis. The geometric structure of temperature and OH isocontours were extracted for fractal analysis. A power-law fractal behavior can be identified in the ensemble-averaged flame length measure. It was found that the inner and outer cut-off scales of OH contours are larger than those of the iso-temperature contours; while the OH images show comparatively smaller fractal dimensions. The joint-pdf ’s between flame temperature and OH LIPF signals at different heights are also derived to evaluate the flame stretch effect on local flame structure. Comparison of image pairs near the extinction limit suggests that Rayleigh thermometry is more adequate to characterize the fine-scale flame front wrinkling in highly stretched turbulent premixed flames. Received: 12 September 1997/Accepted: 19 May 1998     

Transition from Countergradient to Gradient Scalar Transport in Developing Premixed Turbulent Flames

A simple model of turbulent scalar flux developed recently by the present authors is applied to determine the direction of the flux in a statistically planar one-dimensional premixed flame that does not affect turbulence and has self-similar mean structure. Results obtained in the case of statistically stationary turbulence indicate that transition from countergradient to gradient turbulent scalar transport may occur during flame development, as the peak mean rate of product creation moves to the trailing edge of the flame brush. In the case of decaying turbulence, the opposite transition (from gradient to countergradient transport) was simulated in line with available DNS data. In both cases, transition instant depends strongly on turbulence and mixture characteristics. In particular, countergradient transport is suppressed by an increase in the rms turbulent velocity and by a decrease in the laminar flame speed or density ratio, in line with available experimental and DNS data. The obtained results lend qualitative support to the model of turbulent scalar flux addressed in the present work.     

Scale similarity based models and their application to subgrid scale scalar flux modelling in the context of turbulent premixed flames

The performance of a variety of scale similarity (SS) type models for closure of sub-grid scalar flux in the context of Large Eddy Simulations (LES) of premixed turbulent combustion has been assessed. In addition to the well-known SS models, a more recent development by Anderson and Domaradzki (2012) is included in the analysis and also further model extensions and improvements are discussed. The work is based on a priori analysis of two Direct Numerical Simulation (DNS) databases of freely propagating turbulent premixed flames with a range of different Lewis and turbulent Reynolds numbers. Depending on the balance between the effects of flame normal acceleration due to heat release and the effects of turbulent velocity fluctuations, as well as the filter size, the subgrid-scalar flux exhibits both local gradient and counter-gradient transport which presents a considerable modelling challenge. The assessment is based on a correlation analysis and on the magnitude of the model expressions conditional on the Favre averaged reaction progress variable in comparison to the value obtained from DNS. Despite the fact that most of the models have been developed in the context of momentum transport in non-reactive flows they show either comparable or better performance in comparison to more conventional models used for reactive scalar flux closure. It is found that some models are sensitive to the test filter width and recommendations are provided in this regard. Further it is observed that the use of a Favre test filter substantially increases the correlation strength in direction of mean flame propagation where effects of heat release are most pronounced.     

Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation

The effects of varying turbulence intensity and turbulence length scale on premixed turbulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation.

     

Measurements of turbulent premixed flame dynamics using cinema stereoscopic PIV

A new experimental method is described that provides high-speed movies of turbulent premixed flame wrinkling dynamics and the associated vorticity fields. This method employs cinema stereoscopic particle image velocimetry and has been applied to a turbulent slot Bunsen flame. Three-component velocity fields were measured with high temporal and spatial resolutions of 0.9 ms and 140 μm, respectively. The flame-front location was determined using a new multi-step method based on particle image gradients, which is described. Comparisons are made between flame fronts found with this method and simultaneous CH-PLIF images. These show that the flame contour determined corresponds well to the true location of maximum gas density gradient. Time histories of typical eddy–flame interactions are reported and several important phenomena identified. Outwardly rotating eddy pairs wrinkle the flame and are attenuated at they pass through the flamelet. Significant flame-generated vorticity is produced downstream of the wrinkled tip. Similar wrinkles are caused by larger groups of outwardly rotating eddies. Inwardly rotating pairs cause significant convex wrinkles that grow as the flame propagates. These wrinkles encounter other eddies that alter their behavior. The effects of the hydrodynamic and diffusive instabilities are observed and found to be significant contributors to the formation and propagation of wrinkles.     

Measurements and scaling of wall shear stress fluctuations

Measurements of velocity and wall shear stress fluctuations were made in an external turbulent boundary layer developed over a towed surface-piercing flat plate. An array of eight flush-mounted wall shear stress sensors was used to compute the space-time correlation function. A methodology for in situ calibration of the sensors for ship hydrodynamic applications is presented. The intensity of the wall shear stress fluctuations, τ _rms/τ _avg was measured as 0.25 and 0.36 for R _θ =3,150 and 2,160 respectively. The probability density is shown to exhibit positive skewness, and lack of flow reversals at the wall. Correlations between velocity and wall shear stress fluctuations are shown to collapse with outer boundary layer length and velocity scales, verifying the existence of large-scale coherent structures which convect and decay along the wall at an angle of inclination varying from 10 to 13° over the range of Reynolds numbers investigated. The wall shear stress convection velocity determined from narrow band correlation measurements is shown to scale with outer variables. The space-time correlation of the wall shear is shown to exhibit a well-defined convective ridge, and to decay 80% over approximately 3δ for R _θ =3,150. Published online: 7 November 2002     

A Posteriori Assessment of Algebraic Scalar Dissipation Models for RANS Simulation of Premixed Turbulent Combustion

This paper examines the effects of scalar dissipation rate modelling on mean reaction rate predictions in turbulent premixed flames. The sensitivity of the mean reaction rate is explored by using different closures for scalar dissipation and the sensitivity of the scalar dissipation models themselves is also examined with respect to their defining constants. The influence of different scalar dissipation models on the flame location and mean velocities is reported and compared with experimental results. The predicted reaction rate is found to be sensitive to the choice of closure used for scalar dissipation and also to the respective constants used in the scalar dissipation models. It is also found that the scalar dissipation models involving chemical and turbulent time scales yield a more physically plausible reaction rate when compared with the scalar dissipation models relying only on the turbulent time scale.     

Digital image analysis of a turbulent flame

Digital image analysis of cine pictures of an unconfined rich premixed turbulent flame has been used to determine structural characteristics of the turbulent/non-turbulent interface of the flame. The results, comprising various moments of the interface position, probability density functions and correlation functions, establish that the instantaneous flame-interface position is essentially a Gaussian random variable with a superimposed quasi-periodical component. The latter is ascribable to a pulsation caused by the convection and the stretching of ring vortices present within the flame. To a first approximation, the flame can be considered similar to a three-dimensional axisymmetric turbulent jet, with superimposed ring vortices, in which combustion occurs.     

Experimental characterization of non-premixed turbulent jet propane flames

This paper reports an experimental study conducted on turbulent jet propane flames aiming at further understanding of turbulent structure in non-premixed slow-chemistry combustion systems. Measurements of mean and fluctuating velocity and temperature fields, mean concentration of major chemical species, correlation between velocity and temperature fluctuations, and dissipation of temperature fluctuations are reported in a turbulent round jet non-premixed propane flame, Re=20 400 and 37 600, issuing vertically in still air. The experimental conditions were designed to provide a complete definition of the upstream boundary conditions in the measurement domain for the purpose of validating computational models. The measured data depicts useful flow field information for describing turbulent non-premixed slow-chemistry flames. Velocity–temperature correlation measurements show turbulent heat fluxes tended to be restricted to the mixing layer where large temperature gradients occurred. Observations of non-gradient diffusion of heat at x/D=10 were verified. Temperature fluctuation dissipation, χ, showed the highest values in the shear layer, where the variance of temperature fluctuations was maximum and combustion occurred. The isotropy between the temperature dissipation in the radial and tangential directions was confirmed. By contrast, the observed anisotropy between axial and radial directions of dissipation suggests the influence of large structures in the entrainment shear layer on the production of temperature fluctuations in the flame region. The value of the normalized scalar dissipation at the stoichiometric mixture fraction surface, χ_st, was calculated, and ranges between 2 and 4 s⁻¹. The measured data were used to estimate the budgets in the balance equations for turbulent kinetic energy, Reynolds shear stresses, turbulent heat flux and temperature variance, quantifying the mechanisms involved in the generation of turbulence as well as in the transport of the temperature.