基于浓度敏感性分析和遗传算法的甲烷燃烧机理简化与优化

采用良搅拌反应器模型和层流预混火焰模型计算甲烷/空气燃烧过程,通过元素流通法和浓度敏感性分析法,对甲烷燃烧详细化学动力学机理GRIMECH 3.0进行简化。利用遗传算法,以甲烷/空气详细机理获得的组分浓度和一维层流火焰速度为目标,对简化机理进行优化。结果表明,相比于优化之前的简化机理,优化后的简化机理在描述甲烷/空气燃烧反应的组分浓度、层流火焰速度以及反应物和产物的时空分布方面,具有更高的精度。     

Experimental study of the effect of nitrogen addition on gas explosion

Flame propagation and combustion characteristics of methane/air mixed gas in gas explosion were studied in a constant volume combustion bomb. Stretched flame propagation velocity, unstretched laminar flame propagation velocity, unstretched laminar combustion velocity and Markstein length were obtained at various ratios of nitrogen to gas mixture. Combustion stability at various ratios of nitrogen to gas mixture was analyzed by analyzing the pictures of flame propagation. Furthermore, the effect of initial pressure on the flame propagation and combustion characteristics of methane/air mixed gas in gas explosion was analyzed. The results show that the unstretched laminar flame propagation velocity, the unstretched laminar combustion velocity, Markstein length, flame stability, and the maximum combustion pressure decrease distinctly with the increase of nitrogen fraction in the gas mixture. At the same ratios of nitrogen to gas mixture, Markstein length, unstretched laminar flame propagation velocity and unstretched laminar combustion velocity decrease and the maximum combustion pressure increase with the increase of initial pressure of the gas mixture. When nitrogen fraction in the gas mixture is over 20%, the flame will be unstable and is easy to exterminate.     

A constant pressure dust explosion experiment

A new apparatus has been designed for investigating flame propagation in turbulent dust clouds at near constant pressure conditions. The experimental approach is inspired by the classical soap bubble method for measuring burning velocities in gaseous mixtures. Combustible dust is dispersed with pressurised air to form an explosive mixture inside a transparent latex balloon. After a certain delay time, the turbulent dust cloud is ignited by a 40 J chemical igniter. A digital high-speed video camera records the propagating flame and the expansion of the balloon. Experiments were performed with two types of dust, Lycopódium spores and maize starch, as well as with propane–air mixtures under initially quiescent or turbulent conditions. Although the results are primarily qualitative in nature, they nevertheless demonstrate fundamental differences between premixed combustion of gaseous mixtures, and ‘premixed combustion with non-premixed substructures' in mechanical suspensions of solid particles dispersed in air. The discussion highlights some fundamental challenges for future dust explosion research.     

2,5-二甲基呋喃-空气燃烧过程研究

利用球形发展火焰研究了常温常压下不同当量比,不同相态时2,5-二甲基呋喃-空气的层流燃烧速度和马克斯坦长度,分析了火焰拉伸对火焰传播速度的影响。研究结果表明:随着当量比的增加,2,5-二甲基呋喃-空气混合气的马克斯坦长度减少,火焰的稳定性减弱。并且分别计算出当量比为1.25和1.5的层流燃烧速度,分别为:1.189m/s,1.135m/s.。对于同一当量比1.5的情况下,不同相态的2,5-二甲基呋喃-空气混合物,在相同时刻的气液两相混合物的火焰半径已经拉伸火焰传播速度远远大于纯气相的混合物。     

Reaction zone structures and propagation mechanisms of flames in stearic acid particle clouds

Reaction zone structures and propagation mechanisms of two representative flames established in stearic acid (CH₃(CH₂)₁₆CO₂H) particle clouds have been investigated. The reacting zone structure was examined by using a micro-electrostatic probe and a high-speed schlieren system. A distinct difference was observed in the ion current fluctuations recorded across the two representative flames propagating through the clouds of the same total mass density of particles and different mass densities of the particles smaller than 60 μm in diameter. When the mass density of smaller particles was high, a single peak was recorded in the ion current fluctuation. On the other hand, when the mass density of smaller particles was low, multi-peaks of various heights and widths were recorded. In the former case, the single peak was considered to be attributable to a unitary and a relatively thin flame started burning in vapor generated by the evaporation of smaller particles in the preheat zone. The flame propagation mechanism in this case was inferred to be similar to that of a usual hydrocarbon–air premixed flame, although the reaction zone thickness is much larger than that of the premixed flame. In the latter case, the multi peaks of various shapes were considered to be attributable to strong combustion at blue spots far behind the schlieren front. The flame propagation in this case was inferred to be supported by the heat release due to combustion at the blue spots.     

非线性拉伸对甲烷-空气预混火焰燃烧特性的影响

使用定容燃烧弹与高速纹影照相系统研究了不同当量比下甲烷-空气预混气体的层流火焰燃烧特性。实验数据同时应用传统线性模型和非线性模型分析了不同当量比对球形扩展火焰的传播速率和马克斯坦长度的影响。结果显示:随着当量比的增加,层流燃烧速率先增大后减小,直到当量比为1.1时,火焰速率达到最大值。马克斯坦长度始终为正值,且随着当量比的增大而增大。在所有当量比条件下,线性和非线性方法计算的火焰速率大致相同,差值小于0.01 m/s;线性方法得到的马克斯坦长度均大于非线性模型计算的结果,并随着当量比的增大,两种方法得到的马克斯坦长度的差值更加显著。     

Modeling the initial flame acceleration in an obstructed channel using large eddy simulation

The propagation and acceleration of a flame surface past obstructions in a closed square channel was investigated using large eddy simulation. The dynamic Smagorinsky–Lilly subgrid model and the Boger flame surface density combustion model were used. The geometry is essentially two-dimensional with fence-type obstacles distributed on the top and bottom surfaces, equally spaced along the channel length at the channel height. Flame propagation, however, is three dimensional as ignition occurs at a point at the center of the channel cross-section. The effect of obstacle blockage ratio on the development of the flame structure was investigated by varying the obstacle height. Three-dimensional cases were simulated from the initiation of a combustion kernel through spark ignition to the acceleration of the flame front at speeds up to 80 m/s. The transition from laminar flame propagation to turbulent flame propagation within the “thin reaction zone” regime was observed in the simulations. By analyzing the development of the three dimensional flame surface and unburned gas flow field, the formation of several flame structures observed experimentally are explained. Global quantities such as the total flame area and centerline flame velocity were ascertained and compared to the experimental data. High amplitude oscillations in the centerline flame velocity were found to occur from a combination of the unburned gas flow field and fluctuations in the volumetric burning rate.     

丙烷/空气拉伸火焰传播稳定性实验研究

为了揭示空气中丙烷火焰传播特性,利用纹影系统记录了预混气体小能量点火条件下火焰形成与传播过程,得到了火焰表面的微观结构特征,分析了混合气体火焰的稳定性及其影响因素。结果表明:丙烷/空气混合物火焰发展过程及其表面微观特征与浓度直接相关;当混合物浓度接近爆炸上下限时,火焰扩展速率整体不大于0.5 m/s,燃烧区域向上漂浮,浮力成为影响火焰失稳的主导因素;当混合物浓度靠理论配比时,火焰呈规则球形扩展,火焰稳定性按照先减弱后增强的趋势发展,火焰表面褶皱的形成及演化规律是热扩散不稳定性和流体力学不稳定性共存与竞争的作用结果。     

管道内甲烷/空气预混火焰加速传播特性研究

搭建了基于可视化燃烧管道的预混火焰精细结构实验台,通过理论分析和实验等手段对预混火焰的瞬态传播过程及加速传播特性进行了研究,分析得出了不同甲烷含量的预混火焰加速传播规律及预混火焰由层流向湍流的转变规律.     

Behavior of flames propagating through lycopodium dust clouds in a vertical duct

The structure of flame propagating through lycopodium dust clouds has been investigated experimentally. Upward propagating laminar flames in a vertical duct of 1800 mm height and 150×150 mm square cross-section are observed, and the leading flame front is also visualized using by a high-speed video camera. Although the dust concentration decreases slightly along the height of duct, the leading flame edge propagates upwards at a constant velocity. The maximum upward propagating velocity is 0.50 m/s at a dust concentration of 170 g/m³. Behind the upward propagating flame, some downward propagating flames are also observed. Despite the employment of nearly equal sized particles and its good dispersability and flowability, the reaction zone in lycopodium particles cloud shows the double flame structure in which isolated individual burning particles (0.5–1.0 mm in diameter) and the ball-shaped flames (2–4 mm in diameter; the combustion time of 4–6 ms) surrounding several particles are included. The ball-shaped flame appears as a faint flame in which several luminous spots are distributed, and then it turns into a luminous flame before disappearance. In order to distinguish these ball-shaped flames from others with some exceptions for merged flames, they are defined as independent flames in this study. The flame thickness in a lycopodium dust flame is observed to be 20 mm, about several orders of magnitude higher than that of a premixed gaseous flame. From the microscopic visualization, it was found that the flame front propagating through lycopodium particles is discontinuous and not smooth.     

On the determination of the laminar burning velocity from closed vessel gas explosions

A methodology to determine the laminar burning velocity from closed vessel gas explosions is explored. Unlike other methods which have been used to measure burning velocities from closed vessel explosions, this approach belongs to the category which does not involve observation of a rapidly moving flame front. Only the pressure–time curve is required as experimental input. To verify the methodology, initially quiescent methane–air mixtures were ignited in a 20-l explosion sphere and the equivalence ratio was varied from 0.67 to 1.36. The behavior of the pressure in the vessel was measured as a function of time and two integral balance models, namely, the thin-flame and the three-zone model, were fitted to determine the laminar burning velocity. Data on the laminar burning velocity as a function of equivalence ratio, pressure and temperature, measured by a variety of other methods have been collected from the literature to enable a comparison. Empirical correlations for the effect of pressure and temperature on the laminar burning velocity have been reviewed and two were selected to be used in conjunction with the thin-flame model. For the three-zone model, a set of coupled correlations has been derived to describe the effect of pressure and temperature on the laminar burning velocity and the laminar flame thickness. Our laminar burning velocities are seen to fall within the band of data from the period 1953–2003. A comparison with recent data from the period 1994–2003 shows that our results are 5–10% higher than the laminar burning velocities which are currently believed to be the correct ones for methane–air mixtures. Based on this observation it is concluded that the methodology described in this work should only be used under circumstances where more accurate methods can not be applied.     

Estimation of turbulent flame speed during DME/air premixed gaseous explosions

DME is thought to be a good alternative fuel due to its cleanliness and more excellent fuel economy. Although the prediction and loss prevention of flammability hazard is very important for safety of DME installations, the evaluation method with sufficient accuracy has not been established. In this study, a numerical combustion model is constructed and a 3-dimensional computational fluid dynamics (CFD) simulation of a premixed DME/air explosion in a large-scale domain is conducted. The main feature of the numerical model is the solution of a transport equation for the reaction progress variable using a function for turbulent flame velocity which characterizes the turbulent regime of propagation of free flames derived by introducing the fractal theory. The model enables the calculation of premixed gaseous explosion without using fine mesh of the order of micrometer, which would be necessary to resolve the details of all instability mechanisms. The value of the empirical constant contained in the function for turbulent flame velocity is evaluated by analyzing the experimental data of LPG/air and DME/air premixed explosions. The comparison of flame behavior between the experimental result and numerical simulation shows good agreement.     

Structure and flame speed of dilute and dense layered coal-dust explosions

Multidimensional time-dependent simulations were performed to study the interaction of a stoichiometric methane–air detonation with layers of coal dust. The simulations solved equations representing a Eulerian kinetic-theory-based granular multiphase model applicable to dense and dilute particle volume fractions. These equations were solved using a high-order Godunov-based method for compressible fluid dynamics. Two dust layer concentrations were considered: loose with an initial volume fraction of 1%, and dense with an initial volume fraction of 47%. Each layer was simulated with two types of dust: reactive coal and inert ash. Burning of the coal particles results in a coupled complex consisting of an accelerating shock leading a coal-dust flame. The overall structure of the shock–flame complex resembles that of a premixed fast flame with length scales on the order of several meters. The large length scales are direct results of time needed to lift, mix, heat, and autoignite the particle. The flame speeds are large and much larger than the gas-phase velocity. Large spikes of flame speed are characteristic of the 47% case. These spikes and high flame speed are caused by pockets of coal dust autoigniting ahead of the flame. The flame is choked in the 1% case due to the gas-phase products exceeding the sonic velocity with respect to the flame. The 47% case is choked due to attenuation of pressure waves as they propagate through particles. Inert layers of dust substantially reduce the overpressure, impulse, and speed produced by propagating blast wave. The results also show that loose layers of dust are far more dangerous than dense layers. The shock and flame are more strongly coupled for loose layers, propagate at higher velocity, and produce large overpressures and impulses.     

基于重要性分析的甲烷/乙烯混合燃烧框架机理构筑

建立了基于重要性分析的详细化学机理分析平台,利用耦合组分化学存活时间和敏感性系数的重要性参数,确定详细机理中的准稳态组分,通过移除这些组分及其相关反应,得到了计算精度较高的框架机理。针对目前燃烧学界较为关注的混合燃烧问题,以甲烷、乙烯这两种典型低碳碳氢燃料为研究对象,对其详细化学反应机理进行了分析,利用重要性分析法构筑框架机理,并对甲烷/空气和甲烷/乙烯/空气预混火焰进行了数值计算。与详细机理相比,框架机理所涉及的组分数与基元反应数都得到了大幅度的降低,计算时间明显减少,但对火焰温度及反应物、生成物、中间组分浓度的预测与采用详细机理得到的结果吻合良好,证明了重要性分析法的有效性与可靠性。     

3m直径煤油池火灾火焰特性的数值研究

为了预测油池火灾的火焰特性,采用CFD模拟技术开展静风状态下3 m直径煤油液池的火灾场景模拟,探讨火焰温度、火焰羽流速度、辐射热通量、燃烧产物质量分数等油池火焰特性参数随高度的变化关系;并结合火焰形态分布,提出一种4区域模型,即将湍流扩散火焰划分为油气混合燃烧区、燃烧火焰区、烟尘区和热烟气区来分析燃烧气流在不同高度的实际物理化学特性。此外,通过经验公式和CFD模拟2种方法分别计算出3 m直径煤油池火灾的火焰高度、火焰表面的辐射通量及热辐射破坏半径,并对计算结果进行比较分析,结果表明:2种方法可互相补充完善,有助于池火灾的热辐射危害性评估。     

Numerical simulation of premixed methane–air deflagration in a semi-confined obstructed chamber

In this paper, large eddy simulation coupled with a turbulent flame speed cloure (TFC) subgrid combustion model has been utilized to simulate premixed methane–air deflagration in a semi-confined chamber with three obstacles mounted inside.The computational results are in good agreement with published experimental data, including flame structures, pressure time history and flame speed. The attention is focused on the flame flow field interaction, pressure dynamics, as well as the mechanism of obstacle-induced deflagration. It is found that there is a positive feedback mechanism established between the flame propagation and the flow field. The pressure time history can be divided into four stages and the pseudo-combustion concept is proposed to explain the pressure oscillation phenomenon. The obstacle-induction mechanism includes direct effect and indirect effect, but do not always occur at the same time.     

Effects of particle characteristics on flame propagation behavior during organic dust explosions in a half-closed chamber

To reveal the effects of particle characteristics on the mechanisms of flame propagation during organic dust explosions clearly, three long chain monobasic alcohols which are solids at room temperature and have similar physical–chemical properties were chosen to carry out experiments in a half-closed small chamber. A high-speed video camera was used to record the flame propagation process and to obtain the direct light emission photographs. Flame temperature was detected by a fine thermocouple. Based on the experimental results above, analysis was conducted on flame propagation characteristics and temperature profiles of organic particle cloud. As a result, it was found that the particle materials, especially volatility, strongly affected the flame propagation behavior. Particle concentration also affects the combustion zone propagation process significantly. With increasing the particle concentration, the maximum temperature of the combustion zone increases at the lower concentration, reaches a maximum value, and then decreases at the higher concentration. The propagation velocity of the combustion zone has a linear relationship with the maximum temperature, which implies conductive heat transfer is dominant in the flame propagation process of the three different volatile dusts.     

Simulation of dust explosions in complex geometries with experimental input from standardized tests

The present paper describes the development of a new CFD-code (DESC) for the assessment of accidental hazards arising from dust explosions in complex geometries. The approach followed entails the estimation of the laminar burning velocity of dust clouds from standardized laboratory-scale tests, and its subsequent use as input to the combustion model incorporated in DESC. The methodology used to obtain the laminar burning velocities is demonstrated by igniting turbulent propane-air mixtures to deflagration in a standard 20-litre USBM-vessel, and extracting the laminar burning velocity from the pressure–time curves; the results are compared with literature data. Laminar burning velocities for clouds of maize starch dust in air were estimated following the same procedure, and the resulting empirical model was used to simulate dust explosions in a 236-m³ silo.     

多分支管道油气爆炸特性大涡模拟

为了研究油库常见的分支结构空间内发生油气爆炸时火焰和压力的传播特性,建立了基于WALE湍流模型及Zimont预混火焰模型的油气爆炸模型;模拟了6种不同分支管道结构空间内汽油/空气混合物爆炸发生发展过程;研究了分支管道数量及相对设置位置对爆炸超压的影响规律,以及分支管道对火焰传播形态和速度的影响规律;模拟结果与前人相关实验规律进行对比。研究结果表明:分支管道对汽油/空气混合气预混爆炸具有明显的强化激励作用;火焰锋面传播经过分支管道时,经历规则—褶皱—规则的变化过程;主管道内火焰传播速度,在分支管道对流场的突扩作用和湍流作用的共同影响下呈震荡变化的规律。