单只水平浓缩煤粉燃烧器在1MW燃烧试验台上的试验研究

单只水平浓缩煤粉燃烧器在１ＭＷ燃烧试验台上的试验研究孙绍增，吴少华，李争起，杨明新，王新雷，陈力哲，庞丽君，邢春礼，朱彤，孙恩召，秦裕琨（哈尔滨工业大学动力工程系哈尔滨１５０００１）关键词：水平浓淡燃烧，煤粉燃烧器，稳燃，低ＮＯ_ｘ一、引言电力工业对?..     

Steam-enhanced regime for liquid hydrocarbons combustion: velocity distribution in the burner flame

The lab-scale burner device with proprietary design was used for combustion of diesel fuel in a steam-enhanced regime. This operation mode ensures drastic intensification of liquid hydrocarbon combustion due to supply of superheated steam jet to the combustion zone. The particle image velocimetry technique was used for study of velocity field in the burner flame. The method of seeding of flow zone with new kind of tracers (micro-sized silica particles produced from silicon oil added to liquid fuel) was tested.     

Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Combustion characteristics of a dual-mode scramjet combustor with cavity flameholder

Combustion characteristics of a laboratory dual-mode ramjet/scramjet combustor were studied experimentally. The combustor consists of a sonic fuel jet injected into a supersonic crossflow upstream of a wall cavity pilot flame. These fundamental components are contained in many dual-mode combustor designs. Experiments were performed with an isolator entrance Mach number of 2.2. Air stagnation temperatures were varied from 1040 to 1490 K, which correspond to flight Mach numbers of 4.3–5.4. Both pure hydrogen and a mixture of hydrogen and ethylene fuels were used. High speed imaging of the flame luminosity was performed along with measurements of the isolator and combustor wall pressures. For ramjet mode operation, two distinct combustion stabilization locations were found for fuel injection a sufficient distance upstream of the cavity. At low T₀, the combustion was anchored at the leading edge of the cavity by heat release in the cavity shear layer. At high T₀, the combustion was stabilized a short distance downstream of the fuel injection jet in the jet-wake. For an intermediate range of T₀, the reaction zone oscillated between the jet-wake and cavity stabilization locations. Wall pressure measurements showed that cavity stabilized combustion was the steadiest, followed by jet-wake stabilized, and the oscillatory case. For fuel injection close to the cavity, a hybrid stabilization mode was found in which the reaction zone locations for the two stabilization modes overlapped. For this hybrid stabilization, cavity fueling rate was an important factor in the steadiness of the flow field. Scramjet mode combustion was found to only exist in the cavity stabilized location for the conditions studied.     

Large Eddy Simulations of a piloted lean premix jet flame using finite-rate chemistry

A Large Eddy Simulation (LES) model capable of accurately representing finite-rate chemistry effects in turbulent premixed combustion is presented. The LES computations use finite-rate chemistry and implicit LES combustion modelling to simulate an experimentally well-documented lean-premixed jet flame stabilized by a stoichiometric pilot. The validity of the implicit LES assumption is discussed and criteria are expressed in terms of subgrid scale Damköhler and Karlovitz numbers. Simulation results are compared to experimental data for velocity, temperature and species mass fractions of CH₄, CO and OH. The simulation results highlight the validity and capability of the present approach for the flame and in general the combustion regime examined. A sensitivity analysis to the choice of the finite-rate chemistry mechanism is reported, this analysis indicates that the one and two-step global reaction mechanisms evaluated fail to capture the reaction layer with sufficient accuracy, while a 20-species skeletal mechanism reproduces the experimental observations accurately including the key finite-rate chemistry indicators CO and OH. The LES results are shown to be grid insensitive and that the grid resolution within the bounds examined is far less important compared to the sensitivity of the finite-rate chemistry representation. The results are analyzed in terms of the flame dynamics and it is shown that intense small scale mixing (high Karlovitz number) between the pilot and the jet is an important mechanism for the stabilization of the flame.     

气体回流区分级着火试验研究

就回流区分级着火燃烧方式在气体燃料燃烧上的应用进行了一系列冷热态模拟试验和工业试验研究．比较了直流、钝体和开缝钝体喷口的出口速度、温度分布，以及在不同气体燃料热值时的稳焰能力。试验表明，开缝钝体喷口中缝小股射流进入喷口后低速回流区能稳定着火，进一步点燃主流，有效地形成回流区分级着火；开维钝体喷口的稳焰能力最强，特别适合低热值煤气燃烧。在大容量锅炉上进行了低热值高炉煤气燃烧的改造性工业试验，试验证明，开缝钝体燃烧器强化了低热值煤气的燃烧，有效地解决了锅炉煤气段燃烧强度不足以至锅炉尾部超温问题。     

Combustion of partially premixed spray jets

Gas turbines, liquid rocket motors, and oil-fired furnaces utilize the spray combustion of continuously injected liquid fuels. In most cases, the liquid spray is mixed with an oxidizer prior to combustion, and further oxidizer is supplied from the outside of the spray to complete diffusion combustion. This rich premixed spray is called “partially premixed spray.” Partially premixed sprays have not been studied systematically although they are of practical importance. In the present study, the burning behavior of partially premixed sprays was experimentally studied with a newly developed spray burner. A fuel spray and an oxidizer, diluted with nitrogen, was injected into the air. The overall equivalence ratio of the spray jet was set larger than unity to establish partially premixed spray combustion. In the present burner, the mean droplet diameter of the atomized liquid fuel could be varied without varying the overall equivalence ratio of the spray jet. Two combustion modes with and without an internal flame were observed. As the mean droplet diameter was increased or the overall equivalence ratio of the spray jet was decreased, the transition from spray combustion only with an external group flame to that with the internal premixed flame occurred. The results suggest that the internal flame was supported by flammable mixture through the vaporization of fine droplets, and the passage of droplet clusters deformed the internal flame and caused internal flame oscillation. The existence of the internal premixed flame enhanced the vaporization of droplets in the post-premixed-flame zone within the external diffusion flame.     

Flame characterisation of gas-assisted pulverised coal combustion using FPV-LES

A multiphase flamelet/progress variable (FPV) model for the large eddy simulation (LES) of gas-assisted pulverised coal combustion (PCC) is developed. The target of the simulation is the Darmstadt turbulent gas-assisted swirling solid fuel combustion chamber. The coal particles are treated as Lagrangian point particles, the position, momentum and energy of which are tracked. The gas phase is described by the low-Mach Navier-Stokes equations alongside the Eulerian transport equations of the governing variables for the FPV model. The set of chemical states of the PCC flame is pre-tabulated in a six-dimensional flamelet table and determined by the mixing of the primary fuel stream, volatiles and char off-gases with the oxidising air, the progress of chemical reactions, the interphase heat transfer, as well as sub-grid scale variations. A presumed $\beta$-PDF approach for the total mixture fraction is applied to capture sub-grid scale effects. The discrete ordinate method (DOM) with the weighted sum of grey gases model (WSGGM) is employed to model radiation. The FPV-LES results are validated against the experimental evidence and a good agreement of the predicted mean and RMS velocities, as well as the mean gas temperature between experiments and simulations is obtained. The contributions of the pilot, volatile and char off-gas fuel streams to the coal flame are analysed. It is found that most regions of the furnace are dominated by either pilot or volatile combustion, while char conversion only occurs in the far downstream and outer furnace regions. The pilot gas dominates the near-wall region inside the quarl, whereas the volatile gas mainly released from small particles dominates a first volatile combustion zone in the interior of the internal recirculation zone. Larger particles heat up more slowly and release their volatile content further downstream, leading to a secondary volatile combustion zone.     

The effect of fuel inlet turbulence intensity on H2/CH4 flame structure of MILD combustion using the LES method

Large eddy simulations (LES) are employed to investigate the effect of the inlet turbulence intensity on the H₂/CH₄ flame structure in a hot and diluted co-flow stream which emulates the (Moderate or Intense Low-oxygen Dilution) MILD combustion regime. In this regard, three fuel inlet turbulence intensity profiles with the values of 4%, 7% and 10% are superimposed on the annular mixing layer. The effects of these changes on the flame structure under the MILD condition are studied for two oxygen concentrations of 3% and 9% (by mass) in the oxidiser stream and three hot co-flow temperatures 1300, 1500 and 1750 K. The turbulence-chemistry interaction of the numerically unresolved scales is modelled using the (Partially Stirred Reactor) PaSR method, where the full mechanism of GRI-2.11 represents the chemical reactions. The influences of the turbulence intensity on the flame structure under the MILD condition are studied by using the profile of temperature, CO and OH mass fractions in both physical and mixture fraction spaces at two downstream locations. Also, the effects of this parameter are investigated by contours of OH, HCO and CH₂O radicals in an area near the nozzle exit zone. Results show that increasing the fuel inlet turbulence intensity has a profound effect on the flame structure particularly at low oxygen mass fraction. This increment weakens the combustion zone and results in a decrease in the peak values of the flame temperature and OH and CO mass fractions. Furthermore, increasing the inlet turbulence intensity decreases the flame thickness, and increases the MILD flame instability and diffusion of un-burnt fuel through the flame front. These effects are reduced by increasing the hot co-flow temperature which reinforces the reaction zone.     

The impact of N2 micro-jets on the V-to-M flame shape transition in a premixed swirl burner

Injection of N₂ through micro-jets located on the dump plane of a lean premixed swirl stabilized combustor is investigated as a new method for mitigating combustion instabilities. This study focuses on the chemical and fluid dynamic processes by which the N₂ micro-jets impact the flame dynamics. An experimental and numerical investigation is performed to characterize the combustion instability during the V-to-M flame shape transition in a swirl burner fueled with premixed CH₄/air, at an equivalence ratio of 0.62. Reasonable agreements have been found between the experimental measurements and simulation results. Both of them present that the flame changes from V-shape to M-shape periodically, and a low-frequency instability around 10 Hz is observed accordingly. It is confirmed that intermittent flame extinction in the outer recirculation zone (ORZ) is the source of the combustion instability. Furthermore, injection of N₂ through micro-jets located on the combustor dump plane, into the outer recirculation zone, results in a stable V shape flame. It is clearly seen that the ORZ dilution can eliminate the combustion instability without inhibiting the combustion efficiency. A special focus is placed on the impact of the diluent injection on the local flame-flow interaction. The nitrogen micro-jets increase the local nitrogen concentration by 7% on average, lowering the flame speed and extinction strain rates by 27% and 17% respectively. Moreover, the micro-jets increase the turbulence intensity in the ORZ, leading to a significant increase in the Karlovitz number and transferring the local combustion regime from the thin reaction zone regime to the broken reaction zone regime. Hence, the nitrogen micro-jets impact on both the turbulence and the chemical reaction rates prevents flame propagation into the ORZ and results in a stable flame.     

Laminar flame propagation on a horizontal fuel surface: Verification of classical Emmons solution

This work analyses the classical Emmons (1956) solution of flat plate laminar flame combustion on a film of liquid fuel. A two-dimensional (2D) numerical model developed for this purpose has been benchmarked with experimental results available in the literature for methanol. In the parametric study, numerical predictions have been compared with Emmons classical solution. The study shows that the Emmons solution is valid in a range of Reynolds numbers where flame anchors near the leading edge of the methanol pool and the combustion zone is confined around the hydrodynamic and thermal boundary layers. However, in cases of low free stream velocities the combustion zone is beyond the boundary layer zone and the Emmons solution deviates. In cases of very high free stream velocities, the flame moves away from the leading edge and anchors at a location downstream. The Emmons solution is not applicable in this case as well. For the fuel considered in this study (methanol), accounting for thermal radiation, employing an optically thin radiation model, allows better agreement between experimental and numerical temperature profiles but does not affect the mass burning rates.     

An investigation of oscillating instabilities in deflagrations with edges

We examine the Lewis-number-greater-than-1 stability of a deflagration sitting on a porous-plug burner with an inert coflow. The flame edges generated by the coflow influence the stability, and this influence is examined. Very wide flames display the same stability characteristics as unbounded flames (flames sans edges), but for moderately wide flames the instability is suppressed. A new two-dimensional instability can occur for narrow flames. There is a range of mass fluxes for which a monotonic decrease in burner (flame) width generates a transition from unstable flames to stable flames, to unstable flames, to quenching. The insertion of a cold probe into the combustion field can stabilize an unstable flame or destabilize a stable flame, depending on the point of insertion.     

Effect of fuel mixture fraction and velocity perturbations on the flame transfer function of swirl stabilized flames

A novel methodology is developed to decompose the classic Flame Transfer Function (FTF) used in the thermo-acoustic stability analysis of lean premix combustors into contributions of different types. The approach is applied, in the context of Large Eddy Simulation (LES), to partially-premixed and fully-premixed flames, which are stabilized via a central recirculation zone as a result of the vortex breakdown phenomenon. The first type of decomposition is into contributions driven by fuel mixture fraction and dynamic velocity fluctuations. Each of these two contributions is further split into the components of turbulent flame speed and flame surface area. The flame surface area component, driven by the pure dynamic velocity fluctuation, which is shown to be a dominant contribution to the overall FTF, is also additionally decomposed over the coherent flow structures using proper orthogonal decomposition. Using a simplified model for the dynamic response of premixed flames, it is shown that the distribution of the FTF, as obtained from LES, is closely related to the characteristics of the velocity field frequency response to the inlet perturbation. Initially, the proposed method is tested and validated with a well characterized laboratory burner geometry. Subsequently, the method is applied to an industrial gas turbine burner.     

Experimental and numerical study of a conical turbulent partially premixed flame

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 CH₂O 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 CH₂O 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 CH₂O structure is somewhat thinner than the experiments.     

双层多孔介质燃烧器的数值模拟

考虑多孔固体构架的辐射换热以及气固相间的对流换热，引入弥散效应及相间对流换热系数，使用GRI3．0机理和双通量辐射模型，数值求解双层多孔介质内燃烧过程．分析了双层多孔介质燃烧器内火焰稳定性和污染物排放，并与单层的进行比较．结果表明，双层多孔介质燃烧器能够在较宽的流量范围内将火焰稳定在它的交界面附近．     

NH3/H2预混旋流火焰稳定性及燃烧极限实验研究

采用叶轮型旋流燃烧器,选取氢气作为燃料添加剂,研究了掺氢比对氨气旋流火焰稳定性的影响,分析了不同旋流数、叶片数、当量比以及预混气总流量条件下,旋流火焰形态变化。测定并分析了不同参数对旋流火焰燃烧极限范围的影响。结果表明,随掺氢比的增大,火焰逐渐由“V”型转化为稳定的“M”型,燃烧反应愈发充分;高旋流数(1.27)或低叶片数(6片)相比低旋流数(0.42)或高叶片数(8片)更有利于旋流火焰的稳定和燃烧的充分进行;相比富燃,贫燃有利于形成稳定的旋流火焰;预混气总流量较大时,火焰高度较高.对于燃烧极限,掺氢比越高,极限范围越大;总流量的变化对贫燃极限影响较小,对富燃极限影响较大;高旋流数(1.27)条件下,燃烧极限范围较大。     

Vortex dynamics in different combustion regions of a cavity-based scramjet

A hybrid RANS/LES study of a cavity-based scramjet was performed and reasonable agreements were found between simulation results and experimental measurements. In the current case, the flame was stabilized by the subsonic cavity shear layer and propagated downstream into the supersonic flow. The vortex dynamic in the flow, mixing, and combustion regions was comparatively investigated. The averaged vorticity in the combustion regions was lower by 55% compared to the mixing region, primarily due to dilatation as a result of the heat release. Furthermore, the combustion zone was decomposed into four regions based on premixed/diffusion flame and subsonic/supersonic combustion. Then the vorticity and its transport in the four regions were compared. The averaged vorticity in the premixed combustion regions was only slightly larger than that in the diffusion combustion regions. However, the averaged heat release rate was nearly 3 times larger in the premixed regions, leading to higher contributions of dilatation and baroclinic torque in the premixed regions, with an overall weak positive impact on the vorticity generation. In the subsonic combustion regions, the vorticity was three times larger than that in the supersonic combustion regions, despite similar heat release rates on average. It could be explained by the relatively large magnitude of dilatation and baroclinic torque in the supersonic flow. Vortex stretching and dilatation were comparable in the supersonic flame but the former became two times larger than the latter in the subsonic flame. Moreover, the baroclinic torque had larger contributions than diffusion in the supersonic flame whereas the opposite trend was found in the subsonic flame. The results highlight that the subsonic combustion regions in the cavity shear layer and near the walls significantly contribute to the vortex dynamics and mixing process, in addition to flame stabilization.     

Large-eddy simulations of a stratified-charge direct-injection spark-ignition engine: Comparison with experiment and analysis of cycle-to-cycle variations

Multiple-cycle large-eddy simulations (LES) have been performed for an optically accessible, single-cylinder, four-stroke-cycle, spray-guided direct-injection spark-ignition (SG-DISI) engine operating in a stratified globally fuel-lean mode. The simulations combine a standard Smagorinsky turbulence model, a stochastic Lagrangian parcel method for liquid fuel injection and fuel spray modeling, a simple energy-deposition spark-ignition model, and a modified thickened flame model for turbulent flame propagation through highly stratified reactant mixtures. Comparisons between simulations and experiments include individual-cycle and ensemble-average pressure and apparent-heat-release-rate traces, individual-cycle and ensemble-average indicated mean effective pressures (IMEP), and instantaneous two-dimensional vapor-equivalence-ratio contours. Although the number of LES cycles is small (35), the results show that the simulations are able to capture the global combustion behavior that is observed in the experiments, including cycle-to-cycle variations. The simulation results are then analyzed further to provide insight into the conditions that lead to misfire versus robust combustion. As has been reported in earlier experimental and LES studies for homogeneous-charge SI engines, local conditions in the vicinity of the spark gap at the time of ignition largely determine the subsequent flame development. However, in contrast to homogeneous-charge engines, no single local or global quantity correlates as strongly with the eventual peak pressure or IMEP for each cycle. Rather, it is the interplay among the early flame kernel, the velocity field that it experiences, and the fuel distribution that it encounters that ultimately determines the fate of each combustion event. Deeper analysis and quantitative statistical comparisons between experiments and simulations will require the simulation of larger numbers of engine cycles.     

Soot formation and temperature field structure in co-flow laminar methane-air diffusion flames at pressures from 10 to 60 atm

The effects of pressure on soot formation and the structure of the temperature field were studied in co-flow methane-air laminar diffusion flames over a wide pressure range, from 10 to 60 atm in a high-pressure combustion chamber. The selected fuel mass flow rate provided diffusion flames in which the soot was completely oxidized within the visible flame envelope and the flame was stable at all pressures considered. The spatially resolved soot volume fraction and soot temperature were measured by spectral soot emission as a function of pressure. The visible (luminous) flame height remained almost unchanged from 10 to 100 atm. Peak soot concentrations showed a strong dependence on pressure at relatively lower pressures; but this dependence got weaker as the pressure is increased. The maximum conversion of the fuel’s carbon to soot, 12.6%, was observed at 60 atm at approximately the mid-height of the flame. Radial temperature gradients within the flame increased with pressure and decreased with flame height above the burner rim. Higher radial temperature gradients near the burner exit at higher pressures mean that the thermal diffusion from the hot regions of the flame towards the flame centerline is enhanced. This leads to higher fuel pyrolysis rates causing accelerated soot nucleation and growth as the pressure increases.     

双层多孔介质内甲烷富燃制氢过程的研究

采用一维双温度体积平均模型和详细的甲烷化学反应机理GRI 3.0,对双层泡沫陶瓷多孔介质内甲烷富燃燃烧过程进行数值模拟,研究在双层多孔介质交界面附近稳定燃烧时的火焰稳定传播范围、火焰温度和组分分布及氢气的产量和能量转换效率.结果表明,双层多孔介质燃烧器能有效拓宽甲烷在空气中的富燃极限;在当量比大于1.6时,燃烧产物中氢气含量较多,氢气产生分为甲烷部分氧化和水煤气反应两个阶段;当量比在1.6～1,8之间时,能量转换效率较大,最大值约为46%.