首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 26 毫秒
1.
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.  相似文献   

2.
Heat release rate in combustion systems must be understood in order to control thermoacoustic instabilities, flame extinction, and heat losses. Traditionally OH chemiluminescence (OH*) is used to trace heat release rate (HRR) in H2/air flames, but its accuracy as a tracer has not been assessed. Lean premixed H2/air cellular tubular flames are a good test case to evaluate HRR tracers due to the presence of highly reactive flame cells surrounded by regions of near extinction. Comparing the calculated heat release rate to OH* concentration, one finds that [OH*] profiles correlate with the regions of high reactivity (flame cells) but the correlation fails in the low reactivity regions where the HRR is much higher than the [OH*] value indicates. Alternate HRR tracers including [H] and pixel-by-pixel products of [O2]x[H], [OH]x[H2], and [O]x[H2] are analyzed with detailed numerical simulations. The chosen products derive from the main chain reaction steps that contribute to overall HRR in lean, premixed H2/air flames. Findings suggest that [H] is an accurate yet simple way of tracking HRR. Planar measurements of HRR are possible if LIF measurements of [H] are improved.  相似文献   

3.
Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.  相似文献   

4.
Thermo-acoustic instabilities are problematic in the design of continuous-combustion propulsion systems such as gas turbine engines, rocket motors, jet engine afterburners, and ramjets. Conceptually, the coupling between acoustics and flame dynamics can be divided into two categories: flame area fluctuations and changes in the local flame speed. The latter can be caused by the thermodynamic fluctuations that accompany an acoustic wave. This coupling is the focus of the present work. In this paper, we are concerned with the dynamics of laminar premixed flames involving large hydrocarbon species. Through high-fidelity numerical simulations, we investigate the flame response for a wide range of fuels and acoustic frequencies. The combustion of hydrogen and methane is considered for verification purposes and as baseline cases for comparison with two large hydrocarbon fuels, n-heptane and n-dodecane. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. For all fuels, we observe a local peak in the heat release gain. At high frequencies, we find that the fluctuations of the different species mass fractions decrease with the inverse of the acoustic frequency, leading to chemistry being “frozen” in the high-frequency limit. This allows us to predict the flame behavior directly from the steady-state solution.  相似文献   

5.
Stationary combustion regimes, their linear stability and extinction limits of stretched premixed flames in a narrow gap between two heat conducting plates are studied by means of numerical simulations in the framework of one-dimensional thermal-diffusion model with overall one-step reaction. Various stationary combustion modes including normal flame (NF), near-stagnation plane flame (NSF), weak flame (WF) and distant flame (DF) are detected and found to be analogous to the same-named regimes of conventional counterflow flames. For the flames stabilized in the vicinity of stagnation plane at moderate and large stretch rates (which are NF, NSF and WF) the effect of channel walls is basically reduced to additional heat loss. For distant flame characterized by large flame separation distance and small stretch rates intensive interphase heat transfer and heat recirculation are typical. It is shown that in mixture content / stretch rate plane the extinction limit curve has ε-shape, while for conventional counterflow flames it is known to be C-shaped. This result is quite in line with recent experimental findings and is explained by extension of extinction limits at small stretch rates at the expense of heat recirculation. Analysis of the numerical results makes possible to reveal prime mechanisms of flame quenching on different branches of ε-shaped extinction limit curve. Namely, two upper limits are caused by stretch and heat loss. These limits are direct analogs of the upper and lower limits on conventional C-shaped curve. Two other limits are related with weakening of heat recirculation and heat dissipation to the burner. Thus, the present study provides a satisfactory explanation for the recent experimental observations of stretched flames in narrow channel.  相似文献   

6.
The structure and extinction of low strain rate nonpremixed methane–air flames was studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of normal gravity (1-g) and zero gravity (0-g) flames. For comparison with the 2D modelling results, a one-dimensional (1D) flamelet computation using a previously developed numerical code was exercised to provide information on the 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. Photographic images of flames undergoing the process of extinction were compared with model calculations. The axisymmetric numerical model was validated by comparing flame shapes, temperature profiles, and extinction limits with experiments and with the 1D computational results. The 2D computations yielded insight into the extinction mode and flame structure. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. The contribution by each term in the energy equation to the heat release rate was evaluated to investigate the multi-dimensional structure and radiative extinction of the 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame and flow structure were also found. At low strain rates in 1-g flames (‘regime A’), the flame is extinguished from the weak outer edge of the flame, which is attributed to a multi-dimensional flame structure and flow field. At high strain rates, (‘regime B’), the flame extinction initiates near the flame centreline owing to an increased diluent concentration in the reaction zone, similar to the extinction mode of 1D flames. These two extinction modes can be clearly explained by consideration of the specific maximum heat release rate.  相似文献   

7.
Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.  相似文献   

8.
Self-excited combustion instabilities in a mesoscale multinozzle array, also referred to as a micromixer-type injector, have been experimentally investigated in a lean-premixed tunable combustor operating with preheated methane and air. The injector assembly consists of sixty identical swirl injectors of 6.5 mm inner diameter, which are evenly distributed across the combustor dump plane. Their flow paths are divided into two groups – inner and outer stages – to form radially stratified reactant stoichiometry for the control of self-excited instabilities. OH PLIF measurements of stable flames reveal that the presence of radial staging has a remarkable influence on stabilization mechanisms, reactant jet penetration/merging, and interactions between adjacent flame fronts. In an inner enrichment case, two outer (leaner) streams merge into a single jet structure, whereas the inner (richer) reactant jets penetrate far downstream without noticeable interactions between neighboring flames. The constructed stability map in the 〈?i, ?o〉 domain indicates that strong self-excited instabilities occur under even split and outer enrichment conditions at relatively high global equivalence ratios. This is attributed to large-scale flame surface deformation in the streamwise direction, as manifested by vigorous detachment/attachment movements. The use of the inner fuel staging method was found, however, to limit the growth of large-amplitude heat release rate fluctuations, because the center flames are securely anchored during the whole period of oscillation, giving rise to a moderate lateral motion. We demonstrate that the collective motion of sixty flames – rather than the individual local flame dynamics – play a central role in the development of limit cycle oscillations. This suggests that the distribution pattern of the injector array, in combination with the radial fuel staging scheme, is the key to the control of the instabilities.  相似文献   

9.
The present experimental investigation demonstrates important trends and offers physical insights into self-excited combustion instabilities in mesoscale multinozzle flames composed of sixty small injectors. Here we focus on the response of a prototypical micromixer-type injector assembly, fabricated using an additive manufacturing technique, in comparison with the behavior of conventional large-scale swirl-stabilized flames. Our results highlight that the development of self-excited instabilities in unconventional mesoscale flames is fundamentally different from that in large-scale swirl flames, in terms of the onset of instabilities, nonlinear modal dynamics, and amplitude/frequency of limit cycle oscillations under the same operating conditions. These differences are attributable to the alteration in local flow/flame structures and the resulting flame-to-flame/flame-wall interaction mechanisms. An integrated analysis of large datasets reveals that the two interacting swirl-stabilized flames tend to couple strongly with a low-frequency L1 mode at about 220 Hz, whereas the sixty-injector small-scale flames are capable of triggering multiple higher-frequency instabilities at ~ 310, ~ 470, and ~ 600 Hz. That is, the use of the micromixer-type injector assembly in a lean-premixed system causes the occurrence of combustion instabilities to shift toward a higher equivalence ratio. However, due to the absence of a large recirculation zone near the primary reaction region, the combustion system equipped with the small-scale multinozzle injectors was found to suffer from lean blowoff phenomena at low equivalence ratio.  相似文献   

10.
A model is presented for a one-dimensional laminar premixed flame, propagating into a rich, off-stoichiometric, fresh homogenous mixture of water-in-fuel emulsion spray, air and inert gas. Due to its relatively large latent heat of vaporisation, the water vapour acts to cool the flame that is sustained by the prior release of fuel vapour. To simplify the inherent complexity that characterises the analytic solution of multi-phase combustion processes, the analysis is restricted to fuel-rich laminar premixed water-in-fuel flames, and assumes a single-step global chemical reaction mechanism. The main purpose is to investigate the steady-state burning velocity and burnt temperature as functions of parameters such as initial water content in the emulsified droplet and total liquid droplet loading. In particular, the influence of micro-explosion of the spray’s droplets on the flame’s characteristics is highlighted for the first time. Steady-state analytical solutions are obtained and the sensitivity of the flame temperature and the flame propagating velocity to the initial water content of the micro-exploding emulsion droplets is established. A linear stability analysis is also performed and reveals the manner in which the micro-explosions influence the neutral stability boundaries of both cellular and pulsating instabilities.  相似文献   

11.
Flame shape is an important observed characteristic of flames that can be used to scale flame properties such as heat release rates and radiation. Flame shape is affected by fuel type, oxygen levels in the oxidiser, inverse burning and gravity. The objective of this study is to understand the effect of high oxygen concentrations, inverse burning, and gravity on the predictions of flame shapes. Flame shapes are obtained from recent analytical models and compared with experimental data for a number of inverse and normal ethane flame configurations with varying oxygen concentrations in the oxidiser and under earth gravity and microgravity conditions. The Roper flame shape model was extended to predict the complete flame shapes of laminar gas jet normal and inverse diffusion flames on round burners. The Spalding model was extended to inverse diffusion flames. The results show that the extended Roper model results in reasonable predictions for all microgravity and earth gravity flames except for enhanced oxygen normal diffusion flames under earth gravity conditions. The results also show trends towards cooler flames in microgravity that are in line with past experimental observations. Some key characteristics of the predicted flame shapes and parameters needed to describe the flame shape using the extended Roper model are discussed.  相似文献   

12.
The detailed flame structure of laminar premixed cellular flames in the tubular domain is simulated in 2D using a fully-implicit primitive variable finite difference formulation that includes multicomponent transport and detailed chemical kinetics. Numerical results for H2/air flames are presented and compared against spatially resolved experimental measurements of temperature and chemical species including atomic H and OH. The experimental results compare well for flame structure and cell number, despite the numerical model under-predicting the peak temperature by 200 K. Numerical experiments were performed to assess the ability for cellular tubular flames to impact experimental and numerical investigations of practical flames. The cellular flame structure is found to provide a highly sensitive geometry that is useful for validating diffusive transport modelling approximations. This capability is exemplified through the development of a simple and accurate approximation for thermal diffusion (i.e. the Soret effect) that is suitable for practical combustion codes.  相似文献   

13.
14.
This work is devoted to the investigation and subgrid-scale modeling of intrinsic flame instabilities occurring in the propagation of a deflagration wave. Such instabilities, of hydrodynamic and thermodiffusive origin, are expected to be of particular relevance in recent technological trends such as in the use of hydrogen as a clean energy carrier or as a secondary fuel in hydrogen enriched combustion. A dedicated set of direct numerical simulations is presented and used, in conjunction with coherent literature results, in order to develop scaling arguments for the propagation speed of self-wrinkled flames which are also supported by the outcomes of a weakly non-linear model, namely the Sivashinsky equation. The observed scaling is based on the definition of the number of unstable wavelengths in a reference hydrodynamic lengthscale, in other words the ratio between the neutral or cutoff lengthscale of intrinsic instabilities and the lateral domain of a planar flame. The scalings are then employed to develop an algebraic model for the wrinkling factor in the context of a flame surface density closure approach. An a-priori analysis shows that the model correctly captures the flame wrinkling caused by intrinsic instability at sub grid level. A strategy to include the developed self-wrinkling model in the context of a turbulent combustion model is finally discussed on the basis of the turbulence induced cut-off concept.  相似文献   

15.
Pilot-ignited dual fuel combustion involves a complex transition between the pilot fuel autoignition and the premixed-like phase of combustion, which is challenging for experimental measurement and numerical modelling, and not sufficiently explored. To further understand the fundamentals of the dual fuel ignition processes, the transient ignition and subsequent flame development in a turbulent dimethyl ether (DME)/methane-air mixing layer under diesel engine-relevant conditions are studied by direct numerical simulations (DNS). Results indicate that combustion is initiated by a two-stage autoignition that involves both low-temperature and high-temperature chemistry. The first stage autoignition is initiated at the stoichiometric mixture, and then the ignition front propagates against the mixture fraction gradient into rich mixtures and eventually forms a diffusively-supported cool flame. The second stage ignition kernels are spatially distributed around the most reactive mixture fraction with a low scalar dissipation rate. Multiple triple flames are established and propagate along the stoichiometric mixture, which is proven to play an essential role in the flame developing process. The edge flames gradually get close to each other with their branches eventually connected. It is the leading lean premixed branch that initiates the steady propagating methane-air flame. The time required for the initiation of steady flame is substantially shorter than the autoignition delay time of the methane-air mixture under the same thermochemical condition. Temporal evolution of the displacement speed at the flame front is also investigated to clarify the propagation characteristics of the combustion waves. Cool flame and propagation of triple flames are also identified in this study, which are novel features of the pilot-ignited dual fuel combustion.  相似文献   

16.
Understanding the distinguishing physical properties of multi-element lean-premixed high hydrogen content flames is expected to be integral to the development of carbon-neutral, and ultimately carbon-free, gas turbine combustion systems. Despite their fundamental importance, the thermoacoustic and emission-related characteristics of such small-scale flame ensembles are not thoroughly understood, particularly for the full range of 0 to 100% hydrogen content blended with methane fuel. Here we investigate the structure and collective behavior of a multi-element lean-premixed hydrogen/methane/air flame ensemble using measurements of nitrogen oxides emissions and self-excited instability, combined with OH* and OH PLIF flame visualizations. Our results indicate that the system's responses can be classified into several distinctive stages according to their static and dynamic stability, including flame blowoff and thermoacoustically stable regions under relatively low hydrogen concentration conditions, low-frequency self-excited instabilities in intermediate hydrogen concentration, and triggering of intense pressure perturbations at about 1.7 kHz under high- or pure hydrogen combustion conditions. While the low-frequency combustion dynamics are dominated by axisymmetric translational movements of parallel flame fronts, the higher frequency response originates from significant lateral modulations accompanied by small-scale vortical rollup and flame surface annihilation due to front merging and pinch-off. Longitudinal-to-transverse dynamic transition is observed to play a mechanistic role in kinematically accommodating higher-frequency heat release rate fluctuations, and this newly identified mechanism suggests the possibility of high-frequency transverse modes, if such lateral motions are strong enough to induce inter-element flame interactions. In contrast to the substantial differences in thermoacoustic properties for different fuel compositions, the total nitrogen oxides emissions are found to depend primarily on adiabatic flame temperature; the influence of fuel composition is limited to approximately 20% under the inlet conditions considered.  相似文献   

17.
Given the experimental difficulties, most of the available flame speed database is for relatively reduced thermodynamic conditions and for non-simultaneous variations of pressure and temperature. This limitation may be overpassed by using spherically expanding flames with the constant volume method. This methodology, introduced in the 30 s by Lewis and von Elbe, requires the knowledge of the pressure evolution in the combustion chamber. It has been penalized for a long time because of the underlying assumptions and problems in flame instability detection. This method has been greatly renewed recently by Egolfopoulos following a coupled experimental/numerical approach integrating the effects of radiation and dissociation while maintaining moderate computing costs. In parallel with this study, we have worked on an alternative method providing a maximum of information for each test minimizing uncertainties. The current study uses a new experimental device allowing simultaneous recording of pressure and flame radius inside the chamber during the full combustion process. The direct use of these data over the whole flame propagation allows testing kinetic schemes over large pressure and temperature domains with good accuracy. These new experimental targets allowed the identification of key reactions needing improvements.  相似文献   

18.
This paper describes an experimental investigation of the feasibility of using “slow” active control approaches, which “instantaneously” change liquid fuel spray properties, to suppress combustion instabilities. The objective of this control approach was to break up the feedback between the combustion process heat release and combustor pressure oscillations that drive the instability by changing the characteristics of the combustion process (e.g., the characteristic combustion time). To demonstrate the feasibility of such control, this study used a proprietary fuel injector (NanomiserTM), which can vary its fuel spray properties, to investigate the dependence of acoustics–combustion process coupling, i.e., the driving of combustion instabilities, upon the fuel spray properties. This study showed that by changing the spray characteristics it is possible to significantly damp combustion instabilities. Furthermore, using combustion zone chemiluminescence distributions, which were obtained by Abel’s deconvolution synchronized with measured acoustic data, it has been shown that the instabilities were mostly driven midway between the combustor centerline and wall, a short distance downstream from the flame holder, where the mean axial flow velocity is approximately zero in the vortex near the flame holder. The results of this study strongly suggest that a “slow” active control system that employs controllable fuel injectors could be effectively used to prevent the onset of detrimental combustion instabilities.  相似文献   

19.
A numerical investigation of the interaction between a spray flame and an acoustic forcing of the velocity field is presented in this paper. In combustion systems, a thermoacoustic instability is the result of a process of coupling between oscillations in heat released and acoustic waves. When liquid fuels are used, the atomisation and the evaporation process also undergo the effects of such instabilities, and the computational fluid dynamics of these complex phenomena becomes a challenging task. In this paper, an acoustic perturbation is applied to the mass flow of the gas phase at the inlet and its effect on the evaporating fuel spray and on the flame front is investigated with unsteady Reynolds averaged Navier-Stokes numerical simulations. Two flames are simulated: a partially premixed ethanol/air spray flame and a premixed pre-vaporised ethanol/air flame, with and without acoustic forcing. The frequencies used to perturb the flames are 200 and 2500 Hz, which are representative for two different regimes. Those regimes are classified based on the Strouhal number St = (D/U)ff: at 200 Hz, St = 0.07, and at 2500 Hz, St = 0.8. The exposure of the flame to a 200 Hz signal results in a stretching of the flame which causes gas field fluctuations, a delay of the evaporation and an increase of the reaction rate. The coupling between the flame and the flow excitation is such that the flame breaks up periodically. At 2500 Hz, the evaporation rate increases but the response of the gas field is weak and the flame is more stable. The presence of droplets does not play a crucial role at 2500 Hz, as shown by a comparison of the discrete flame function in the case of spray and pre-vaporised flame. At low Strouhal number, the forced response of the pre-vaporised flame is much higher compared to that of the spray flame.  相似文献   

20.
Combustion instabilities in annular combustors are of great interest because of their industrial relevance. Azimuthal acoustic modes, which involve transverse acoustic forcing to flames, have become a key process related to annular combustor instabilities. Transverse mean flow may be a factor that affects azimuthal oscillations. This paper provides an analytical model for a transversely forced two-dimensional Bunsen flame under transverse mean flow. The model is established using a low-amplitude perturbation assumption applied to a G-equation formulation. Forced flame displacement and flame transfer functions (FTFs) are calculated. The results are verified based on numerical solutions of the G-equation. Effects of frequency, transverse mean flow velocity and vertical mean flow velocity on the FTFs are discussed. The symmetric flame without transverse mean flow has a vanishing response to transverse acoustic forcing, while asymmetric flames, which are formed with transverse mean flow, have a bandpass response to transverse forcing. The response at very low and high forcing frequencies is small, with higher transfer function gains only in a certain frequency range. This bandpass response, which is inherently linked to the asymmetry of the flame, is an important factor to account for when considering the flame dynamics related to transverse acoustic effects.  相似文献   

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

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