燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性;热通量辐射模型估算辐射换热,分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解,液相采用Lagrange处理,并用PSIC算法对其进行求解．通过实验结果和计算结果的比较,表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟,3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程,而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

Investigation of Two-Fluid Methods for Large Eddy Simulation of Spray Combustion in Gas Turbines

An extension of the large eddy simulation (LES) technique to two-phase reacting flows, required to capture and predict the behavior of industrial burners, is presented. While most efforts reported in the literature to construct LES solvers for two-phase flow focus on Euler–Lagrange formulation, the present work explores a different solution (‘two-fluid’ approach) where an Eulerian formulation is used for the liquid phase and coupled with the LES solver of the gas phase. The equations used for each phase and the coupling terms are presented before describing validation in two simple cases which gather some of the specificities of real combustion chamber: (1) a one-dimensional laminar JP10/air flame and (2) a non-reacting swirled flow where solid particles disperse (Sommerfeld and Qiu, Int. J. Multiphase Flow 19(6):1093–1127, 1993). After these validations, the LES tool is applied to a realistic aircraft combustion chamber to study both a steady flame regime and an ignition sequence by a spark. Results bring new insights into the physics of these complex flames and demonstrate the capabilities of two-fluid LES.     

AN ANALYTICAL DESIGN METHODOLOGY FOR AN ANNULAR COMBUSTOR

SUMMARY

This paper describes a computational procedure for the optimization of the performance parameters of a simulated annular combustor. This method has been applied to analyze the influence of the performance parameters and geometries on the annular combustor characteristics and provide a good understanding of combustor internal flow fields, and therefore it can be used for guiding the combustor design process. The approach is based on the solution of governing nonlinear, elliptic partial differential equations for 3-D axisymmetric recirculating turbulent reacting swirling flows and the modelling of turbulence, combustion, thermal radiation and pollutant formation. The turbulence effects are introduced via the modified two-equation κ-ε model. Turbulent combustion is modelled using the κ-ε-g model and a two-step turbulent combustion model is employed for the excess emission of carbon monoxide CO. For the evaluation of the NO pollutant formation rate, the NO pollutant formation model, which takes into account the influence of turbulence, presented here. The radiative heat transfer is handled by the heat flux model. The predictions of the combustor character-istics and performance parameters are made using the present approach.

Predictions of velocity, length of the recirculation zone, combustion efficiency and wall temperature are compared with measurements. Agreement between the predictions and experimental data is very satisfactory.     

大速差射流燃烧室中烟煤与贫煤燃烧的数值模拟

本文基于颗粒相的轨道模型,对大速差射流燃烧室中烟煤粉与贫煤粉的二维流动,混合及燃烧进行了数值模拟,模拟结果从两相耦合的角度,阐明了煤粉颗粒在燃烧室中运动的规律,煤粉与大速差射流诱导的中心气体逆流之间的混合及其对煤粉火焰稳定的影响,指出此种燃烧室中煤粉火焰稳定的回流区燃烧机理,气相流场及回流区的预报结果与实验符合良好。     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 2. Flows with heat transfer and combustion

The objective of this part of the paper is to summarize the information concerning the authors' works in the field of simulation of two-phase gas-particle turbulent flows with heat transfer and combustion. A kinetic equation had been derived for the probability density function (PDF) of the particle velocity, temperature, and mass distributions in turbulent flows. This PDF equation is used for the construction of the governing conservation equations of mass, momentum, and heat transfer in the dispersed particle phase.The numerical scheme incorporates two-phase fluid dynamics, convective and radiative heat transfer, and combustion. The proposed models have been applied to the calculation of various particle-laden turbulent flows in jets, combustion and gasification chambers, and furnaces.     

Pressure-based method for combustion instability analysis

An improved pressure-based method has been applied to predict the two-dimensional instability analysis of liquid-fuelled rocket engines. This method is non-iterative for transient flow calculations and applicable to all-speed flows. Validation cases include the shock-tube problem, the blast flow field and unsteady spraycombusting flows. Computations for the combustion instability analysis were carried out for various combustion parameters such as spray initial conditions and combustor geometries. Unsteady behaviours of the stable and unstable spray flame fields and effects of acoustic oscillations on the fuel droplet vaporization and combustion process are studied in detail. The present numerical model successfully demonstrates the capability of predicting combustion instability as well as fast transient compressible flows at all speeds.     

On the surface equations in two-phase flows and reacting single-phase flows

This paper summarizes the mathematical surface equations which are useful in two-phase flows and single-phase reacting flows. The connection between the interfacial area concentration transport equation for two-phase flows and the flame surface density transport equation for turbulent reacting flows is established. Several analytical examples are given to clarify the physical significance of the different quantities involved in the different transport equations. An introduction to the mathematical treatment of anisotropic interfaces is also given. This theory is illustrated on two different numerical examples: a single inclusion in a simple shear and a single inclusion in an uni-axial elongation.     

Experimental Characterization of Premixed Flame Instabilities of a Model Gas Turbine Burner

In recent years, the NO _x emissions of heavy duty gas turbine burners have been significantly reduced by introducing premixed combustion. These highly premixed burners are known to be prone to combustion oscillations. In this paper, investigations of a single model gas turbine burner are reported focusing on thermo-acoustic instabilities and their interaction with the periodic fluctuations of the velocity and pressure. Phase-locked optical measurement techniques such as LDA and LIF gave insight into the mechanisms.Detailed investigations of a gas turbine combustor rig revealed that the combustor as well as the air plenum oscillate in Helmholtz modes. These instabilities could be attributed to the phase lag of the pressure oscillations between the air plenum and the combustor, which causes an acceleration and deceleration of the air flow through the burner and, therefore, alternating patterns of fuel rich and lean bubbles. When these bubbles reach the reaction zone, density fluctuations are generated which in turn lead to velocity fluctuations and, hence, keep up the pressure oscillations.With increasing the equivalence ratio strong combustion oscillations could be identified at the same frequency. Similarly as with weak oscillations, Helmholtz mode pressure fluctuations are present but the resulting velocity fluctuations in the combustor can be described as a pumping motion of the flow. By the velocity fluctuations the swirl stabilization of the flame is disturbed. At the same time, the oscillating pressure inside the combustor reaches its minimum value. Shortly after the flame expands again, the pressure increases inside the combustor. This phenomenon which is triggered by the pressure oscillations inside the air plenum seems to be the basic mechanism of the flame instability and leads to a significant increase of the pressure amplitudes.     

Efficient simulation of particle-laden turbulent flows with high mass loadings using LES

The paper is concerned with the simulation of particle-laden two-phase flows based on the Euler–Lagrange approach. The methodology developed is driven by two major requirements: (i) the necessity to tackle complex turbulent flows by eddy-resolving schemes such as large-eddy simulation; (ii) the demand to predict dispersed multiphase flows at high mass loadings. First, a highly efficient particle tracking algorithm was developed working on curvilinear, block-structured grids. Second, to allow the prediction of dense two-phase flows, the fluid–particle interaction (two-way coupling) as well as particle–particle collisions (four-way coupling) had to be taken into account. For the latter instead of a stochastic collision model, in the present study a deterministic collision model is considered. Nevertheless, the computational burden is minor owing to the concept of virtual cells, where only adjacent particles are taken into account in the search for potential collision partners. The methodology is applied to different test cases (plane channel flow, combustion chamber flow). The computational results are compared with experimental measurements and good agreement is found.     

Large Eddy Simulation of Reactive Two-Phase Flow in an Aeronautical Multipoint Burner

Because of compressibility criteria, fuel used in aeronautical combustors is liquid. Their numerical simulation therefore requires the modeling of two-phase flames, involving key phenomena such as injection, atomization, polydispersion, drag, evaporation and turbulent combustion. In the present work, particular modeling efforts have been made on spray injection and evaporation, and their coupling to turbulent combustion models in the Large Eddy Simulation (LES) approach. The model developed for fuel injection is validated against measurements in a non-evaporating spray in a quiescent atmosphere, while the evaporation model accuracy is discussed from results obtained in the case of evaporating isolated droplets. These models are finally used in reacting LES of a multipoint burner in take-off conditions, showing the complex two-phase flame structure.     

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.     

Large-eddy Simulation of Triangular-stabilized Lean Premixed Turbulent Flames: Quality and Error Assessment

In this numerical study, an algebraic flame surface wrinkling (AFSW) reaction submodel based on the progress variable approach is implemented in the large-eddy simulation (LES) context and validated against the triangular stabilized bluff body flame configuration measurements i.e. in VOLVO test rig. The quantitative predictability of the AFSW model is analyzed in comparison with another well validated turbulent flame speed closure (TFC) combustion model in order to help assess the behaviour of the present model and to further help improve the understanding of the flow and flame dynamics. Characterization of non-reacting (or cold) and reacting flows are performed using various subgrid scale models for consistent grid size variation with 300,000 (coarse), 1.2 million (intermediate) and 2.4 million (fine) grid cells. For non-reacting flows at inlet velocity of 17?m/s and inlet temperature 288?K, coarse grid leads to over prediction of turbulence quantities due to low dissipation at the early stage of flow development behind the bluff body that convects downstream eventually polluting the resulting solution. The simulated results with the intermediate (and fine) grid for mean flow and turbulence quantities, and the vortex shedding frequency (f_s) closely match experimental data. For combusting flows for lean propane/air mixtures at 35?m/s and 600?K, the vortex shedding frequency increase threefold compared with cold scenario. The predicted results of mean, rms velocities and reaction progress variable are generally in good agreement with experimental data. For the coarse grid the combustion predictions show a shorter recirculation region due to higher turbulent burning rate. Finally, both cold and reacting LES data are analyzed for uncertainty in the solution using two quality assessment techniques: two-grid estimator by Celik, and model and grid variation by Klein. For both approaches, the resolved turbulent kinetic energy is used to estimate the grid quality and error assessment. The quality assessment reveals that the cold flows are well resolved even on the intermediate mesh, while for the reacting flows even the fine mesh is locally not sufficient in the flamelet region. The Klein approach estimates that depending on the recirculation region in cold scenario both numerical and model errors rise near the bluff-body region, while in combusting flows these errors are significant behind the stabilizing point due to preheating of unburned mixture and reaction heat release. The total error mainly depends on the numerical error and the influence of model error is low for this configuration.     

微小尺度下氢气预混燃烧的实验研究

对内径为1.66mm的不锈钢管燃烧室的氢气预混燃烧实验进行了描述，采用 红外测温仪测量了燃烧室壁面的温度场分布，获得了不同燃烧热功率下的运行界限．在突扩 段内高温回流区的作用下，在带有5mm长突扩段的燃烧室内可以实现完全预混燃 烧，最高运行界限可达1.415．由于较高的进气速度和较大的燃烧室壁面散热，在不带突扩 段的不锈钢管内无法实现完全预混燃烧．结果表明突扩段对微小尺度燃烧具有稳定火焰、拓 宽燃烧运行界限的作用．通过对火焰形状和结构的观察，结合突扩段燃烧流场的分析，合理 解释了燃烧室壁面温度场随过量空气系数的变化规律．     

Dispersion of Entropy Perturbations Transporting through an Industrial Gas Turbine Combustor

In the context of combustion noise and combustion instabilities, the transport of entropy perturbations through highly simplified turbulent flows has received much recent attention. This work performs the first systematic study into the transport of entropy perturbations through a realistic gas turbine combustor flow-field, exhibiting large-scale hydrodynamic flow features in the form of swirl, separation, recirculation zones and vortex cores, these being ubiquitous in real combustor flows. The reacting flow-field is simulated using low Mach number large eddy simulations, with simulations validated by comparison to available experimental data. A generic artificial entropy source, impulsive in time and spatially localized at the flame-front location, is injected. The conservation equation describing entropy transport is simulated, superimposed on the underlying flow-field simulation. It is found that the transport of entropy perturbations is dominated by advection, with both thermal diffusion and viscous production being negligible. It is furthermore found that both the mean flow-field and the large-scale unsteady flow features contribute significantly to advective dispersion — neither can be neglected. The time-variation of entropy perturbation amplitude at combustor exit is well-modelled by a Gaussian profile, whose dispersion exceeds that corresponding to a fully-developed pipe mean flow profile roughly by a factor of three. Finally, despite the attenuation in entropy perturbation amplitude caused by advective dispersion, sufficient entropy perturbation strength is likely to remain at combustor exit for entropy noise to make a meaningful contribution at low frequencies.     

基于液体碳氢燃料的旋转爆轰燃烧特性研究

为了探索液体碳氢燃料参与旋转爆轰所产生的不完全燃烧现象,采用守恒元与求解元方法,开展柱坐标系下的汽油/空气两相旋转爆轰燃烧室三维数值模拟研究,针对燃料喷注压力和反应物当量比对旋转爆轰流场结构及燃烧室性能的影响进行分析。分析结果表明:保持总当量比为1.00,随着燃料喷注压力的上升,燃烧室内燃料不均匀分布增强,产生局部富燃区,燃料在燃烧室未能完全反应,导致燃烧室燃料比冲下降;保持喷注压力不变,减小当量比,在贫燃工况下依然存在局部富燃区,导致燃烧室内出现不完全燃烧现象,降低燃烧室比冲性能。由此可知,反应物喷注方案对气液两相旋转爆轰的不完全燃烧有显著影响。     

Experimental Study on Stabilization of Lifted Swirl Flames in a Model GT Combustor

This study reports on experimental investigations on isothermal and reacting swirled non-premixed flows under varying pressure conditions. In this configuration, a central high speed fuel jet was surrounded by a heated swirling air flow. For the reacting case natural gas served as fuel whereas for isothermal conditions fuel was replaced by a mixture of helium and air to achieve Reynolds-similarity. The optically accessible combustor allowed for application of laser diagnostics. Here we report on Laser Doppler Anemometry and planar laser-induced fluorescence (PLIF) experiments used to characterize the flow field and visualize selected scalars, respectively. Acetone served as a fluorescence marker for mixture fraction investigations. The hydroxyl radical was used to provide general features of the reaction zone such as flame shape and mean stabilization. To expose the influence of pressure on the flame structure three different operating points were investigated varying the combustor pressure between 2 and 6 bar while the inflow bulk velocities remained the same. Striking features of the present configuration are a detached flame, multiple recirculation zones, and complex coherent flow structures.     

Confinement-Induced Instabilities in a Jet-Stabilized Gas Turbine Model Combustor

Self-sustained jet flapping is observed in a confined, premixed and preheated methane-air turbulent flame, generated in a single-nozzle jet-stabilized gas turbine model combustor designed based on the FLOX ^? concept. The flapping frequency and its complex motion within the confinement of the combustor are characterized in detail using proper orthogonal decomposition (POD) of the flow fields measured by particle imaging velocimetry (PIV). The influence of jet flapping on combustion stability is examined using simultaneous PIV/OH chemiluminescence imaging and PIV/planar laser-induced fluorescence of OH radicals (OH PLIF) at 5 kHz repetition rate. By influencing the size and location of the recirculation zones, jet flapping modifies the flame shape and flame lift-off height. It also controls the amount of hot gas entrainment into the recirculation zones. In extreme cases, jet flapping is found to cause temporary local extinction of the flame, due to jet impingement on the combustor wall and partial blockage of burned gas entrainment. The flame is only able to recover after the jet detaches from the wall and reopens the back flow channel. The results suggest that jet flapping could play a key role in the stabilization mechanisms in similar jet-stabilized combustors.     

Measurement and simulation of the two-phase velocity correlation in sudden-expansion gas-particle flows

In this paper the present authors measured the gas-particle two-phase velocity correlation in sudden expansion gas-particle flows with a phase Doppler particle anemometer(PDPA) and simulated the system behavior by using both a Reynolds-averaged Navier-Stokes(RANS) model and a large-eddy simulation(LES).The results of the measurements yield the axial and radial time-averaged velocities as well as the fluctuation velocities of gas and three particle-size groups(30 μ m,50 μ m,and 95 μ m) and the gas-particle velocity correlation for 30 μ m and 50 μ m particles.From the measurements,theoretical analysis,and simulation,it is found that the two-phase velocity correlation of sudden-expansion flows,like that of jet flows,is less than the gas and particle Reynolds stresses.What distinguishes the two-phase velocity correlations of sudden-expansion flow from those of jet and channel flows is the absence of a clear relationship between the two-phase velocity correlation and particle size in sudden-expansion flows.The measurements,theoretical analysis,and numerical simulation all lead to the above-stated conclusions.Quantitatively,the results of the LES are better than those of the RANS model.     

The Effect of Partial Premixing and Heat Loss on the Reacting Flow Field Prediction of a Swirl Stabilized Gas Turbine Model Combustor

This work addresses the prediction of the reacting flow field in a swirl stabilized gas turbine model combustor using large-eddy simulation. The modeling of the combustion chemistry is based on laminar premixed flamelets and the effect of turbulence-chemistry interaction is considered by a presumed shape probability density function. The prediction capabilities of the presented combustion model for perfectly premixed and partially premixed conditions are demonstrated. The effect of partial premixing for the prediction of the reacting flow field is assessed by comparison of a perfectly premixed and partially premixed simulation. Even though significant mixture fraction fluctuations are observed, only small impact of the non-perfect premixing is found on the flow field and flame dynamics. Subsequently, the effect of heat loss to the walls is assessed assuming perfectly premixing. The adiabatic baseline case is compared to heat loss simulations with adiabatic and non-adiabatic chemistry tabulation. The results highlight the importance of considering the effect of heat loss on the chemical kinetics for an accurate prediction of the flow features. Both heat loss simulations significantly improve the temperature prediction, but the non-adiabatic chemistry tabulation is required to accurately capture the chemical composition in the reacting layers.     

Effects of operating pressure on flame oscillation and emission characteristics in a partially premixed swirl combustor

The influence of varying combustor pressure on flame oscillation and emission characteristics in the partially premixed turbulent flame were investigated. In order to investigate combustion characteristics in the partially premixed turbulent flame, the combustor pressure was controlled in the range of −30 to 30 kPa for each equivalence ratio (Φ = 0.8-1.2). The r.m.s. of the pressure fluctuations increased with decreasing combustor pressure for the lean condition. The combustor pressure had a sizeable influence on combustion oscillation, whose dominant frequency varied with the combustor pressure. Combustion instabilities could be controlled by increasing the turbulent intensity of the unburned mixture under the lean condition. An unstable flame was caused by incomplete combustion; hence, EICO greatly increased. Furthermore, EINO_x simply reduced with decreasing combustor pressure at a rate of 0.035 g/10 kPa. The possibility of combustion control on the combusting mode and exhaust gas emission was demonstrated.