障碍物对预混火焰特性影响的大涡数值模拟

障碍物在预混气体火焰传播过程中对其结构及传播特性造成较大影响,对火焰的加速和爆燃转爆轰过程(deflagration-to-detonation transition, DDT)起到直接的促进作用。通过障碍物条件下可视管道中甲烷/空气预混火焰传播实验,捕获其火焰微观结构变化。采用三维物理模型,采用壁面自适应局部涡黏模型(wall-adapting local eddy-viscosity, WALE)的大涡模拟(large eddy simulation, LES),并用火焰增厚化学反应模型(thickened flame model, TFM)对实验过程进行重现。分析开口管道中预混火焰翻越障碍物后的复杂流场变化,并分析层流向湍流转变过程的特点。揭示了在障碍物影响下预混火焰扰动失稳现象的直接原因,是由障碍物引发的3个气流涡团同时作用而形成Kelvin-Helmholtz不稳定及Rayleigh-Taylor不稳定现象耦合作用所导致。     

含双侧分支受限空间油气爆炸火焰行为与超压特性大涡模拟

为研究含分支结构狭长受限空间油气爆炸特性规律，基于大涡模拟WALE模型和Zimont预混火焰模型，对横截面为100 mm×100 mm的含双侧分支管道受限空间油气泄压爆炸特性进行了数值模拟。通过对火焰形态、火焰传播速度和动态超压3个物理量的对比，验证了所建立模型对于含分支结构受限空间油气爆炸计算的适用性。基于数值模拟结果，对爆炸过程中的流场结构、火焰形态和超压变化规律进行了分析，指出了“浪花状”火焰的形成原因。结果表明:（1）火焰传播进入分支管道前，在主管道和分支管道交界处会产生旋转方向相反的对称涡旋结构，并随着火焰传播不断向分支管道内部发展；（2）当火焰传播进入分支管道后，分支管道内部前期已建立流场决定了火焰的形态，火焰锋面在涡旋结构作用下呈“浪花状”，此后火焰和流场相互影响，流场向湍流转捩，火焰锋面褶皱变形；（3）爆炸超压升压过程可划分为4个阶段，受到火焰锋面面积和分支管道泄压共同作用，表明爆炸流场、火焰行为和动态超压呈现出显著耦合性。     

方管内汽油-空气混合气体密闭爆炸和泄爆特性研究

为研究汽油-空气混合气体密闭爆炸和泄爆特性，采用可视化方管进行了两种爆炸模式实验研究，并基于壁面自适应局部涡黏(wall-adapting local eddy-viscosity,WALE)模型和Zimont预混火焰模型进行了数值模拟研究。结果表明:（1）泄爆工况超压-时序曲线峰值数量多于密闭爆炸工况，且泄爆工况超压-时序曲线存在剧烈的类似简谐振动的振荡，而密闭爆炸工况的爆炸超压特征参数显著高于泄爆工况；（2）密闭爆炸工况最大火焰传播速度明显小于泄爆工况，但前者在火焰传播初期即达到最大值，而后者在火焰传播末期才达到最大值；（3）密闭爆炸工况出现郁金香形火焰，而泄爆工况出现蘑菇形火焰，郁金香火焰的形成与管道内火焰锋面、流场和流场动压三者之间耦合效应相关，蘑菇形火焰由外部流场湍流和斜压效应的共同作用引起。     

多孔障碍物对预混火焰传播的影响

以甲烷为代表性气体，研究了半封闭管道中设置多孔障碍物对可燃气体爆炸火焰传播的影响，基于大涡模拟对实验进行了重现，对比了实验与模拟中火焰传播过程的形状、位置及速度，分析了模拟结果中火焰穿过障碍物前后的流场和表面积变化，给出了衡量火焰褶皱程度的指标及算法。结果表明:大涡模拟结果与实验结果有较好的一致性；火焰在存在障碍物的管道内传播，经历层流快速膨胀、受阻回流、湍流快速发展和脉动减速4个阶段，各阶段火焰依次分别呈现加速、减速、二次加速、二次减速的波动变化；当可燃气体在开口与点火位置同端的管道内爆炸，火焰在接近障碍物时，受管道封闭端和障碍物约束显著，而出现脉动回流现象；火焰穿过多孔障碍物后，传播速度骤升至峰值，较未穿过障碍物前的最大速度可增加58.7%；障碍物是导致火焰面破碎以及面积褶皱率增大的直接原因，火焰褶皱率最大可达44.8%，比未穿过障碍物前的最大褶皱率增大39.27%。     

连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟

采用WALE模型和Zimont预混火焰模型对内置圆孔障碍物油气泄压爆炸火焰特性进行了大涡模拟,并将大涡模拟计算结果和RNG k-ε湍流模型计算结果以及实验结果进行对比分析,验证了大涡模拟的精确性。结果表明:（1）大涡模拟在预测油气爆炸超压、火焰传播速度以及火焰形态变化等方面比RNG k-ε湍流模型精确度更高,且能表现出更多流场的精细化结构;（2）障碍物诱导管道内形成湍流度较高的流场区域,导致火焰产生褶皱弯曲变形,增大火焰面积,加速火焰传播;（3）爆炸超压、火焰传播速度和火焰面积内在联系密切,具有显著的耦合性,且随时间的变化趋势存在高度的一致性。     

甲烷/空气预混气体火焰的传播特征

利用高速纹影摄像等技术探讨了密闭管道内不同当量比的甲烷/空气预混气体火焰的传播特征。结果表明,当甲烷含量接近当量值时,预混气体火焰传播中会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程,逐渐由层流燃烧转变成湍流燃烧,并形成Tulip火焰结构;当甲烷含量偏离当量值一定程度时,预混火焰呈现出典型的层流燃烧特征,不会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程。Tulip火焰结构形成于火焰传播速度迅速降低的区间里,且只有当减速阶段的最大加速度的绝对值大于某一数值时才能形成;Tulip火焰结构是预混火焰由层流燃烧向湍流燃烧转变的一个中间过程。     

燃烧室中悬跨圆柱的火焰绕流

利用高速阴影照相系统,对方形管道内火焰流经悬吊圆柱时所产生的变化及其诱导的流场进行了实验研究,并利用高精度PPM格式对上述现象进行了数值模拟,其中,湍流采用大涡模拟(LES)方法,燃烧采用旋涡破碎(EBU)模型,圆柱边界采用沉浸边界法(IBM).实验结果与计算结果比较吻合,在一定程度上揭示了火焰阵面变形和加速的规律. 根据研究结果,对火焰与悬吊圆柱相互作用过程中火焰的三维形状、流场中涡量和湍流强度的分布、火焰与流场的相互影响以及火焰传播轨迹等进行了讨论.      

湍流强度对预混燃烧影响的大涡模拟研究

本文对甲烷和空气在突扩燃烧室内的预混燃烧进行了大涡模拟(LES)研究,主要关注于预混燃料的湍流强度对甲烷-空气预混燃烧效率的影响.LES数值模拟给出了不同入口湍流强度的条件下燃烧室内湍流预混燃烧反应流场的速度场,温度场和浓度场分布.分析的结果表明:随着预混燃料入口湍流强度的增大,燃烧反应更加剧烈,燃烧效率提高,火焰长度也随之变短,同时,一氧化碳的生成区域明显变小,整体生成量减少.研究的结果对有效地提高燃气轮机的燃烧效率和减少污染物的排放提供一定的参考依据.     

波纹管道阻火器内火焰传播的实验与数值模拟研究

对乙烯-空气预混火焰在波纹管道阻火器中的传播与淬熄过程进行了实验和数值模拟研究,实验结果显示:当乙烯接近当量浓度时,预混气体爆炸压力变化过程可分为4个阶段,等压燃烧阶段、缓慢上升阶段、快速上升阶段和压力振荡阶段;在爆炸过程中,由于反射压力波和火焰相互作用的影响,超压值出现多次振荡,压力振荡阶段一般可以持续数十毫秒;乙烯-空气火焰传播速度随管径增加、阻火单元波纹高度减小呈递增趋势,而且随着阻火单元厚度的增加,阻火器的阻火能力明显提高,可以更有效地使火焰淬熄。数值模拟结果显示:在管道封闭端点火后,火焰面呈半球形并以层流扩散的方式向四周传播;当火焰传播到管道壁面时,在管道壁面的约束作用下,火焰面发生变形,壁面附近的火焰逐渐超过了管道轴线附近的火焰,最后形成了“郁金香”状的火焰结构;当爆燃火焰经过阻火单元时,高温已燃气体被其吸收大量热量,同时在反应区产生的稀疏波作用下,气体温度逐渐降低、化学反应速率迅速减小,最终导致火焰被熄灭。通过模拟计算结果可以看出,在整个爆炸过程中,火焰传播速度与爆炸压力波动均较为明显。并提出了孔隙率和阻火单元厚度对火焰传播的影响机制。基于传热学理论模型,并结合实验数据,得出了爆燃火焰速度与爆炸压力之间的关系,为工业装置阻火器的设计和选型提供更为准确的参考依据。     

无约束条件下甲烷/空气预混气体燃爆特性研究

为研究无约束条件下甲烷(CH₄)/空气(air)预混气体的燃爆特性,以乳胶气球为反应容器开展了甲烷爆炸实验,结合Chemkin模拟和改进的比色测温技术,研究了不同当量比下甲烷/空气预混气体的火焰传播速度、爆炸超压及温度场分布等特征以及静置时间对预混气体燃爆特性的影响。实验结果表明：甲烷/空气预混气体的爆炸火焰传播速度呈振荡分布,当量比为0.83、1.06、1.30和1.55时的平均火焰传播速度分别为1.554m·s^-1、2.122m·s^-1、1.892m·s^-1和1.428m·s^-1;峰值超压随当量比的增加呈先增大后减小的趋势,当量比为1.06时基元反应CH₃·+O₂?O·+CH₃O·的敏感性系数最大,从而加速了生成二氧化碳(CO₂)的链反应,使得燃烧化学反应最彻底,峰值超压值最大;静置时间对火焰传播速度和爆炸峰值压力影响显著,最佳静置时间为6min;随着当量比的增大,爆炸火焰的平均温度呈现...     

Large Eddy Simulation and PIV Measurements of Unsteady Premixed Flames Accelerated by Obstacles

In gas explosions, the unsteady coupling of the propagating flame and the flow field induced by the presence of blockages along the flame path produces vortices of different scales ahead of the flame front. The resulting flame–vortex interaction intensifies the rate of flame propagation and the pressure rise. In this paper, a joint numerical and experimental study of unsteady premixed flame propagation around three sequential obstacles in a small-scale vented explosion chamber is presented. The modeling work is carried out utilizing large eddy simulation (LES). In the experimental work, previous results (Patel et al., Proc Combust Inst 29:1849–1854, 2002) are extended to include simultaneous flame and particle image velocimetry (PIV) measurements of the flow field within the wake of each obstacle. Comparisons between LES predictions and experimental data show a satisfactory agreement in terms of shape of the propagating flame, flame arrival times, spatial profile of the flame speed, pressure time history, and velocity vector fields. Computations through the validated model are also performed to evaluate the effects of both large-scale and sub-grid scale (SGS) vortices on the flame propagation. The results obtained demonstrate that the large vortical structures dictate the evolution of the flame in qualitative terms (shape and structure of the flame, succession of the combustion regimes along the path, acceleration-deceleration step around each obstacle, and pressure time trend). Conversely, the SGS vortices do not affect the qualitative trends. However, it is essential to model their effects on the combustion rate to achieve quantitative predictions for the flame speed and the pressure peak.     

方形管内楔形障碍物对火焰结构与传播的影响

通过实验与数值模拟方法对CH4/空气预混火焰在有楔形障碍物的卧式燃烧方管内的传播进行了研究。采用多镜头Cranz Schardin高速摄像机和压力传感器等实验设备获得了高清晰度的障碍物诱导火焰失稳的分幅时序照片以及障碍物背风表面压力变化曲线。数值模拟则基于RANS方法与EDU-Arrhenius燃烧模型,计算结果与实验结果基本相符,反映了火焰在管内传播与变形的详细过程。通过综合分析实验与计算结果,得到了由楔形障碍物导致的火焰加速与变形的内在机理,揭示了火焰传播过程中由层流燃烧向湍流燃烧转捩的本质。     

Experimental investigation of flame/solid interactions in turbulent premixed combustion

An experimental study has been carried out to investigate the interaction between propagating turbulent premixed flames and solid obstacles. The experimental rig was configured specifically to allow detailed measurements with laser-based optical diagnostics. A wall-type solid obstacle was mounted inside a laboratory-scale combustion chamber with rectangular cross-section. The flame was initiated, by igniting a combustible mixture of methane in air at the center of the closed end of the combustion chamber. The flame front development was visualized by a high-speed (9000 frame/s) digital video camera and flame images were synchronized with ignition timing and chamber pressure data. The tests were carried out with lean, stoichiometric and rich mixtures of methane in air. The images were used to calculate highly resolved temporal and spatial data for the changes in flame shape, speed, and the length of the flame front. The results are discussed in terms of the influence of mixture equivalence ratio on the flame structure and resulting overpressure. The reported data revealed significant changes in flame structure as a result of the interaction between the propagating flame front and the transient recirculating flow formed behind the solid obstacle. Combustion images show that the flame accelerates and decelerates as it impinges on the obstacle wall boundaries. It is also found that the mixture concentrations have a significant influence on the nature of the flame/solid interactions and the resulting overpressure. The highest flame speed of 40 m/s was obtained with the unity fuel–air equivalence ratio. Burning of trapped mixture behind the solid obstruction was found to be highly correlated with the flame front length and the rate of pressure rise.     

Modelling of a Premixed Swirl-stabilized Flame Using a Turbulent Flame Speed Closure Model in LES

This paper proposes a combustion model based on a turbulent flame speed closure (TFC) technique for large eddy simulation (LES) of premixed flames. The model was originally developed for the RANS (Reynolds Averaged Navier Stokes equations) approach and was extended here to LES. The turbulent quantities needed for calculation of the turbulent flame speed are obtained at the sub grid level. This model was at first experienced via an test case and then applied to a typical industrial combustor with a swirl stabilized flame. The paper shows that the model is easy to apply and that the results are promising. Even typical frequencies of arising combustion instabilities can be captured. But, the use of compressible LES may also lead to unphysical pressure waves which have their origin in the numerical treatment of the boundary conditions.     

Resolution Requirements in Stochastic Field Simulation of Turbulent Premixed Flames

The spatial resolution requirements of the Stochastic Fields probability density function approach are investigated in the context of turbulent premixed combustion simulation. The Stochastic Fields approach is an attractive way to implement a transported Probability Density Function modelling framework into Large Eddy Simulations of turbulent combustion. In premixed combustion LES, the numerical grid should resolve flame-like structures that arise from solution of the Stochastic Fields equation. Through analysis of Stochastic Fields simulations of a freely-propagating planar turbulent premixed flame, it is shown that the flame-like structures in the Stochastic Fields simulations can be orders of magnitude narrower than the LES filter length scale. The under-resolution is worst for low Karlovitz number combustion, where the thickness of the Stochastic Fields flame structures is on the order of the laminar flame thickness. The effect of resolution on LES predictions is then assessed by performing LES of a laboratory Bunsen flame and comparing the effect of refining the grid spacing and filter length scale independently. The usual practice of setting the LES filter length scale equal to grid spacing leads to severe under-resolution and numerical thickening of the flame, and to substantial error in the turbulent flame speed. The numerical resolution required for accurate solution of the Stochastic Fields equations is prohibitive for many practical applications involving high-pressure premixed combustion. This motivates development of a Thickened Stochastic Fields approach (Picciani et al. Flow Turbul. Combust. X, YYY (2018) in order to ensure the numerical accuracy of Stochastic Fields simulations.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.