SIMULATION OF GASEOUS COMBUSTION AND HEAT TRANSFER IN A VORTEX COMBUSTOR

This article presents the numerical simulation of gaseous combustion and heat transfer in a novel vortex combustor (VC) recently developed for commercial heating applications. A new algebraic Reynolds stress model for strongly swirling turbulence, the eddy break up model for turbulent combustion, and the four-heat-flux model for thermal radiation were employed in the present calculations. Different thermal boundary conditions were specified on the combustor walls. The calculations were conducted on a 22-cm ID VC firing gaseous fuel at 58-, 41-, and 38-kW thermal inputs. Detailed effects of firing rate on the gas axial and tangential velocities, streamlines, temperature, species mass concentrations and wall heat removal are presented, from which the distinct features of the VC flow, combustion, and heat transfer are delineated.     

A NEW ALGEBRAIC MASS FLUX MODEL FOR SIMULATING TURBULENT MIXING IN SWIRLING FLOW

A new algebraic turbulent mass flux model (AFM), which properly accounts for swirl–turbulence interactions, is proposed and formulated in the present article. A simplified explicit form of this new AFM is derived. By incorporating this model with the new algebraic Reynolds stress model proposed previously, a new algebraic stress/flux model (ASM/AFM) is obtained. The new ASM/AFM is applied to the simulation of swirling turbulent flow and mixing in a combustor with helium/air jet and swirling air stream. The calculated gas axial and tangential velocities, gas axial and tangential fluctuating velocities, and helium concentration are in agreement with the measured test data. Comparing with the k–? model for turbulent mass transport, the helium concentration predicted by the new ASM/AFM is much improved.     

Role of entrainment in the stabilisation of jet-in-hot-coflow flames

The aim of the research on the Delft jet-in-hot coflow (DJHC) burner is to gain understanding in the interplay of turbulence and chemistry in conditions as encountered in devices operating in flameless combustion mode, and to test the validity of numerical models when applied to these flameless combustion conditions. Datasets on velocities, temperatures and qualitative OH data of several Dutch natural gas flames in the DJHC burner have been obtained and are discussed in this paper. It was found that the mean velocity and turbulent stresses are not significantly affected by the chemical reactions, which is in line with the very moderate increase of mean temperatures in the flames. Even at heights where flame structures are present, peak temperatures do not always approach the adiabatic flame temperature. With both flame luminescence and OH-PLIF measurements, it is seen that chemical reactions begin to occur at a lower location when the jet velocity (and thereby the jet Reynolds number) is increased. By analysing the velocity and temperature data in the near-nozzle region, the entrainment of coflow fluid into the turbulent jet has been quantified. The increased entrainment of the higher Reynolds number jet, in combination with the positive temperature gradient in radial direction in the near field of the jet, is shown to be responsible for the decrease of the height where reactions start to occur.     

旋流燃烧室内甲烷湍流燃烧的数值模拟

为合理考虑湍流-复杂化学反应的相互作用,建立了甲烷湍流四步反应的温度脉动简化概率密度函数(PDF)模型.应用该模型对TECFLAM燃烧室内的甲烷湍流旋流燃烧进行了数值模拟,得到了与实验相符合的气体轴向、径向与切向速度、温度、温度脉动均方根值及甲烷、氧气、二氧化碳与水蒸气质量分数分布.得到的一氧化碳和氢气质量分数分布与实验基本符合.     

旋流燃烧室内煤粉多相湍流流动与燃烧的数值模拟

提出了考虑湍流-颗粒反应相互作用的颗粒随机轨道模型,以此为基础建立煤粉燃烧综合理论模型并应用于旋流燃烧室内煤粉多相湍流流动与燃烧的数值模拟.模拟结果给出了气相温度场、速度场与温度脉动均方根值分布、颗粒相温度场、速度场与表观密度场以及颗粒瞬时温度与质量随时间的变化.研究表明,考虑湍流-颗粒反应相互作用对气相与颗粒相温度场的模拟结果有一定的影响,使气相温度分布与实验数据更为接近.     

Studies of the effect of a coal concentrator on NO formation in swirling coal combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas–particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas–particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.     

Large-eddy simulation of two-phase spray combustion for gas turbine combustors

In three-dimensional arbitrary curvilinear coordinates, an Eulerian–Lagrangian formulation is applied to large-eddy simulation (LES) of instantaneous gas–liquid two-phase turbulent combustion flows in gas turbine combustors. Three dimensional block-structured grids are generated by zone method and solving a system of elliptic partial differential equations. The k-equation sub-grid scale model is used to simulate the sub-grid eddy viscosity and the EBU combustion sub-grid scale model is employed to predict the chemical reaction rate. The gas-phase governing equations are solved with SIMPLE algorithm and hybrid scheme in non-staggered grid system. A stochastic separated flow formulation is used to track the droplet trajectories velocities, size and temperature history by Lagrangian equations of motion and thermal balance. Multi-zone coupling method is employed to transport data between interfaces. The influences of two different primary hole positions and three different fuel–air ratios on turbulent two-phase reacting flows are calculated. Predictions are in reasonable agreement with the measured velocity using PIV system and temperature, species concentration measurements at the exit. It is shown that the present approach may be used to study spray combustion flow fields for guiding the design of advanced gas turbine combustors.     

Effects of strain rate,turbulence, reactant stoichiometry and heat losses on the interaction of turbulent premixed flames with stoichiometric counterflowing combustion products

Intense strain, turbulence, heat transfer, and mixing with combustion products can affect premixed flames in practical combustion devices. These effects are systematically studied in turbulent premixed CH₄/N₂/O₂ flames using a reactant versus product counterflow system and independently varying bulk strain rate, turbulent Reynolds number, equivalence ratio of the reactant mixture, and temperature of the stoichiometric counterflowing combustion products. The flow field and the turbulent flames are investigated using particle image velocimetry (PIV) measurements and laser-induced fluorescence (LIF) imaging of OH. The OH-LIF images are used to identify the interface between the counterflowing streams, referred to here as the gas mixing layer interface (GMLI). The flame response for different flow conditions is compared in terms of the probability of localized extinction along the GMLI, the turbulent flame brush thickness, and flame position relative to the GMLI, by using an OH-LIF-based progress variable. The probability of localized extinction at the GMLI increases as the separation between the turbulent flame brush and the GMLI decreases. Flame fronts in the vicinity of the GMLI are more likely to extinguish as a result of heat losses, dilution of the reaction zone by the product stream, and large local strain rates. A higher probability of localized extinction at the GMLI is induced by either a larger bulk strain rate or a slower flame speed. As the turbulent Reynolds number increases, the corresponding increase in turbulent flame brush thickness enhances the interactions of the flame fronts with the GMLI. Heat losses are substantially less significant for cases in which the turbulent flame brush is sufficiently separated from the GMLI. For flames in close proximity to the GMLI, the effects of the product stream on the flame front differ for lean and rich reactant mixtures. These disparities are attributed in part to differences in the ignitibility of the reactant mixtures by the hot product stream.     

SIMULATION OF SWIRLING TURBULENT FLOWS AND HEAT TRANSFER IN AN ANNULAR DUCT

The article presents a numerical simulation of swirling turbulent flows and heat transfer in an annular duct. The time-averaged governing equations are solved, which are closed by a new algebraic Reynolds stress model (ASM). The simulation is performed under different flow conditions. The calculated results of gas axial and tangential velocities, turbulent kinetic energy, temperature, and local heat transfer coefficients on the inner and outer walls of the annulus are provided. They illustrate the effect of swirl number, inlet axial velocity, and ratio of inner to outer radius on the mean flow and turbulence properties, as well as on enhancing heat transfer in the annular duct.     

入流条件对同轴射流旋流燃烧室内湍流流动模拟的影响

应用一种新的代数Reynolds应力模型对同轴射流旋流燃烧室内两股射流为同向旋转和反向旋转条件下的湍流旋流流动进行了数值模拟。为得到合理的流场分布结果，研究了入流条件对同轴射流旋流燃烧室内湍流流动模拟结果的影响。计算中对旋流燃烧室进口处两股旋转射流的轴向与切向速度采用了均匀分布和实验测量分布两种方式来给定。将两种进口速度分布条件下得到的燃烧室内气体轴向与切向速度分布计算结果与实验数据进行了比较。     

The numerical computation of the evaporative cooling of falling water film in turbulent mixed convection inside a vertical tube

An analysis has been developed for studying the evaporative cooling of liquid film falling inside a vertical insulated tube in turbulent gas stream is presented. Heat and mass transfer characteristics in air–water system are mainly considered. A low Reynolds number turbulence model of Launder and Sharma is used to simulate the turbulent gas stream and a modified Van Driest model suggested by Yih and Liu is adopted to simulate the turbulent liquid film. The model predictions are first compared with available experimental data for the purpose of validating the model. Parametric computations were performed to investigate the effects of Reynolds number, inlet liquid temperature and inlet liquid mass flow rate on the liquid film cooling mechanism. Results show that significant liquid cooling results for the system with a higher gas flow Reynolds number Re, a lower liquid flow rate Γ₀ or a higher inlet liquid temperature T_L0.     

Large eddy simulation of a 660 MW utility boiler under variable load conditions

Large eddy simulation (LES) has become a promising tool for pulverized coal combustion with the development of computational fluid dynamics (CFD) technologies in recent years. LES can better capture the unsteady features and turbulent structures of coal jet flame than Reynolds averaged Navier Stokes (RANS). The coal-fired power plants in China are now required to be operated in a wide load range and quickly respond to the electric grid. The boiler performance of variable loads should be evaluated in terms of flow, heat transfer, and combustion processes. In this paper, LES was applied to simulate a 660 MW ultra-supercritical boiler under BMCR (boiler maximum continue rate), 75%THA-100, and 50%THA-100 conditions. The predicted gas velocities agree well with the thermal calculation and the temperature error is less than 130 K. The simulation results show that the operation load has significant effects on the boiler performance. It is also proved that LES can provide guidance for the design and operation of advanced coal-fired boilers.     

湿度对HAT循环燃烧室旋流扩散燃烧特性的影响

为研究湿度对燃烧特性的影响,采用湍流雷诺应力模型和层流小火焰模型,对湿空气透平(HAT)循环燃气轮机带有旋流器的燃烧室内甲烷扩散燃烧过程进行了数值模拟对比了在4种不同空气含湿量(0、100、200、300g／kg(DA))情况下的燃烧室内部温度场、速度场以及NO组分分布的情况,分析了湿度对HAT循环燃烧室扩散燃烧特性的影响结果表明,加湿降低了整个燃烧室的温度,并使其内部温度分布更加均匀;加湿使燃烧室的NO浓度大大降低;加湿减小了回流区长度。     

Numerical modeling of two-phase flows in the presence of dynamic and thermal interactions

The study of turbulent flow in mixtures consisting of gas and particulate matter is a current area of research among modern scientists. The main objective of this area is to prevent the reduction of burnout in flow-through systems. Therefore, the purpose of this research is to develop a mathematical model of turbulent gas flow with solids to explore the structure of the combustion region in a channel with a partially open boundary to determine the interfacial interaction of the injection of a flat heterogeneous flow in a confined area and the influence of regime parameters on the ignition process. For this purpose, a modeling method was used, which involved using the fundamental laws of continuity of flow, conservation of energy, and the quantity of motion and matter, and an experiment was conducted to reconcile the proposed model with theoretical data. As a result of this work, a model based on the averaged Reynolds Navier–Stokes equations and an algorithm for solving a turbulent subsonic two-dimensional flow in a heterogeneous medium based on an Euler–Lagrange representation have been developed. Thus, the effects of the parameters of entrainment, temperature, and component concentration on the combustion range and mixing intensity were determined. The proposed methodology can be applied to solve application problems of various kinds, for example, in the design of thermal power plants. The scope of the mathematical and numerical models developed is sufficiently broad that they are universal.     

Simulation of Swirling Turbulent Flow and Combustion in a Combustor

An algebraic concentration moment (ACM)-PDF turbulent combustion model is proposed. It is formulated by jointly utilizing the explicit algebraic expressions for the second-order-moment of concentration fluctuation and the presumed probability density function (PDF) of gas instantaneous temperature. A set of analytical expressions for the time-averaged temperature relevant quantity is obtained for the closure of the time-averaged reaction rate. The model is applied to the simulation of turbulent flow and combustion in a swirl combustor. The calculated gas velocity, temperature, species concentrations, and turbulent fluctuating velocity are in agreement with the measured data. They are much improved over those obtained by the EBU-Arrhenius model.     

耦合天然气详细反应机理的三维湍流预混火焰结构数值预测

以具有279机理的天然气燃烧为例,分别采用涡耗散概念EDC湍流燃烧模型、修正的涡旋破碎EBU湍流燃烧模型以及基于时均值的Arrhenius关系,对燃烧室内复杂的湍流反应流进行了三维数值模拟,并对预测结果进行了分析．结果表明,EDC模型可以较好地反映湍流化学作用,并且能够较好地描述各基元反应,从而为工程实际复杂燃烧情况下其有害排放、中间物质、自由基和痕迹物质生成机理的研究提供基础。     

Evaluation and modification of gas-particle covariance models by Large Eddy Simulation of a particle-laden turbulent flows over a backward-facing step

Three-dimensional numerical investigation of a low speed particle-laden turbulent flow over a backward-facing step has been carried out. An assumption of incompressibility of the flow is used due to low Mach number of the flow. The gas phase is performed by Large Eddy Simulation (LES) and the particle phase is solved by a Lagrangian particle tracking model. The simulation results such as mean streamwise velocities and fluctuation velocities for the both phase are validated by experimental results performed by Fessler and Eaton (1999) [1]. Reynolds number of the gas phase over the backward-facing step with an expansion ratio of 5:3 is 18,400, based on the maximum inlet velocity and step height. The flow is considered as dilute. Hence a one-way coupling method is applied, in which we only consider the effect of fluid on the particle. Particle–particle collisions are also neglected. The success of simulation in predicting a particle-laden turbulent flow using LES and Lagrangian trajectory model provides a numerical basis for revisiting the gas-particle correlations models. Four second-order closure models for gas-particles covariance are evaluated in the present study. A modified better gas-particle covariance model is proposed in this paper.     

当量比对涡轮叶间燃烧性能影响的数值模拟

为探究涡轮叶间燃烧性能,设计了4种不同当量比的工况,利用 FLUENT 软件的 Realizable k-ε湍流模型、PDF 燃烧模型、DO 辐射模型和离散相模型对燃烧室的流动及燃烧进行数值模拟.结果表明:燃烧室能在广泛的当量比(2.59~0.81)下保持性能稳定,燃烧效率保持在96%以上、总压损失低于2.4%,气体温度提高650,K 左右;降低当量比,能够提高燃烧效率,降低 CO、UHC、NOx 等污染物排放,改善温度分布,但会造成更大的总压损失;最优当量比等于1.00,此时燃烧效率在99.95%以上,总压损失相对低(1.5%),出口径向温度呈抛物线型分布,最适合燃烧室设计.与文献对比发现,选取的工况合理,其结果对涡轮叶间燃烧室设计具有参考价值     

Analysis of the interaction between turbulent combustion and thermal radiation using unsteady coupled LES/DOM simulations

Radiation exchanges must be taken into account to improve the prediction of heat fluxes in turbulent combustion. The strong interaction with turbulence and its role on the formation of polluting species require the study of unsteady coupled calculations using Large Eddy Simulations (LESs) of the turbulent combustion process. Radiation is solved using the Discrete Ordinate Method (DOM) and a global spectral model. A detailed study of the coupling between radiative heat transfer and LES simulation involving a real laboratory flame configuration is presented. First the impact of radiation on the flame structure is discussed: when radiation is taken into account, temperature levels increase in the fresh gas and decrease in the burnt gas, with variations ranging from 100 K to 150 K thus impacting the density of the gas. Coupling DOM and LES allows to analyze radiation effects on flame stability: temperature fluctuations are increased, and a wavelet frequency analysis shows changes in the flow characteristic frequencies. The second part of the study focuses on the Turbulence Radiation Interaction (TRI) using the instantaneous radiative fields on the whole computational domain. TRI correlations are calculated and are discussed along four levels of approximation. The LES study shows that all the TRI correlations are significant and must be taken into account. These correlations are also useful to calculate the TRI correlations in the Reynolds Averaged Navier–Stokes (RANS) approach.     

一维湍流预混火焰的数值模拟

建立湍流燃烧的“双流体”数学模型,用于一维湍流预混稳态火焰的描述,假定燃烧火焰由冷的反应物（预混气体）和热的生成物（燃烧产物）组成,它们既有各自的属性,又相互作用,进行热量,质量和动量的交换,采用Patankar和Spalding的Phoenics计算程序来求解该数学模型,成功地模拟了一维湍流参混火焰的压力场,密度场,速度场。