Local curvature measurements of a lean,partially premixed swirl-stabilised flame

A swirl-stabilised, lean, partially premixed combustor operating at atmospheric conditions has been used to investigate the local curvature distributions in lifted, stable and thermoacoustically oscillating CH₄-air partially premixed flames for bulk cold-flow Reynolds numbers of 15,000 and 23,000. Single-shot OH planar laser-induced fluorescence has been used to capture instantaneous images of these three different flame types. Use of binary thresholding to identify the reactant and product regions in the OH planar laser-induced fluorescence images, in order to extract accurate flame-front locations, is shown to be unsatisfactory for the examined flames. The Canny-Deriche edge detection filter has also been examined and is seen to still leave an unacceptable quantity of artificial flame-fronts. A novel approach has been developed for image analysis where a combination of a non-linear diffusion filter, Sobel gradient and threshold-based curve elimination routines have been used to extract traces of the flame-front to obtain local curvature distributions. A visual comparison of the effectiveness of flame-front identification is made between the novel approach, the threshold binarisation filter and the Canny-Deriche filter. The novel approach appears to most accurately identify the flame-fronts. Example histograms of the curvature for six flame conditions and of the total image area are presented and are found to have a broader range of local flame curvatures for increasing bulk Reynolds numbers. Significantly positive values of mean curvature and marginally positive values of skewness of the histogram have been measured for one lifted flame case, but this is generally accounted for by the effect of flame brush curvature. The mean local flame-front curvature reduces with increasing axial distance from the burner exit plane for all flame types. These changes are more pronounced in the lifted flames but are marginal for the thermoacoustically oscillating flames. It is concluded that additional fuel mixture fraction and velocimetry studies are required to examine whether processes such as the degree of partial-premixedness close to the burner exit plane, the velocity field and the turbulence field have a strong correlation with the curvature characteristics of the investigated flames.     

Estimation of three-dimensional flame surface densities from planar images in turbulent premixed combustion

Turbulence motions are, by nature, three-dimensional while planar imaging techniques, widely used in turbulent combustion, give only access to two-dimensional information. For example, to extract flame surface densities, a key ingredient of some turbulent combustion models, from planar images implicitly assumes an instantaneously two-dimensional flow, neglecting the unresolved flame front wrinkling. The objective here is to estimate flame surface densities from two-dimensional measurements assuming that (1) the flow is statistically two dimensional; (2) the measuring plane is a plane of symmetry of the mean flow, either by translation (homogeneous third direction as in slot burners for example) or by rotation (axi-symmetrical flows such as jets) and (3) flame movements in transverse directions are similar. The unknown flame front wrinkling is then modelled from known quantities. An excellent agreement is achieved against direct numerical simulation (DNS) data where all three-dimensional quantities are known, but validations in other conditions (larger DNS, experiments) are required.     

Investigation of bed shear stresses subject to external turbulence

The fluctuating bed shear stress has largely been investigated only for uniform channel flows and boundary layers. In practical engineering, the flow conditions are often modified due to the presence of various hydraulic structures. To what extent the modification affects the known characteristics of the bed shear stress is not clear. This study presents experimental results of the bed shear stress fluctuations, which are obviously subjected to external turbulence generated by superimposing artificial structures in the open channel flows. The statistical analysis of the measurements shows that the probability density function of the bed shear stress can be described by the lognormal function. Some associated relations concerning higher-order moments, skewness and kurtosis, which are derived from the lognormal function, are further compared with the experimental data. Physical implication of the skewed probability density distribution is finally discussed.     

Flame Curvature Distribution in High Pressure Turbulent Bunsen Premixed Flames

The flame curvature statistics of turbulent premixed Bunsen flames have been analysed in this paper using a Direct Numerical Simulation (DNS) database of turbulent Bunsen flames at ambient and elevated pressures. In order to be able to perform a large parametric study in terms of pressure, heat release parameter, turbulence conditions and nozzle diameter, a single step Arrhenius type irreversible chemistry has been used for the purpose of computational economy, where thermo-chemical parameters are adjusted to match the behavior of stoichiometric methane-air flames. This analysis focuses on the characterization of the local flame geometry in response to turbulence and hydro-dynamic instability. The shape of the flame front is found to be consistent with existing experimental data. Although the Darrieus Landau instability promotes cusp formation, a qualitatively similar flame morphology can be observed for hydro-dynamically stable flames. A criterion has been suggested for the curvature PDF to become negatively skewed.     

The Scalar Gradient Alignment Statistics of Flame Kernels and its Modelling Implications for Turbulent Premixed Combustion

The effects of mean flame curvature on reaction progress variable gradient, $\nabla c$ , alignment with local turbulent strain rate are studied based on three-dimensional Direct Numerical Simulation (DNS) data of turbulent premixed flame kernels with different initial radii under decaying turbulence. A statistically planar flame is also considered in order to compare the results obtained from the kernels with a flame of zero mean curvature. It is found that the dilatation rate effects diminish with decreasing kernel radius due to defocusing of heat in the positively curved regions. This gives rise to a decrease in the extent of reaction progress variable gradient alignment with most extensive principal strain rate with decreasing kernel radius. The modelling implications of the statistics of the alignment of $\nabla c$ with local strain rate have been studied in terms of scalar dissipation rate transport. A new modelling methodology for the contribution of the scalar-turbulence interaction term in the transport equation for the mean scalar dissipation is suggested addressing the reduced effects of dilatation rate for flame kernels and the diminished value of turbulent straining at the small length scales at which turbulence interacts with small flame kernels. The performance of the new models is found to be satisfactory while comparing to DNS results. The existing models for the dilatation contribution and the combined chemical reaction and molecular dissipation contributions to the transport of mean scalar dissipation, which were originally proposed for statistically planar flames, are found to satisfactorily predict the corresponding quantities for turbulent flame kernels.     

Low Reynolds number axisymmetric turbulent boundary layer on a cylinder

In the present study, an axisymmetric turbulent boundary layer growing on a cylinder is investigated experimentally using hot wire anemometry. The combined effects of transverse curvature as well as low Reynolds number on the mean and turbulent flow quantities are studied. The measurements include the mean velocity, turbulence intensity, skewness and flatness factors in addition to wall shear stress. The results are presented separately for the near wall region and the outer region using dimensionless parameters suitable for each case. They are also compared with the results available in the open literature.The present investigation revealed that the mean velocity in near wall region is similar to other simple turbulent flows (flat plate boundary layer, pipe and channel flows); but it differs in the logarithmic and outer regions. Further, for dimensionless moments of higher orders, such as skewness and flatness factors, the main effects of the low Reynolds number and the transverse curvature are present in the near wall region as well as the outer region.     

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.

     

Extension of Stoney's formula to non-uniform temperature distributions in thin film/substrate systems. The case of radial symmetry

Current methodologies used for the inference of thin film stress through curvature measurements are strictly restricted to stress and curvature states which are assumed to remain uniform over the entire film/substrate system. By considering a circular thin film/substrate system subject to non-uniform, but axisymmetric temperature distributions, we derive relations between the film stresses and temperature, and between the plate system's curvatures and the temperature. These relations featured a “local” part which involves a direct dependence of the stress or curvature components on the temperature at the same point, and a “non-local” part which reflects the effect of temperature of other points on the location of scrutiny. Most notably, we also derive relations between the polar components of the film stress and those of system curvatures which allow for the experimental inference of such stresses from full-field curvature measurements in the presence of arbitrary radial non-uniformities. These relations also feature a “non-local” dependence on curvatures making full-field measurements of curvature a necessity for the correct inference of stress. Finally, it is shown that the interfacial shear tractions between the film and the substrate are proportional to the radial gradients of the first curvature invariant and can also be inferred experimentally.     

Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

Investigation of the Syngas Flame Characteristics at Elevated Pressures Using Optical and Laser Diagnostic Methods

The effect of pressure on the characteristics of syngas flames is investigated under gas turbine relevant conditions using planar laser induced fluorescence of OH radicals and OH^* chemiluminescence imaging. An optically accessible combustor fitted with a swirl burner was operated with two different syngas mixtures, preheated air at 700?K, and pressures ranging from 5 to 20?bars. The thermal load varied from 15 to 25?kW/bar at an equivalence ratios 0.5. The OH-PLIF measurements show that the flames under all conditions exhibited two reaction fronts, one at the shear layer between the inner recirculation zone and the fuel inlet, and one between the fuel inlet and the air nozzle. The more or less continuous reaction front at low pressure turned into a highly corrugated flame front at higher pressures, with isolated regions of ignition and extinction. The probability density distribution of the flame curvature for the mixtures studied showed that the inner and outer flame responded differently to the pressure increase, with the mean curvature magnitude also depending on the mixture composition and thermal load. The measurements clearly shows the limitations associated with the use of OH^* chemiluminescence images as a marker for the heat release rate especially in case of syngas mixtures.     

Non-uniform, axisymmetric misfit strain： in thin films bonded on plate substrates/substrate systems： the relation between non-uniform film stresses and system curvatures

Current methodologies used for the inference of thin film stress through curvature measurements are strictly restricted to stress and curvature states which are assumed to remain uniform over the entire film/substrate system. By considering a circular thin film/substrate system subject to non-uniform, but axisymmetric misfit strain distributions in the thin film, we derived relations between the film stresses and the misfit strain, and between the plate system’s curvatures and the misfit strain. These relations feature a ‘‘local’’ part which involves a direct dependence of the stress or curvature components on the misfit strain at the same point, and a ‘‘non-local’’ part which reflects the effect of misfit strain of other points on the location of scrutiny. Most notably, we also derived relations between the polar components of the film stress and those of system curvatures which allow for the experimental inference of such stresses from full-field curvature measurements in the presence of arbitrary radial non-uniformities. These relations also feature a ‘‘non-local’’ dependence on curvatures making a full-field measurement a necessity. Finally, it is shown that the interfacial shear tractions between the film and the substrate are proportional to the radial gradients of the first curvature invariant and can also be inferred experimentally.     

Three-dimensional temporally resolved measurements of turbulence–flame interactions using orthogonal-plane cinema-stereoscopic PIV

A new orthogonal-plane cinema-stereoscopic particle image velocimetry (OPCS-PIV) diagnostic has been used to measure the dynamics of three-dimensional turbulence–flame interactions. The diagnostic employed two orthogonal PIV planes, with one aligned perpendicular and one aligned parallel to the streamwise flow direction. In the plane normal to the flow, temporally resolved slices of the nine-component velocity gradient tensor were determined using Taylor’s hypothesis. Volumetric reconstruction of the 3D turbulence was performed using these slices. The PIV plane parallel to the streamwise flow direction was then used to measure the evolution of the turbulence; the path and strength of 3D turbulent structures as they interacted with the flame were determined from their image in this second plane. Structures of both vorticity and strain-rate magnitude were extracted from the flow. The geometry of these structures agreed well with predictions from direct numerical simulations. The interaction of turbulent structures with the flame also was observed. In three dimensions, these interactions had complex geometries that could not be reflected in either planar measurements or simple flame–vortex configurations.     

On spatial resolution issues related to time-averaged quantities using hot-wire anemometry

The effect of spatial resolution on streamwise velocity measurements with single hot-wires is targeted in the present study, where efforts have been made to distinguish between spatial resolution and Reynolds number effects. The basis for measurements to accurately determine the mean velocity and higher order moments is that the probability density distribution is measured correctly. It is well known that the turbulence intensity is increasingly attenuated with increasing wire length. Here, it is also shown (probably for the first time) that besides the probability density distribution and hence the higher order moments, even the mean velocity is affected, albeit to subtle extent, but with important consequences in studies of concurrent wall-bounded turbulence.     

Influences of acoustic instabilities on turbulent-flame propagation

Influences of acoustic instabilities on premixed turbulent-flames have been studied experimentally in a Taylor–Couette (TC) combustor for downward flame propagation in a turbulent flow-field generated in the annulus between two cylinders. Flow-field velocities were measured at a fixed location upstream of the propagating flame using laser-doppler velocimetry, while flame speeds were determined from video-recorded images. It is found that the existence of pre-ignition turbulence in the combustor (generated by rotation of the combustor-cylinder walls) does not eliminate acoustic instabilities, however as the level of pre-ignition turbulence is increased the influence of the secondary acoustic instability on the turbulent-flame speed becomes insignificant. For low intensities of pre-ignition turbulence the flame is found to accelerate during the latter stages of flame propagation, while for high levels of pre-ignition turbulence the flame propagates at a statistically constant speed, even though velocity fluctuations have been substantially amplified by the time the flame reaches the bottom end of the combustor as a result of acoustic instability.     

The lateral movement of the three-dimensional cellular flame at low Lewis numbers

The lateral movement of the three-dimensional (3-D) cellular flame at low Lewis numbers is numerically investigated. The equation used is the compressible Navier–Stokes equation including a one-step irreversible chemical reaction. We superimpose the hexagonal disturbance with the peculiar wave number on the stationary plane flame and calculate the evolution of the disturbed flame. When the Lewis number is unity, i.e., only the hydrodynamic effect has an influence on the flame instability, the stationary cellular flame is formed. When the Lewis number is lower than unity, i.e., the diffusive-thermal and hydrodynamic effects have an influence, the laterally moving cellular flame is formed. With a decrease in the Lewis number, the laterally moving velocity of the cell increases. The laterally moving velocity of the three-dimensional cellular flame is much larger than that of the two-dimensional (2-D) cellular flame. Because, the increment of local temperature at the convex flame front toward the unburned gas in the three-dimensional flame is great compared with that in the two-dimensional flame.     

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.     

障碍物对预混火焰特性影响的大涡数值模拟

障碍物在预混气体火焰传播过程中对其结构及传播特性造成较大影响,对火焰的加速和爆燃转爆轰过程(deflagration-to-detonation transition, DDT)起到直接的促进作用。通过障碍物条件下可视管道中甲烷/空气预混火焰传播实验,捕获其火焰微观结构变化。采用三维物理模型,采用壁面自适应局部涡黏模型(wall-adapting local eddy-viscosity, WALE)的大涡模拟(large eddy simulation, LES),并用火焰增厚化学反应模型(thickened flame model, TFM)对实验过程进行重现。分析开口管道中预混火焰翻越障碍物后的复杂流场变化,并分析层流向湍流转变过程的特点。揭示了在障碍物影响下预混火焰扰动失稳现象的直接原因,是由障碍物引发的3个气流涡团同时作用而形成Kelvin-Helmholtz不稳定及Rayleigh-Taylor不稳定现象耦合作用所导致。     

Modification of near-wall turbulence structure in a shear-driven three-dimensional turbulent boundary layer

 Most high Reynolds number flows of engineering interest are three-dimensional in nature. Key features of three-dimensional turbulent boundary layers (3DTBLs) include: non-colateral shear stress and strain rate vectors, and decreasing ratio of the shear stresses to the turbulent kinetic energy with increasing three-dimensionality. These are indicators that the skewing has a significant effect on the structure of turbulence. In order to further investigate the flow physics and turbulence structure of these complex flows, an innovative method for generating a planar shear-driven 3DTBL was developed. A specialized facility incorporating a relatively simple geometry and allowing for varying strengths of crossflow was constructed to facilitate studies where the skewing is decoupled from the confounding effects of streamwise pressure gradient and curvature. On-line planar particle image velocimetry (PIV) measurements and flow visualization results indicate that the experimental configuration generates the desired complex flow, which exhibits typical characteristics associated with 3DTBLs. Furthermore, spanwise shear results in modification of the near-wall turbulence structure. Analysis of near-wall flow visualization photographs revealed a reduction of mean streak length with increasing spanwise shear, while streak spacing remained relatively constant. In the most strongly sheared case, where the belt velocity is twice that of the freestream velocity, the mean streak length was reduced by approximately 50%. Received: 28 October 1997/Accepted: 4 February 1998     

Pulsed hot-wire measurements in two- and three-dimensional wall jets

Pulsed Hot-Wire Anemometry (PHWA) measurements are performed in well defined two- and three-dimensional turbulent wall jets. For the two-dimensional wall jet the objective is to study reported differences between conventional Hot-Wire Anemometry (HWA) and Laser Doppler Anemometry (LDA) results. In the three dimensional wall jet, new improved data are provided, employing a measuring technique suitable for highly turbulent flows. This, since only hot-wire results previously have been published for this flow. The pulsed wire results show good agreement with existing Laser Doppler anemometer data in the two-dimensional wall jet, both reporting significantly higher turbulence levels in the outer region of the flow than hot-wires do. The hot-wire anemometer errors generally increase with increasing local turbulence intensity and since the three-dimensional wall jet has a higher turbulence level than its two-dimensional equivalent, the new pulsed hot-wire results improve the information available for the turbulence field in this flow significantly. Received: 29 January 1998/ Accepted: 19 February 1999