Linear-Eddy Model Formulated Probability Density Function and Scalar Dissipation Rate Models for Premixed Combustion

A tabulated, pseudo-turbulent Probability Density Function (PDF) model for premixed combustion is proposed. The Linear-Eddy Model (LEM) is used to construct the PDFs for a temperature-based progress variable in a premixed, turbulent methane/air V-flame produced by the Cambridge slot burner. As a second case study, the LEM PDFs are similarly compared to PDFs extracted from Direct Numerical Simulations (DNS) of a turbulent premixed flame. LEM demonstrates the ability to reproduce the salient features from experimental and DNS PDFs; moreover, it is able to better capture turbulent effects than previously suggested laminar flamelet PDF models. The Scalar Dissipation Rate (SDR) for premixed combustion is likewise investigated. The stochastic nature of LEM enables it to mimic the overall behaviors of turbulent reactions inexpensively and qualitatively. Crucially, LEM appears to be well suited for the preprocessing tabulation of PDF and SDR models for a number of premixed combustion simulation strategies.     

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature,Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO₂. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.     

Joint PDF Closure of Turbulent Premixed Flames

In this paper, a novel model for turbulent premixed combustion in the corrugated flamelet regime is presented, which is based on transporting a joint probability density function (PDF) of velocity, turbulence frequency and a scalar vector. Due to the high dimensionality of the corresponding sample space, the PDF equation is solved with a Monte-Carlo method, where individual fluid elements are represented by computational particles. Unlike in most other PDF methods, the source term not only describes reaction rates, but accounts for “ignition” of reactive unburnt fluid elements due to propagating embedded quasi laminar flames within a turbulent flame brush. Unperturbed embedded flame structures and a constant laminar flame speed (as expected in the corrugated flamelet regime) are assumed. The probability for an individual particle to “ignite” during a time step is calculated based on an estimate of the mean flame surface density (FSD), latter gets transported by the PDF method. Whereas this model concept has recently been published [21], here, a new model to account for local production and dissipation of the FSD is proposed. The following particle properties are introduced: a flag indicating whether a particle represents the unburnt mixture; a flame residence time, which allows to resolve the embedded quasi laminar flame structure; and a flag indicating whether the flame residence time lies within a specified range. Latter is used to transport the FSD, but to account for flame stretching, curvature effects, collapse and cusp formation, a mixing model for the residence time is employed. The same mixing model also accounts for molecular mixing of the products with a co-flow. To validate the proposed PDF model, simulation results of three piloted methane-air Bunsen flames are compared with experimental data and very good agreement is observed.     

A Comprehensive Study of Effects of Mixing and Chemical Kinetic Models on Predictions of n-heptane Jet Ignitions with the PDF Method

The composition Probability Density Function (PDF) model is coupled with a Reynolds-averaged k???ε turbulence model and three computationally efficient, yet widely used chemical mechanisms to simulate transient n-heptane spray injection and ignition in a high temperature and high density ambient fluid. Molecular diffusion is modelled by three mixing models, namely the interaction by exchange with the mean (IEM), modified Curl (MC) and Euclidean minimum spanning trees (EMST) models. The liquid phase is modelled by a discrete phase model (DPM). This represents among the first applications of the PDF method in practical diesel engine conditions. A non-reacting case is first considered, with the focus on the ability of the model to capture the spray structure, e.g., vapour penetration and liquid length, fuel mixture fraction and its variance. Reacting cases are then investigated to compare and evaluate the three different chemical mechanisms and the three mixing models. It is concluded that the EMST mixing model in conjunction with a reduced chemical kinetic model (Lu et al., Combust Flame 156(8):1542–1551, 2009) performs the best among the options considered. The sensitivity of the results to the choice of the mixing constant is also studied to understand its effect on the flame ignition and stabilisation. Finally, the PDF model is compared to a well-mixed model that assumes turbulent fluctuations are negligible, which has been widely used in the diesel spray combustion community. Significant structural differences in the modelled flame are revealed comparing the PDF method with the well-mixed model. Quantitatively, the PDF model exhibits excellent agreement with the measurements and shows much better results than the well-mixed model in all ambient O₂ and temperature conditions.     

Nonpremixed bluff-body burner flow and flame imaging study

A nonpremixed bluff-body burner flow and flame have been studied using planar flow visualization and species concentration imaging techniques. The burner consists of a central jet of CH ₄ in a cylindrical bluff-body and an outer coflowing-air stream. Planar flow visualization, using Mie scattering from seed particles added to the fuel jet, Raman scattering from CH ₄ and laser-induced fluorescence of CH combined with Raman scattering of CH ₄ provided information on turbulent flow, mixing and combustion. The CH ₄ imaging system utilized two cameras, which enhanced the dynamic range of the diagnostic system by a factor of 10 over a single-camera system. It was observed that the fuel jet stagnated on the axis due to interaction with the high velocity air flow. The flow and mixing were found to have significant coherent and noncoherent, large-scale, time-varying structures. The detailed CH ₄ Raman and CH fluorescence measurements of an air-dominated bluff-body flame revealed that the stagnation zone governs mixing and flame stability. Through large-scale mixing, the stagnated jet feeds the recirculation zone and also creates a favorable condition to stabilize the flame detached from the bluff-body. The instantaneous flame zone, as defined by CH, was found to be narrow and concentrated in an envelope around the stagnation zone. This narrow flame characteristic is consistent with that observed for jet flames. Although the internal structure of the flame envelops have not yet been defined, these results suggest that this bluff-body flame can be modeled by a flame sheet type approach, where the reaction front is captured by the large-scale structures. This should simplify the development of modeling approaches for these flows since molecular mixing and chemical reaction, which occur within the flame sheet, can be separated from the large-scale mixing process.     

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.     

Transported PDF Modeling of Ethanol Spray in Hot-Diluted Coflow Flame

This paper presents a numerical modeling study of one ethanol spray flame from the Delft Spray-in-Hot-Coflow (DSHC) database, which has been used to study Moderate or Intense Low-oxygen Dilution (MILD) combustion of liquid fuels (Correia Rodrigues et al. Combust. Flame 162(3), 759–773, 2015). A “Lagrangian-Lagrangian” approach is adopted where both the joint velocity-scalar Probability Density Function (PDF) for the continuous phase and the joint PDF of droplet properties are modeled and solved. The evolution of the gas phase composition is described by a Flamelet Generated Manifold (FGM) and the interaction by exchange with the mean (IEM) micro-mixing model. Effects of finite conductivity on droplet heating and evaporation are accounted for. The inlet boundary conditions starting in the dilute spray region are obtained from the available experimental data together with the results of a calculation of the spray including the dense region using ANSYS Fluent 15. A method is developed to determine a good estimation for the initial droplet temperature. The inclusion of the “1/3” rule for droplet evaporation and dispersion models is shown to be very important. The current modeling approach is capable of accurately predicting main properties, including mean velocity, droplet mean diameter and number density. The gas temperature is under-predicted in the region where the enthalpy loss due to droplet evaporation is important. The flame structure analysis reveals the existence of two heat release regions, respectively having the characteristics of a premixed and a diffusion flame. The experimental and modeled temperature PDFs are compared, highlighting the capabilities and limitations of the proposed model.     

Development and evaluation of a new turbulence generator for atomization research

Planar Mie scattering visualizations in compressible mixing layers are used to compute the probability density function of a passive scalar. Mixing layer flows with relative Mach numbers of 0.63 and 1.49 are studied. Ethanol condensation is used to generate both scalar transport seeding and product formation seeding. All PDFs exhibit a marching behavior. The condensation process in the product formation seeding is modeled to provide an estimate of the error embedded in the scalar transport PDFs. The mixing efficiency is found to be 0.56 in the product formation experiments, and the overprediction of mixing efficiency by the scalar PDFs is estimated to be 11% based on results from the ethanol condensation model.List of Symbols _291-01 Damköhler number based on - J droplet nucleation rate - k Boltzmann constant - m _c molecular mass of ethanol - M _r relative Mach number, M _r = 2U/(a₁ + a₂) - N ^* number of nucleated droplets - p(,) probability density function - P _d internal droplet pressure - P _m total mixed fluid probability - P _sat ethanol saturation partial pressure - P _v ethanol vapor partial pressure - r freestream velocity ratio, r=U ₂/U₁; droplet radius - r ^* critical nucleation radius - R gas constant for air - _291-2 Reynolds number based on - s freestream density ratio, s = ₂/₁ - T local static temperature - U ₁ high speed freestream velocity - U ₂ low speed freestream velocity - U _c large structure convection velocity, - U freestream velocity difference, U=U ₁–U₂ - x streamwise coordinate - y transverse coordinate - mixing layer thickness - _i incompressible mixing layer thickness - mixture fraction - similarity variable, = (y–y ₀)/ - _c condensed phase ethanol density - droplet surface tension     

Stationary modes of combustion in a fine-scale turbulent stream

The majority of models of the turbulent combustion of gases are based mainly on intuitive concepts concerning the processes occurring in the flame. The characteristics of a turbulent flame are estimated from considerations of dimensionality and similarity. A detailed review of works on turbulent combustion is given in [1]. Problems on the calculation of the combustion rate in a turbulent stream as a proper value of the equations of heat and mass transfer and of the corresponding boundary conditions have recently been raised. Here too one must rest on assumptions of a semiempirical nature, which in large measure is connected with the inadequate level of development of turbulence theory. In the present work the equation of propagation of the zone of chemical reactions in the stream is averaged statistically by analogy with studies of turbulent flows. Correct averaging is possible at scales of hydrodynamic disturbances smaller than the flame thickness (fine-scale turbulence). The temperature pulsations are related with the size of the heat flux using the theory of mixing lengths. The main influence is specific to effects arising during averaging of the heat release function. Two stationary modes, distinguished by the normal propagation velocity ₁, are isolated within the framework of the Cauchy problem with a given initial mixture temperature and zero heat flux in the burned gas. A heat conduction mode occurs with a stream velocity > ₁ and an induction mode with < ₁. An expression is found for ₁ which reflects the principal effects in the flame and which in the limit coincides with the equation of Zel'dovich and Frank-Kamenetskii for a laminar flame. In those cases when the distorting effect of the heat release function is small, the turbulence affects the combustion rate through mechanisms of intensification of transport processes.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 118–124, September–October, 1973.     

Effect of velocity ratio on plane mixing layer development: Influence of the splitter plate wake

An experimental study has been conducted to investigate the effect of velocity ratio on the approach of a plane mixing layer to self-similarity. Plane mixing layers with five different velocity ratios (0.5, 0.6, 0.7, 0.8 and 0.9) were generated in a newly designed mixing layer wind tunnel with both initial boundary layers tripped. For each velocity ratio, mean flow and turbulence measurements were obtained at eight streamwise locations with a single cross-wire probe. The results indicate that the splitter plate wake plays a very dominant and, in some cases, a lasting role in the development of the mixing layer. For velocity ratios between 0.5 and 0.7, self-similarity of the mixing layer was observed with the asymptotic states comparable. Mixing layers with the higher velocity ratios failed to achieve a self-similar state within the measurement domain, although a slow approach to it was apparent. The development distance decreased with increasing velocity ratio up to 0.7, after which it appeared to increase. Almost all of the observed effects may be attributed to the presence of the splitter plate wake and its complex interaction with the mixing layer.List of symbols C _f boundary layer skin friction coefficient - H boundary layer shape factor - r velocity ratio of the two streams, (=U ₂/U ₁) - Re _L Reynolds number, (=UL/v) - R correlation coefficient in least squares fit - U, V, W mean velocity in the X, Y, Z directions, respectively - U ^* velocity parameter, [=(U–U ₂)/(U ₁–U ₂)] - U ₀ velocity difference, (=U ₁–U ₂) - U _e free-stream velocity in the wind tunnel - u, , w fluctuating velocity components in the X, Y, Z directions, respectively - u, , w instantaneous velocity in the X, Y, Z directions, respectively, e.g. u=U+u - X ₀ virtual origin of the mixing layer - X, Y, Z cartesian coordinates for streamwise, normal, and spanwise directions, respectively - Y ₀ centerline of mixing layer from error function fit - mixing layer width from error function fit - ₉₉ initial boundary layer thickness - similarity coordinate [=(Y–Y ₀)/] - initial boundary layer momentum thickness - modified velocity ratio [=(1–r)/(1+r)] - _n initial instability wavelength in the mixing layer - spreading parameter [=1/(d/dX)] - ₀ spreading parameter for single-stream mixing layer - - (overbar) Time-averaged quantity - ( )_max maximum value at given X-station - ( )_min minimum value at given X-station - ( )₁ value for high-speed side - ( )₂ value for low-speed side     

Computational Modelling of Turbulent Mixing in Confined Swirling Environment Under Constant and Variable Density Conditions

A high-intensity swirling flow in a model combustor subjected to large density variations has been examined computationally. The focus is on the Favre-averaged Navier–Stokes computations of the momentum and scalar transport employing turbulence models based on the differential second-moment closure (SMC) strategy. An updated version of the basic, high-Reynolds number SMC model accounting for a quadratic expansion of both the pressure–strain and dissipation tensors and a near-wall SMC model were used for predicting the mean velocity and turbulence fields. The accompanied mixing between the annular swirling air flow and the central non-swirling helium jet was studied by applying three scalar flux models differing mainly in the model formulation for the pressure-scalar gradient correlation. The computed axial and circumferential velocities agree fairly well with the reference experiment [So et al., NASA Contractor Report 3832, 1984; Ahmed and So, Exp. Fluids 4 (1986) 107], reproducing important features of such a weakly supercritical flow configuration (tendency of the flow core to separate). Although the length at which the mixing was completed was reproduced in reasonable agreement with the experimental results, the mixing activity in terms of the spreading rate of the shear/mixing layer, that is its thickness, was somewhat more intensive. Prior to these investigations, the model applied was validated by computing the transport of the passive scalar in the non-swirling (Johnson and Bennet, Report NASA CR-165574, UTRC Report R81-915540-9, 1981) and swirling (Roback and Johnson, NASA Contractor Report 168252, 1983) flow in a model combustor.     

Velocity and spanwise vorticity measurements in an excited mixing layer of a plane jet

The mixing layer of a plane jet was subjected to periodic weak excitation at two different frequencies corresponding to shear layer mode (St =0.012) and preferred mode (St _D =0.36). The nozzle exit boundary layer was identical for the unexcited and excited flows. Measurements of mean velocity, longitudinal and lateral velocity fluctuations, Reynolds shear stress and spanwise component of fluctuating vorticity were made over a longitudinal distance x/D of 6 for both the unexcited and the excited flows. Even weak excitation was observed to influence the development of the mixing layer. Under shear layer mode of excitation, the width of the layer and longitudinal turbulence level decrease compared to the naturally developing (unexcited) flow whereas preferred mode of excitation results in increase in the width and turbulence levels. The rms spanwise vorticity showed an increase for shear layer mode of excitation whereas the preferred mode of excitation resulted in a decrease compared to the values in an unexcited flow. Spectra of velocity and vorticity fluctuations exhibited subharmonic peaks, suggesting the possible occurrence of vortex pairing in both unexpected and excited flows. The influence of excitation is found to decrease as x/D increases and is not significant at x/D=6.This study was partly supported by a grant from the Research Grant Council, Hong Kong. The support and hospitality of the Department of Mechanical Engineering, University of Hong Kong are gratefully acknowledged by SR. The authors are grateful the referees for valuable comments.     

Density effects on jet characteristics in confined swirling flow

Density effects on isothermal jet mixing in confined swirling flow are investigated. The experiment is carried out with helium/air as the jet fluid in the same facility as that used by So et al. (1984a, b) and the test conditions are chosen to be the same as before. Contrary to the homogeneous mixing results, the helium jet is preserved up to 40 jet diameters downstream. The behavior of the mean and turbulence field depends highly on the initial jet velocity. Since the jets are fully turbulent and the jet momentum fluxes for inhomogeneous mixing are less than those for homogeneous mixing, the cause of this difference in behavior is directly attributed to the combined action of density difference and swirl. In spite of this, near isotropy of the turbulence field is again observed at 40 jet diameters downstream.     

Investigations of mean and turbulent flow characteristics of a two dimensional plane diffuser

This paper reports the investigation of mean and turbulent flow characteristics of a two-dimensional plane diffuser. Both experimental and theoretical details are considered. The experimental investigation consists of the measurement of mean velocity profiles, wall static pressure and turbulence stresses. Theoretical study involves the prediction of downstream velocity profiles and the distribution of turbulence kinetic energy using a well tested finite difference procedure. Two models, viz., Prandtl's mixing length hypothesis and k- model of turbulence, have been used and compared. The nondimensional static pressure distribution, the longitudinal pressure gradient, the pressure recovery coefficient, percentage recovery of static pressure, the variation of U _max/U _bar along the length of the diffuser and the blockage factor have been valuated from the predicted results and compared with the experimental data. Further, the predicted and the measured value of kinetic energy of turbulence have also been compared. It is seen that for the prediction of mean flow characteristics and to evaluate the performance of the diffuser, a simple turbulence model like Prandtl's mixing length hypothesis is quite adequate.List of symbols C ₁ , C ₂ ,C turbulence model constants - F _x body force - k kinetic energy of turbulence - l _m mixing length - L length of the diffuser - u, v, w rms value of the fluctuating velocity - u, v, w turbulent component of the velocity - mean velocity in the x direction - _A average velocity at inlet - U _bar average velocity in any cross section - U _max maximum velocity in any cross section - V mean velocity in the y direction - W local width of the diffuser at any cross section - x, y coordinates - dissipation rate of turbulence - _m eddy diffusivity - Von Karman constant - mixing length constant - _l laminar viscosity - _eff effective viscosity - v kinematic viscosity - density - _k effective Schmidt number for k - effective Schmidt number for - stream function - non dimensional stream function     

A Simple Approach for Specifying Velocity Inflow Boundary Conditions in Simulations of Turbulent Opposed-Jet Flows

A new methodology is developed to specify inflow boundary conditions for the velocity field at the nozzle exit planes in turbulent counterflow simulations. The turbulent counterflow configuration consists of two coaxial opposed nozzles which emit highly-turbulent streams of varying species compositions depending on the mode considered. The specification of velocity inflow boundary conditions at the nozzle exits in the counterflow configuration is non-trivial because of the unique turbulence field generated by the turbulence generating plates (TGPs) upstream of the nozzle exits. In the method presented here, a single large-eddy simulation (LES) is performed in a large domain that spans the region between the TGPs of the nozzles, and the time series of the velocity fields at the nozzle exit planes are recorded. To provide inflow boundary conditions at the nozzle exit planes for simulations under other conditions (e.g., different stream compositions, bulk velocity, TGP location), transformations are performed on the recorded time series: the mean and r.m.s. (root-mean-square) quantities of velocity, as well as the longitudinal integral length scale on the centerline, at the nozzle exits in simulations are matched to those observed in experiments, thereby matching the turbulent Reynolds number R e _t . The method is assessed by implementing it in coupled large-eddy simulation/probability density function (LES/PDF) simulations on a small cylindrical domain between the nozzle exit planes for three different modes of the counterflow configuration: N ₂ vs. N ₂; N ₂ vs. hot combustion products; and C H ₄/N ₂ vs. O ₂. The inflow method is found to be successful as the first and second moments of velocity from the LES/PDF simulations agree well with the experimental data on the centerline for all three modes. This simple yet effective inflow strategy can be applied to eliminate the computational cost required to simulate the flow field upstream of the nozzle exits. It is also emphasized that, in addition to the predicted time series data, the availability of experimental data close to the nozzle exit planes plays a key role in the success of this method.     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

High-resolution imaging of dissipative structures in a turbulent jet flame with laser Rayleigh scattering

High-resolution 2-D imaging of laser Rayleigh scattering is used to measure the detailed structure of the thermal dissipation field in a turbulent non-premixed CH₄/H₂/N₂ jet flame. Measurements are performed in the near field (x/d = 5–20) of the flame where the primary combustion reactions interact with the turbulent flow. The contributions of both the axial and radial gradients to the mean thermal dissipation are determined from the 2-D dissipation measurements. The relative contributions of the two components vary significantly with radial position. The dissipation field exhibits thin layers of high dissipation. Noise suppression by adaptive smoothing enables accurate determination of the dissipation-layer widths from single-shot measurements. Probability density functions (PDF) of the dissipation-layer widths conditioned on temperature are approximately log-normal distributions. The conditional layer width PDFs are self-similar functions with the layer widths scaling with temperature to the 0.75 power. The high signal-to-noise ratio of the Rayleigh scattering images coupled with an interlacing technique for noise suppression enable fully resolved measurements of the mean power spectral density (PSD) of the temperature gradients. These spectra are used to determine the turbulence microscales by measuring a cutoff wavelength, λ _C , at 2% of the peak PSD. The Batchelor scale is estimated from λ _C , and the results are compared with estimates from scaling laws in non-reacting flows. At x/d = 20, the different approaches to determining the Batchelor scale are comparable on the jet centerline. However, the estimates from non-reacting flow scaling laws are significantly less accurate in off-centerline regions and at locations closer to the nozzle exit. Throughout the near field of the jet flame, the measured ratio of a characteristic dissipation-layer width to the local Batchelor scale is larger than values previously reported for the far field of non-reacting flows.     

A Turbulence Model Sensitivity Study for CH4/H2 Bluff-Body Stabilized Flames

The objective of this work is to assess the performances of different turbulence models in predicting turbulent diffusion flames in conjunction with the flamelet model.The k– model, the Explicit Algebraic Stress Model (EASM) and the k– model withvaried anisotropy parameter C (LEA k– model)are first applied to the inert turbulent flow over a backward-facing step, demonstrating the quality of the turbulence models. Following this, theyare used to simulate the CH₄/H₂ bluff-body flame studied by the University of Sydney/Sandia.The numerical results are compared to experimental values of the mixture fraction, velocity field, temperature and constituent mass fractions.The comparisons show that the overall result depends on the turbulence model used, and indicate that theEASM and the LEA k– models perform better than the k– model and mimic most of the significant flow features.