首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 406 毫秒
1.
We investigated the local flame speed of a two-dimensional, methane-air triple flame in a rectangular burner. The velocity fields and the concentration profiles were measured with particle image velocimetry and the Rayleigh scattering method, respectively. There was a requisite combination of initial velocity and initial concentration gradient for consistency of the local concentration gradient at the leading edge of the flame. In these cases, the flame curvatures were also consistent. Accordingly, the burning velocity, defined as local flow velocity at the triple point, was determined by the flame curvature. The burning velocity increased with increasing flame curvature, when the curvature was near zero. After that, the burning velocity decreased with increasing curvature. The peak value thus exceeded the adiabatic one-dimensional laminar burning velocity. Comparing the effects of the measured flame stretch rate on the flow strain κs and flame curvature κc, κs is larger and increases more rapidly than κc for flame curvatures satisfying 1/Rf < 250 m−1 and then becomes constant while κc still increases for 250 m−1 < 1/Rf, so that κc becomes much larger than κs. There is also a peak in burning velocity at roughly the transition in flame curvature specified above. Therefore, the burning velocity for a low concentration gradient correlates with the flame stretch rate.  相似文献   

2.
Premixed turbulent flames of methane–air and propane–air stabilized on a bunsen type burner were studied using planar Rayleigh scattering and particle image velocimetry. The fuel–air equivalence ratio range was from lean 0.6 to stoichiometric for methane flames, and from 0.7 to stoichiometric for propane flames. The non-dimensional turbulence rms velocity, u′/SL, covered a range from 3 to 24, corresponding to conditions of corrugated flamelets and thin reaction zones regimes. Flame front thickness increased slightly with increasing non-dimensional turbulence rms velocity in both methane and propane flames, although the flame thickening was more prominent in propane flames. The probability density function of curvature showed a Gaussian-like distribution at all turbulence intensities in both methane and propane flames, at all sections of the flame.The value of the term , the product of molecular diffusivity evaluated at reaction zone conditions and the flame front curvature, has been shown to be smaller than the magnitude of the laminar burning velocity. This finding questions the validity of extending the level set formulation, developed for corrugated flames region, into the thin reaction zone regime by increasing the local flame propagation by adding the term to laminar burning velocity.  相似文献   

3.
Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, Pb, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When Pb is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate Pb and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.  相似文献   

4.
DNS is performed to analyse the effects of Lewis number (Le), density ratio and gravity in stagnating turbulent premixed flames. The results show good agreement with those of Lee and Huh (Combustion and Flame, Vol. 159, 2012, pp. 1576–1591) with respect to the turbulent burning velocity, ST, in terms of turbulent diffusivity, flamelet thickness, mean curvature and displacement speed at the leading edge. In all four stagnating flames studied, a mean tangential strain rate resulting in a mean flamelet thickness smaller than the unstretched laminar flame thickness leads to an increase in ST. A flame cusp of positive curvature involves a superadiabatic burned gas temperature due to diffusive–thermal instability for an Le less than unity. Wrinkling tends to be suppressed at a larger density ratio, not enhanced by hydrodynamic instability, in the stagnating flow configuration. Turbulence is produced, resulting in highly anisotropic turbulence with heavier unburned gas accelerating through a flame brush by Rayleigh–Taylor instability. Results are also provided on brush thickness, flame surface density and conditional velocities in burned and unburned gas and on flame surfaces to represent the internal brush structures for all four test flames.  相似文献   

5.
Turbulent burning velocities for methane/air mixtures at pressures ranging from atmospheric pressure up to 1.0 MPa and mixture temperatures of 300 and 573 K were measured, which covers the typical operating conditions of premixed-type gas-turbine combustors. A bunsen-type flame stabilized in a high-pressure chamber was used, and OH-PLIF visualization was performed with the pressure and mixture temperature being kept constant. In addition to a burner with an outlet diameter of 20 mm for the high-pressure experiments, a large-scale burner with an outlet diameter of 60 mm was used at atmospheric pressure to extend the turbulence Reynolds number based on the Taylor microscale, Rλ, as a common parameter to compare the pressure and temperature effects. It was confirmed that Rλ over 100 could be attained and that u′/SL could be extended even at atmospheric pressure. Based on the contours of the mean progress variable c = 0.1 determined using OH-PLIF images, turbulent burning velocity was measured. ST/SL was also found to be greatly affected by pressure for preheated mixtures at 573 K. The bending tendency of the ST/SL curves with u′/SL was seen regardless of pressure and mixture temperature and the Rλ region where the bending occurs corresponded well to the region where the smallest scale of flame wrinkling measured as a fractal inner-cutoff approaches the characteristic flame instability scale and becomes almost constant. A power law of ST/SL with (P/P0)(u′/SL) was clearly seen when ST was determined using c = 0.1 contours, and the exponent was close to 0.4, indicating agreement with the previous results using the mean flame cone method and the significant pressure effects on turbulent burning velocity.  相似文献   

6.
N2 Q-branch CARS spectra have been recorded and evaluated for temperature determination in a turbulent, premixed CH4/air stagnation flame with a burner of 40 mm diameter and 22 kW thermal load. Temperature histograms on the flame axis at different distances from the stagnation plate have been measured. Problems of practical applicability are addressed, including those arising from the limited spatial resolution of the BOXCARS geometry, from an insufficient dynamic range of the diode array detector, and from a memory effect of the detector in the case of measurements in highly turbulent flame areas with strong intermittency. Some information is given on the computerized acquisition and on the evaluation of the large amounts of data that are necessary for extensive investigations in large combustion systems.  相似文献   

7.
The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH2O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH2O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH2O structure is somewhat thinner than the experiments.  相似文献   

8.
The thickness of the instantaneous flamelets in a turbulent flame brush on a weak-swirl burner burning in the thin reaction zones regime has been analysed experimentally, theoretically, and numerically. The experimental flame thickness has been measured correlating two simultaneous Rayleigh images and one OH-image from two closely spaced cross sections in the flame. It appears that the low temperature edge of the flame is thickened by turbulent eddies but that these structures cannot penetrate far enough into the flame front to distort the inner layer for the moderate Karlovitz numbers used. The flame front based on the temperature gradient at the inner layer becomes thinner for lean flames and thicker for rich methane–air flames. This has been explained theoretically and numerically by studying the influence of flame stretch and preferential diffusion on the flame thickness. It appears that the flame front thickness at the inner layer (and mass burning rate) is not influenced by turbulent mixing processes, and it seems that eddies of the size of the inner layer have to be used to change this picture. Experiments closer to the boundary of the broken reaction zones regime have to confirm this in the future.  相似文献   

9.

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localized determination of subgrid scale closure parameters is essential. We formulate a semi-localized model based on the dynamical equation for the subgrid scale turbulence energy k sgs. The turbulent flame speed s t is of the order √2k sgs. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, s t∝√2k sgs, are investigated.  相似文献   

10.
The performance of a dynamic subgrid model for the turbulent burning speed of a premixed flame is investigated for a series of idealized test cases where the flame front is wrinkled by a multiple-scale shear flow; a rigorous asymptotic subgrid model is also implemented for comparison. Explicit formulae for the flame wrinkled shape and turbulent speed are available to generate a reference database. The role of the subgrid wrinkling models is to achieve the same overall flame shape and propagation speed in a simulation where only the largest scales of the flow are explicitly accounted for. Very good results are obtained when the subgrid burning speed enhancement is estimated using the asymptotic subgrid model. On the other hand, the dynamic model attempts to exploit the scaling observable in the simulation to extrapolate the turbulent burning speed enhancement in the original system. The performance of this strategy is adequate for some regimes but poor for others; the source of the problem is traced back to the existence of a scaling transition that occurs as the flame propagating speed is adjusted during the large-eddy simulation. A modification to the scaling of the enhanced burning is implemented to account for the existence of the two distinct scaling ranges; it improves significantly the predictions of the dynamic model away from the transition, but results in the near-critical range remain predictably very poor compared with the rigorous asymptotic model results. These conclusions based on a priori performance for the reference steady data are confirmed by comparing unsteady large-eddy and direct simulations. Results based on rigorous mathematical tools are possible here because of the separation of length scales in the special class of idealized flow fields used in this study: their relevance to more realistic flows is also discussed.  相似文献   

11.
A unique burner was constructed to experimentally realize a one-dimensional unstrained planar non-premixed flame, previously considered only in idealized theoretical models. One reactant, the fuel mixture in the current experiments, is supplied through a porous plug at the bottom of the combustion chamber and flows vertically up towards the horizontal flame. The crux of the design is the introduction of the oxidizer from above in such a way that its diffusion against the upward product flow is essentially one-dimensional, i.e., uniform over the burner cross-section. This feature was implemented by introducing the oxidizer into the burner chamber from the top through an array of 625 closely spaced hypodermic needles, and allowing the hot products to escape vertically up through the space between the needles. Due to the injection of oxidizer through discrete tubes, a three-dimensional “injection layer” exists below the exit plane of the oxidizer supply tubes. Experimental evidence suggests that this layer is thin and that oxidizer is supplied to the flame by 1-D counterdiffusion, producing a nearly unstrained flame. To characterize the burner, flame position measurements were conducted for different compositions and flowrates of H2–CO2 and O2–CO2 mixtures. The measured flame locations are compared to an idealized one-dimensional model in which only diffusion of oxidizer against the product flow is considered. The potential of the new burner is demonstrated by a study of cellular structures forming near the extinction limit. Consistent with previous investigations, cellular instabilities are shown to become more prevalent as the initial mixture strength and/or the Damköhler number are decreased. As the extinction limit is approached, the number of cells was observed to decrease progressively.  相似文献   

12.
An experimental study for 1-butanol single droplet flames in constant and oscillatory flow fields was conducted under microgravity conditions at elevated pressure. In the constant flow experiments, flow velocities from 0 to 40 cm/s were tested. Using obtained data of d2, the burning rate constants were evaluated. The burning rate constant in the quiescent condition was also calculated successfully at high pressure by the extrapolation method based on the Frössling relation. In the oscillatory flow experiments, the flow velocities were varied from 0 to 40 cm/s at the frequencies of 2–40 Hz. Results showed that the burning rate constant during the droplet lifetime varied following the quasi-steady relation at 0.1 MPa; however, in the conditions with higher frequencies at 0.4 MPa, the average burning velocity became larger than that for the constant flow case with the velocity equivalent to the maximum velocity in the oscillatory flow. Under the condition where the burning rate constant increased, it was observed that the flame did not sufficiently move back upstream, leading to enhancement of the heat transfer from the flame to the droplet surface. Therefore, the instantaneous burning rate constant increased. To investigate the mechanism of such flame behavior, the ratio of two characteristic times, τf/τD (τf: flow oscillation characteristic time, τD: diffusion characteristic time), were compared. As the flow oscillatory frequency increased, τf/τD becomes smaller. τf/τD also became smaller at high pressure. If τf/τD is small due to the small mass diffusion rate, the droplet flame could not move back to the appropriate position for the minimum velocity in steady flow, leading to an increase of the burning rate constant, especially in the case of higher frequency at high pressure.  相似文献   

13.
Turbulent vortex rings were investigated in weakly turbulent flow and in three different grid generated turbulent flows to clarify the reciprocal action of the vortex ring with defined external turbulence. Assuming self-similarity and turbulent viscosity as proportional to V0D0 the equations for the ring diameter D(t) and the velocity of propagation V(t) were derived. The time difference Δt between the virtual origins of 1/V(t) respectively D2(t) led to an invariant term. The equation of momentum is fulfilled. – Position and diameter of the vortex rings were determined till their decay by means of an optical system, which did not disturb the vortex rings. The experimental results in weakly turbulent ambient flow obtained by the author and by others confirm the theory very well. The ambient turbulence was nearly constant in the measuring region; its effect could be described by simply adding its viscosity to the vortex ring’s internal turbulent viscosity. The results could be represented in unified non-dimensional diagrams. Moreover, an explanation was found as to why the mean internal turbulent viscosity is constant.  相似文献   

14.
The burning and extinction characteristics of isolated small nonane droplets are examined in a buoyant convective environment and in an environment with no external axial convection (as created by doing experiments at low gravity) to promote spherical droplet flames. The ambience is air and a mixture of 30%O2/70%He to assess the influence of soot formation. The initial droplet diameter (Do) ranges from 0.4 to 0.95 mm. Measurements are reported of the extinction diameter and time to extinction, and of the evolution of droplet diameter, flame diameter, soot shell diameter, burning rate, and broadband radiative emissions.In a buoyancy-free environment for air larger droplets burn slower than smaller droplets for the range of Do examined, which is attributed to the influence of soot. In the presence of a buoyant flow in air, no influence of Do is observed on the burning rate while the buoyant flames are still heavily sooting. The effect of Do is believed to be due to a combination of dominance of the nonluminous, nonsooting, portion of the buoyant flame around the forward half of the droplet on heat transport and the secondary role of the luminous wake portion of the flame. In a non-sooting helium inert at low gravity, no effect of Do is found on the evolution of droplet diameter.Flame extinction is observed only in the 30%O2/70%He ambience. For all of the observations, extinction appears to occur before the disappearance of the droplet which is then followed by a period of evaporation. The extinction diameter and time to extinction increases with Do and an empirical correlation is presented for these two variables.  相似文献   

15.
The present study experimentally investigates the structure and instabilities associated with extremely low-stretch (1 s−1) gaseous diffusion flames. Ultra-low-stretch flames are established in normal gravity by bottom burning of a methane/nitrogen mixture discharged from a porous spherically symmetric burner of large radius of curvature. OH-PLIF and IR imaging techniques are used to characterize the reaction zone and the burner surface temperature, respectively. A flame stability diagram mapping the response of the ultra-low-stretch diffusion flame to varying fuel injection rate and nitrogen dilution is explored. In this diagram, two main boundaries are identified. These boundaries separate the stability diagram into three regions: sooting flame, non-sooting flame, and extinction. Two distinct extinction mechanisms are noted. For low fuel injection rates, flame extinction is caused by heat loss to the burner surface. For relatively high injection rates, at which the heat loss to burner surface is negligible, flame radiative heat loss is the dominant extinction mechanism. There also exists a critical inert dilution level beyond which the flame cannot be sustained. The existence of multi-dimensional flame phenomena near the extinction limits is also identified. Various multi-dimensional flame patterns are observed, and their evolutions are studied using direct chemiluminescence and OH-PLIF imaging. The results demonstrate the usefulness of the present burner configuration for the study of low-stretch gaseous diffusion flames.  相似文献   

16.
This paper describes the simultaneous application of time-sequenced laser-induced fluorescence imaging of OH radicals and stereoscopic particle image velocimetry for measurements of the flame front dynamics in lean and premixed LP turbulent flames. The studied flames could be acoustically driven, to simulate phenomena important in LP combustion technologies. In combination with novel image post processing techniques we show how the data obtained can be used to track the flame front contour in a plane defined by the illuminating laser sheets. We consider effects of chemistry and convective fluid motion on the dynamics of the observed displacements and analyse the influence of turbulence and acoustic forcing on the observed contour velocity, a quantity we term as s d 2D. We show that this quantity is a valuable and sensitive indicator of flame turbulence interactions, as (a) it is measurable with existing experimental methodologies, and (b) because computational data, e.g. from large eddy simulations, can be post processed in an identical fashion. s d 2D is related (to a two-dimensional projection) of the three-dimensional flame displacement speed s d , but artifacts due to out of plane convective motion of the flame surface and the uncertainty in the angle of the flame surface normal have to be carefully considered. Monte Carlo simulations were performed to estimate such effects for several distributions of flame front angle distributions, and it is shown conclusively that s d 2D is a sensitive indicator of a quantity related to s d in the flames we study. s d 2D was shown to increase linearly both with turbulent intensity and with the amplitude of acousting forcing for the range of conditions studied.  相似文献   

17.
This study has been mainly motivated to assess computationally and theoretically the conditional moment closure (CMC) model and the transient flamelet model for the simulation of turbulent nonpremixed flames. These two turbulent combustion models are implemented into the unstructured grid finite volume method that efficiently handles physically and geometrically complex turbulent reacting flows. Moreover, the parallel algorithm has been implemented to improve computational efficiency as well as to reduce the memory load of the CMC procedure. Example cases include two turbulent CO/H2/N2 jet flames having different flow timescales and the turbulent nonpremixed H2/CO flame stabilized on an axisymmetric bluff-body burner. The Lagrangian flamelet model and the simplified CMC formulation are applied to the strongly parabolic jet flame calculation. On the other hand, the Eulerian particle flamelet model and full conservative CMC formulation are employed for the bluff-body flame with flow recirculation. Based on the numerical results, a detailed discussion is given for the comparative performances of the two combustion models in terms of the flame structure and NO x formation characteristics.  相似文献   

18.
In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification.

All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1?φ?1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate.

In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1?φ?1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4?φ?1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1?φ?1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data, thereby rendering this approach promising for application in LES.  相似文献   

19.
Plane Couette flow, the flow between two parallel planes moving in opposite directions, is an example of wall-bounded flow experiencing a transition to turbulence with an ordered coexistence of turbulent and laminar domains in some range of Reynolds numbers [R g, R t] . When the aspect-ratio is sufficiently large, this coexistence occurs in the form of alternately turbulent and laminar oblique bands. As R goes up trough the upper threshold R t, the bands disappear progressively to leave room to a uniform regime of featureless turbulence. This continuous transition is studied here by means of under-resolved numerical simulations understood as a modelling approach adapted to the long time, large aspect-ratio limit. The state of the system is quantitatively characterised using standard observables (turbulent fraction and turbulence intensity inside the bands). A pair of complex order parameters is defined for the pattern which is further analysed within a standard Ginzburg–Landau formalism. Coefficients of the model turn out to be comparable to those experimentally determined for cylindrical Couette flow.  相似文献   

20.
There is a need to better understand particle size distributions (PSDs) from turbulent flames from a theoretical, practical and even regulatory perspective. Experiments were conducted on a sooting turbulent non-premixed swirled ethylene flame with secondary (dilution) air injection to investigate exhaust and in-burner PSDs measured with a Scanning Mobility Particle Sizer (SMPS) and soot volume fractions (fv) using extinction measurements. The focus was to understand the effect of systematically changing the amount and location of dilution air injection on the PSDs and fv inside the burner and at the exhaust. The PSDs were also compared with planar Laser Induced Incandescence (LII) calibrated against the average fv. LII provides some supplemental information on the relative soot amounts and spatial distribution among the various flow conditions that helps interpret the results. For the flame with no air dilution, fv drops gradually along the centreline of the burner towards the exhaust and the PSD shows a shift from larger particles to smaller. However, with dilution air fv reduces sharply where the dilution jets meet the burner axis. Downstream of the dilution jets fv reduces gradually and the PSDs remain unchanged until the exhaust. At the exhaust, the flame with no air dilution shows significantly more particles with an fv one to two orders of magnitude greater compared to the Cases with dilution. This dataset provides insights into soot spatial and particle size distributions within turbulent flames of relevance to gas turbine combustion with differing dilution parameters and the effect dilution has on the particle size. Additionally, this work measures fv using both ex situ and in situ techniques, and highlights the difficulties associated with comparing results across the two. The results are useful for validating advanced models for turbulent combustion.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号