Effects of nitrogen and argon on ammonia-oxygen explosion

The primary task of this work is to clarify ammonia-oxygen explosion characteristics under nitrogen and argon atmosphere. Firstly, the flame behavior and explosion pressure are experimentally obtained. Then their correlation is revealed quantitatively. The thermal, diffusive and chemical analysis is conducted at last. The results demonstrated that the variation tendency of flame propagation velocity (FPV), maximum explosion pressure (MEP) and maximum rate of pressure rise (MRPR) is completely consistent. All index of FPV, MEP and MRPR, becomes increased and decreased with increasing equivalence ratio, continues to decrease with increasing inert gas fraction. All index of FPV, MEP and MRPR under argon atmosphere is totally larger than that under nitrogen atmosphere. By considering the state equation of ideal gas, spherically smooth flame and adiabatic compression, the flame behavior and explosion pressure under nitrogen and argon atmosphere is significantly controlled by laminar burning speed (LBS). As the inert gas fraction and equivalence ratio change, the LBS under nitrogen and argon atmosphere is significantly controlled by adiabatic flame temperature. The joint action of adiabatic flame temperature and thermal diffusivity contributes to the LBS difference.     

Measurements of fuel mixture fraction oscillations of a turbulent jet non-premixed flame

This work describes new type of combustion instability for which the 3-way coupling between mixing, flame heat release, and acoustics is modified by local buoyancy effects. Measurements of fuel mixture fraction are made for a non-premixed jet flame in a combustion chamber to assess the dynamics of mixing under imposed acoustic oscillations (22-55 Hz). Infrared laser absorption and phase resolved acetone-planar laser induced fluorescence are used to measure the fuel mixture fraction and then the degree of fuel/air mixing is calculated by determining the unmixedness. Results show acoustic excitation causes oscillations in the degree of fuel/air mixing at the driving frequency, which results in oscillatory flame behavior. This oscillatory flame behavior couples to the buoyancy and this in turn affects the mixing. Results also show that the mixing becomes less effective when the excitation frequency is increased or when the flame is present, compared to the non-reacting case. This work describes a key coupling mechanism that occurs when buoyancy is a significant factor in the flow field.     

Stabilization performances and mechanisms of a diffusion-like vortex-tube combustor for oxygen-enriched combustion

The stabilization performance and mechanisms in a diffusion-like vortex-tube combustor is investigated for oxygen-enriched combustion. The stability limit, flame configuration, and pressure fluctuation are investigated under various conditions. Results show that a diffusion-like flame structure is established in the combustor and the nonpremixed peculiarity becomes more prominent with the increase of oxygen mole fraction. The steady combustion can be achieved in the range of global equivalence ratio 0.01 to 1.0 with a low-pressure fluctuation amplitude always less than 1300 Pa, indicating a good combustion stability of this combustor. Additionally, the stabilization mechanism is discussed from the time matching and velocity matching. Based on the axial fuel entry method, the Damköhler number (Da) is always less than 1.0 as a whole, which is the principal reason for the tubular flame shape and the steady combustion procedure in this vortex-tube combustor. The intensified combustion under oxygen-enriched combustion can increase the flame speed, and subsequently reduce the mixing quality and make the yellow flame more visible. Besides, the temperature distribution and the flow field structure can explain the corrugation and deformation of the flame front under oxygen-enriched conditions.     

Hydrogen–hydrocarbon turbulent non-premixed flame structure

In this study, the structure of turbulent non-premixed CH₄-H₂/air flames is analyzed with a special emphasis on mixing and air entrainment. The amount of H₂ in the fuel mixture varies under constant volumetric fuel flow. Mixing is described by mixture fraction and its variance while air entrainment is characterized by the ratio of gas mass flow to fuel mass flow at the inlet section. The flow field and the chemistry are coupled by the flamelet assumption. Mixture fraction and its variance are transported by the computational fluid dynamics (CFD) code. The slow chemistry aspect of NO_x is handled by solving an additional transport equation with a source term derived from flamelet library.     

Experimental study on explosion characteristics of syngas with different ignition positions and hydrogen fraction

The influence of different ignition positions and hydrogen volume fractions on the explosion characteristics of syngas is studied in a rectangular half-open tube. Three ignition positions were set at the axis of the tube, which are 0 mm, 600 mm and 1100 mm away from the closed end, respectively. A range of hydrogen volume fraction (φ) from 10% to 90% were concerned. Experimental results show that different ignition positions and hydrogen volume fraction have important influence on flame propagation structure. When ignited at 600 mm from the closed end on the tube axis, distorted tulip flame forms when flame propagates to the closed end. The formations of the tulip flame and the distorted tulip flame are accompanied by a change in the direction of the flame front propagation. The flame propagation structure and pressure are largely affected by the ignition position and the hydrogen volume fraction. At the same ignition position, flame propagation speed increases with the growing of hydrogen volume fraction. And the pressure oscillates more severe as the ignition location is closer to the open end. And pressure oscillations bring two different forms. The first form is that the pressure has a periodic oscillation. The amplitude of the pressure oscillation gradually increases. It takes several cycles from the start of the oscillation to the peak. For the second form, the pressure reaches the peak of the oscillation in the first cycle of the start to the oscillation.     

Soot formation, spray characteristics, and structure of jet spray flames under high pressure

Coaxial jet spray flames of kerosene and oxygen are experimentally studied over a pressure range of 0.1–1.0 MPa to determine the relationship between flame structure, droplet behavior, and soot formation region, which varies with changes in pressure. The direct images and chemiluminescence spectra show that the spray flames have three regions: the blue flame region, which has a peak of CH* and C₂* radical chemiluminescence, luminous flame region caused by soot emission, and blue emission region caused by CO₂ emission. With increase in ambient pressure, the flame length shortens drastically, the luminous flame region envelopes the blue flame region, and the blue emission becomes more intense. The result of phase-Doppler anemometry shows that a large number of small droplets evaporate and disappear near the burner, and the evaporation of large droplets also occurs rapidly under high pressure. The result of temperature measurements shows that high-temperature regions appear near the burner. The flame temperature drastically decreases along the central axis, and a minimum temperature point appears. This point moves upstream with increase in ambient pressure because evaporation of the droplets occurs further upstream. A laser-induced incandescence measurement shows that the soot volume fraction does not monotonously increase or decrease with increase in ambient pressure. The soot volume fraction at the central axis becomes low upstream and high downstream. As pressure increases, the vertical position at which the peak of soot volume fraction appears at the central axis moves upstream.     

Mixing and stabilization study of a partially premixed swirling flame using laser induced fluorescence

A laboratory-scale swirling burner, presenting many similarities with gas turbines combustors, has been studied experimentally using planar laser induced fluorescence (PLIF) on OH radical and acetone vapor in order to characterize the flame stabilization process. These diagnostics show that the stabilization point rotates in the combustion chamber and that air and fuel mixing is not complete at the end of the mixing tube. Fuel mass fraction decays exponentially along the mixing tube axis and transverse profiles show a gaussian shape. However, radial pressure gradients tend to trap the fuel in the core of the vortex that propagates axially in the mixing tube. As the mixing tube vortex enters the combustion chamber, vortex breakdown occurs through a precessing vortex core (PVC). The axially propagating vortex shows a helicoidal trajectory in the combustion chamber which trace is observed with transverse acetone PLIF. As a consequence, the stabilizing point of the flame in the combustion chamber rotates with the PVC structure. This phenomenon has been observed in the present study with a high speed camera recording spontaneous emission of the flame. The stabilization point rotation frequency tends to increase with mass flow rates. It was also shown that the coupling between the PVC and the flame stabilization occurs via mixing, explaining one possible coupling mechanism between acoustic waves in the flow and the reaction rate. This path may also be envisaged for flashback, an issue that will be more completely treated in a near future.     

初始条件对燃气轮机DLE燃烧室性能影响研究

以单头部中心分级旋流干式低排放(Dry Low Emission, DLE)燃烧室为研究对象,以天然气为燃料,针对不同的全局当量比、进口温度、进口压力条件开展试验测试和数值模拟,研究燃烧室的燃烧性能以及污染物排放的变化规律。研究发现：随全局当量比增大,中心回流区长度略有增大、宽度变窄、回流速度增大,燃料量的增加使得高温区面积明显扩大,燃烧室出口温升明显增大,出口温度分布系数变化不大,燃烧室出口CO和NO_x排放摩尔分数明显增大;随进口温度的增大,中心回流区长度先明显增大再减小、宽度变窄、回流速度先增大再减小,进口空气温度的升高使得反应速率加快从而导致燃烧室出口温度升高,但温升、出口温度分布系数变化不大,CO和NO_x排放摩尔分数增大;随进口压力的增大,中心回流区长度、宽度略有增大,回流速度增大,燃烧室内部和燃烧室出口温度无明显变化,出口温度分布系数减小,CO和NO_x排放摩尔分数受影响较小。     

Characterization of knocking combustion in HCCI DME engine using wavelet packet transform

A diesel engine is modified for homogeneous charge compression ignition (HCCI) combustion with dimethyl ether. With and without knock, in-cylinder pressure is acquired, and in-cylinder temperature, rate of heat release (ROHR), pressure rise rate and pressure rise acceleration obtained. Wavelet packet transform is performed to decompose pressure signal into three layers with subsignals obtained. Three wavelet packet quantifiers for seven subsignals, including mean absolute value of coefficients, wavelet packet energy and entropy, are compared. The three quantifiers are correlated with maximum pressure rise rate and pressure rise acceleration, respectively. The analysis shows that the in-cylinder pressure, temperature and ROHR change smoothly in normal combustion. When combustion gets into knock, they have a steep rise and a strong fluctuation; the ROHR peaks increase for both cool flame and hot flame, and heat release advances, especially for hot flame. The pressure rise rate and pressure rise acceleration fluctuate more violently, and their maximums increase remarkably and advance somewhat. Without knock, mean absolute value of coefficients, wavelet packet energy and entropy for the subsignal 1 are much greater than others. As knock occurs, three wavelet packet quantifiers for seven subsignals increase greatly, and for the subsignal 6 becomes the largest. Wavelet packet quantifiers for seven subsignals should be monitored for knock detection. The correlation coefficient similarly increases first, decreases afterwards and increases again through seven subsignals. Among three wavelet packet quantifiers, mean absolute value of coefficients has the maximum correlation coefficient except for the subsignal 6. Its maximum correlation coefficient appears at subsignal 7, whose frequency band is 8.75–10 kHz.     

Numerical analysis of the effect of the hydrogen composition on a partially premixed gas turbine combustor

Large-eddy simulations (LESs) of a hydrogen-enriched 1/3-scale GE7EA gas turbine combustor are conducted. Four different fuel compositions are employed to investigate the role of the CH₄/H₂ syngas composition on the resulting flame structure and pressure oscillations occurring inside the combustor. A comparison with the experimental data is conducted to validate the numerical results. First, imaging processing is performed using an Abel-inversion technique for the accumulated OH mass fraction showing good agreement with the experimental images. Then, the calculated velocity fields are successfully compared to the experimental (particle image velocimetry) results. The results show that the flame structure is readily altered when changing the syngas composition; this strongly affects the flow field and therefore the pressure oscillations inside the combustor. When the hydrogen composition is increased, the flame becomes shorter and thicker, and its effect on the outer recirculation zone is minimized. When the flame length approaches the radial length of the combustor under certain conditions, the flame periodically attaches to the rigid wall and the pressure oscillations inside the combustor become amplified. Overall, the LES combined with the multi-step kinetics successfully predicts the variation in the flow fields due to fuel composition changes and reveals the role of the syngas composition in the combustor.     

PaSR湍流燃烧模型对典型湍流射流火焰的数值模拟

采用PaSR湍流燃烧模型对湍流燃烧研究中典型的甲烷湍流射流火焰进行了数值模拟.计算采用简化的化学反应机理,并将计算得到的平均温度场、速度场和各组分的分布与相应的权威实验数据进行了对比.对反应系统中流动和燃烧的不同时间尺度以及二者之间的关系作了探讨.计算结果表明,PaSR模型能够很好地模拟燃烧过程中流场和组分的变化.在火...     

Effects of hydrogen concentration and film thickness on the vented explosion in a small obstructed rectangular container

The explosion venting is an effective way to reduce hydrogen-air explosion hazards, but the explosion venting has been less touched in an obstructed container. The present study mainly focused on the effects of hydrogen concentration and film thickness on the explosion venting in a small obstructed rectangular container. High speed schlieren photography was employed to obtain the flame fine structure and velocity. Pressure transducers were used to measure the overpressure nearby the obstacle. The experimental results show that the obstacle has a significant effect on the flame shape, tip speed and overpressure. In the process of flame evolution, the flame surface becomes more wrinkled with time after the tulip flame. Compared with the cases without the obstacle, the flame surface becomes more distorted and wrinkled downstream of the obstacle under the influence of obstacle enhanced turbulence and flow instability. Upstream of the obstacle, the lower part of the flame surface becomes concave while the upper part shows convex. The pressure histories show that the maximum overpressure increases with the hydrogen concentration in the range of 11.8%–23.7%. Two main pressure peaks were observed for all hydrogen concentrations in the presence of the obstacle. The Helmholtz oscillations appear after the second pressure peak and its duration increases slightly when the hydrogen concentration increases. The combined effect of the obstacle and hydrogen concentration on the second peak overpressure is more significant than on the first peak overpressure. Moreover, the maximum overpressure shows a monotonic increase with the film thickness.     

Experimental study of premixed syngas/air flame deflagration in a closed duct

The propagation behaviour of a deflagration premixed syngas/air flame over a wide range of equivalence ratios is investigated experimentally in a closed rectangular duct using a high-speed camera and pressure transducer. The syngas hydrogen volume fraction, φ, ranges from 0.1 to 0.9. The flame propagation parameters such as flame structure, propagation time, velocity and overpressure are obtained from the experiment. The effects of the equivalence ratio and hydrogen fraction on flame propagation behaviour are examined. The results indicate that the hydrogen fraction in a syngas mixture greatly influences the flame propagation behaviour. When φ, the hydrogen fraction, is ≥0.5, the prominently distorted tulip flame can be formed in all equivalence ratios, and the minimum propagation time can be obtained at an equivalence ratio of 2.0. When φ < 0.5, the tulip flame distortion only occurs in a hydrogen fraction of φ = 0.3 with an equivalence ratio of 1.5 and above. The minimum flame propagation time can be acquired at an equivalence ratio of 1.5. The distortion occurs when the maximum flame propagation velocity is larger than 31.27 m s^?1. The observable oscillation and stepped rise in the overpressure trajectory indicate that the pressure wave plays an important role in the syngas/air deflagration. The initial tulip distortion time and the plane flame formation time share the same tendency in all equivalence ratios, and the time interval between them is nearly constant, 4.03 ms. This parameter is important for exploring the quantitative theory or models of distorted tulip flames.     

Numerical investigation of high pressure and high Reynolds diffusion flame using Large Eddy Simulation

Today,with nonstop improvement in computational power,Large-Eddy Simulation（LES） is a high demanding research tool for predicting engineering flows.Such flows on high pressure condition like diesel engines is extensively employed in ground and marine transportation,oblige the designer to control and predict toxic pollutants,while maintaining or improving their high thermal efficiency.This becomes one of the main challenging issues in decades.In the present work,numerical investigation of diffusion flame dynamics is performed in the near-field of high-Reynolds jet flow on high pressure condition encountered in diesel engine applications.This work discusses the implementation of Partially Stirred Reactor（PaSR） combustion model by the approaches of large eddy simulation（LES）.The simulation results show that LES,in comparison with Reynolds-Averaged Navier-Stokes（RANS） simulation predicts and captures transient phenomena very well.These phenomena such as unsteadiness and curvature are inherent in the near-field of high Reynolds diffusion flame.The outcomes of this research are compared and validated by other researchers＇ results.Detailed comparisons of the statistics show good agreement with the corresponding experiments.     

Direct numerical simulation of nonpremixed flame-wall interactions

The objective of the present study is to use detailed numerical modeling to obtain basic information on the interaction of nonpremixed flames with cold wall surfaces. The questions of turbulent fuel-air-temperature mixing, flame extinction, and wall-surface heat transfer are studied using direct numerical simulation (DNS). The DNS configuration corresponds to an ethylene-air diffusion flame stabilized in the near-wall region of a chemically inert solid surface. Simulations are performed with adiabatic or isothermal wall boundary conditions and with different turbulence intensities. The simulations feature flame extinction events resulting from excessive wall cooling and convective heat transfer rates up to 90 kW/m². The structure of the simulated wall flames is studied in terms of a classical mass-mixing variable, the fuel-air based mixture fraction, and a less familiar heat loss variable, the excess enthalpy variable, introduced to provide a measure of nonadiabatic behavior due to wall cooling. In addition to the flame structure, extinction events are also studied in detail and a modified flame extinction criterion that combines the concepts of mixture fraction and excess enthalpy is proposed and then tested against the DNS data.     

Study of the performance of three micromixing models in transported scalar PDF simulations of a piloted jet diffusion flame (“Delft Flame III”)

Numerical simulation results are presented for a turbulent nonpremixed flame with local extinction and reignition. The transported scalar PDF approach is applied to the turbulence-chemistry interaction. The turbulent flow field is obtained with a nonlinear two-equation turbulence model. A C₁ skeletal scheme is used as the chemistry model. The performance of three micromixing models is compared: the interaction by exchange with the mean model (IEM), the modified Curl's coalescence/dispersion model (CD) and the Euclidean minimum spanning tree model (EMST). With the IEM model, global extinction occurs. With the standard value of model constant C?=2, the CD model yields a lifted flame, unlike the experiments, while with the EMST model the correct flame shape is obtained. However, the conditional variances of the thermochemical quantities are underestimated with the EMST model, due to a lack of local extinction in the simulations. With the CD model, the flame becomes attached when either the value of C? is increased to 3 or the pilot flame thermal power is increased by a factor of 1.5. With increased value of C? better results for mixture fraction variance are obtained with both the CD and the EMST model. Lowering the value of C? leads to better predictions for mean temperature with EMST, but at the cost of stronger overprediction of mixture fraction variance. These trends are explained as a consequence of variance production by macroscopic inhomogeneity and the specific properties of the micromixing models. Local time stepping is applied so that convergence is obtained more quickly. Iteration averaging reduces statistical error so that the limited number of 50 particles per cell is sufficient to obtain accurate results.     

Experimental study of industrial gas turbine flames including quantification of pressure influence on flow field,fuel/air premixing and flame shape

A commercial swirl burner for industrial gas turbine combustors was equipped with an optically accessible combustion chamber and installed in a high-pressure test-rig. Several premixed natural gas/air flames at pressures between 3 and 6 bar and thermal powers of up to 1 MW were studied by using a variety of measurement techniques. These include particle image velocimetry (PIV) for the investigation of the flow field, one-dimensional laser Raman scattering for the determination of the joint probability density functions of major species concentrations, mixture fraction and temperature, planar laser induced fluorescence (PLIF) of OH for the visualization of the flame front, chemiluminescence measurements of OH^* for determining the lift-off height and size of the flame and acoustic recordings. The results give insights into important flame properties like the flow field structure, the premixing quality and the turbulence–flame interaction as well as their dependency on operating parameters like pressure, inflow velocity and equivalence ratio. The 1D Raman measurements yielded information about the gradients and variation of the mixture fraction and the quality of the fuel/air mixing, as well as the reaction progress. The OH PLIF images showed that the flame was located between the inflow of fresh gas and the recirculated combustion products. The flame front structures varied significantly with Reynolds number from wrinkled flame fronts to fragmented and strongly corrugated flame fronts. All results are combined in one database that can be used for the validation of numerical simulations.     

利用离子信号计算已燃质量分数

在定容燃烧弹中分别用离子信号和压力对已燃质量分数进行了研究,通过对点火电极附近离子信号的研究表明,其信号分别在点火、火焰前锋与后区中存在3个峰值。利用从点火到达火焰前锋与后区两个峰值的时间间隔以及所测量的压力值分别计算了已燃质量分数,其结果表明,两方法所计算出的已燃质量分数比较一致,当相对空燃比等于1时,已燃质量分数曲线较为陡峭;而当其偏离1时,曲线逐渐平坦,并且起始位置也愈加滞后,选取系数为0.98的计算结果表明,当相对空燃比在0.9～1.1附近时,两者计算结果的最大差值小于10%。     

Numerical simulation on the spontaneous ignition of high-pressure hydrogen release through a tube at different burst pressures

Spontaneous ignition induced by high-pressure hydrogen release is one of the huge potential risks in the promotion of hydrogen energy. However, the understanding of the microscopic dynamic characteristics of spontaneous ignition, such as ignition initiation and flame development, remains unresolved. In this paper, the spontaneous ignition caused by high-pressure hydrogen release through a tube is investigated by two-dimensional numerical simulation at burst pressure ranging from 2.67 to 15 MPa. Especially, the thermal and species characteristics in hydrogen shock-induced ignition under different strengths of shock wave are discussed carefully. The results show that the stronger shock wave caused by higher burst pressure leads to larger heating area and higher heating temperature inside the tube, increasing the possibility of spontaneous ignition. The shortening effect of initial ignition time and initial ignition distance will decrease with the increase of the burst pressure. Ignition will be initiated when the temperature is raised to about 1350–1400 K under the heating effect of shock waves. It is also found that the ignition occurs under the lean-fuel condition firstly on the upper and lower walls of the tube. The flame branch after spontaneous ignition is observed in the mixing layer. Two ignition kernels show different characteristics during the process of combustion and flow. The evolution of HRR and mass fraction of key species (OH, H, HO₂) are also compared to identify the flame front. The mass fraction of H has the better trend with HRR. It is suggested that H radical is a more reasonable choice as the indicator of the flame front.