旋流燃烧室内同向和反向旋转射流湍流流动的数值模拟

针对发展高效低污染旋流燃烧技术的需要,对旋流燃烧室内两股同轴旋转射流相互作用的湍流旋流流动进行了数值模拟.计算中采用了一种新的代数Reynolds应力模型和QUICK离散格式.在两股射流同向旋转和反向旋转两种条件下,将模拟得到的燃烧室内湍流旋流流动的时均气体轴向速度场、切向速度场和静压场与实验数据进行了比较.     

旋流燃烧室内湍流燃烧的PDA实验研究

本文建立了采用分级进风的旋流燃烧室实验台，实现了用三维激光粒子动态分析仪（ＰＤＡ）测最湍流旋流燃烧的热态瞬时速度场。在一次风旋流数为０．６９的工况下，衍到了燃烧室内气体时均轴向与切向速度和轴向与切向脉动速度均方根值的分布。     

Stability characteristics and flowfields of turbulent non-premixed swirling flames

A simple, yet representative, burner geometry is used for the investigation of highly swirling turbulent unconfined, non-premixed, flames of natural gas. The burner configuration comprises a ceramic faced bluff-body with a central fuel jet. The bluff-body is surrounded by an annulus that delivers a swirling primary flow of air. The entire burner assembly is housed in a wind tunnel providing a secondary co-flowing stream of air. This hybrid bluff-body/swirl burner configuration stabilizes complex turbulent flames not unlike those found in practical combustors, yet is amenable to modelling because of its well-defined boundary conditions. Full stability characteristics including blow-off limits and comprehensive maps of flame shapes are presented for swirling flames of three different fuel mixtures: compressed natural gas (CNG), CNG–air (1:2 by volume) and CNG–H₂ (1:1 by volume).

It is found that with increased fuel flow, flame blow-off mode may change with swirl number, S_g. At low swirl, the flame remains stable at the base but blows off in the neck region further downstream. At higher swirl numbers, the flames peel off completely from the burner's base. Swirling CNG–air flames are distinct in that they only undergo base blow-off. In the low range of swirl number, increasing S_g causes limited improvement in the blow-off limits of the flames investigated and (for a few cases) can even lead to some deterioration over a small intermediate range of S_g. It is only above a certain threshold of swirl that significant improvements in blow-off limits appear. Six flames are selected for further detailed flowfield and composition measurements and these differ in the combination of swirl number, primary axial velocity through the annulus, U_s, and bulk fuel jet velocity, U_j. Only velocity field measurements are presented in this paper. A number of flow features are resolved in these flames, which resemble those already associated with non-reacting swirling flows of equivalent swirl obtained with the present burner configuration. Additionally, asymmetric flowfields inherent to some flames are revealed where the fluidic centreline of the flow (defined in the two-dimensional (U–W velocity pair) velocity field by the ?ω? = 0 tangential velocity contour), meanders strongly on either side of the geometric centreline downstream by about one bluff-body diameter. Flow structures revealed by the velocity data are correlated to flame shapes to yield a better understanding of how the velocity field influences the flames physical characteristics.     

Simulation of spray combustion in a lean-direct injection combustor

Large-eddy simulation (LES) of a liquid-fueled lean-direct injection (LDI) combustor is carried out by resolving the entire inlet flow path through the swirl vanes and the combustor. A localized dynamic subgrid closure is combined with a subgrid mixing and combustion model so that no adjustable parameters are required for both non-reacting and reacting LES. Time-averaged velocity predictions compare well with the measured data. The unsteady flow features that play a major role in spray dispersion, fuel–air mixing and flame stabilization are identified from the simulation data. It is shown that the vortex breakdown bubble (VBB) is smaller with more intense reverse flow when there is heat release. The swirling shear layer plays a major role in spray dispersion and the VBB provides an efficient flameholding mechanism to stabilize the flame.     

LES-CMC of high-altitude relight in an RQL aeronautical combustor

Large-Eddy Simulations with the Conditional Moment Closure sub-grid combustion model and detailed chemistry for kerosene were performed for the ignition process in an Rich-Quench-Lean aviation gas turbine combustor at high-altitude conditions. The simulations used realistic boundary conditions for the flow inlet and spray droplet size distributions and velocity. Due to the large droplets, the Central Recirculation Zone (CRZ) is filled with fuel, mostly in liquid form. The first phase of the ignition process is critical and the results show that the spark kernel must provide enough energy to evaporate the spray and pyrolyse the fuel for the flame to grow and establish in the corner of the combustor. The second phase is characterised by the flame burning the mixture in the scorner and propagating around the Inner Shear Layer. This phase is also critical, as the flame needs the prevaporised fuel and smaller droplets in the corner to sufficiently increase the temperature and be able to propagate inside the CRZ, filled with liquid fuel and cold air. If this propagation inside the CRZ is achieved, phase three is accomplished and the burner is fully ignited. The simulations demonstrate the particular importance of detailed chemistry and proper boundary conditions for flame ignition simulations in high-altitude relight conditions.     

Combustion characteristics of a dual-mode scramjet combustor with cavity flameholder

Combustion characteristics of a laboratory dual-mode ramjet/scramjet combustor were studied experimentally. The combustor consists of a sonic fuel jet injected into a supersonic crossflow upstream of a wall cavity pilot flame. These fundamental components are contained in many dual-mode combustor designs. Experiments were performed with an isolator entrance Mach number of 2.2. Air stagnation temperatures were varied from 1040 to 1490 K, which correspond to flight Mach numbers of 4.3–5.4. Both pure hydrogen and a mixture of hydrogen and ethylene fuels were used. High speed imaging of the flame luminosity was performed along with measurements of the isolator and combustor wall pressures. For ramjet mode operation, two distinct combustion stabilization locations were found for fuel injection a sufficient distance upstream of the cavity. At low T₀, the combustion was anchored at the leading edge of the cavity by heat release in the cavity shear layer. At high T₀, the combustion was stabilized a short distance downstream of the fuel injection jet in the jet-wake. For an intermediate range of T₀, the reaction zone oscillated between the jet-wake and cavity stabilization locations. Wall pressure measurements showed that cavity stabilized combustion was the steadiest, followed by jet-wake stabilized, and the oscillatory case. For fuel injection close to the cavity, a hybrid stabilization mode was found in which the reaction zone locations for the two stabilization modes overlapped. For this hybrid stabilization, cavity fueling rate was an important factor in the steadiness of the flow field. Scramjet mode combustion was found to only exist in the cavity stabilized location for the conditions studied.     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

Self-oscillations in a jet flow and gaseous flame with strong swirl

Investigation results on unsteady flow dynamics in a gaseous jet flame with strong swirl, vortex breakdown, and precession of a vortex core obtained by panoramic optical methods are presented, as well as the results of theoretical analysis of the fastest growing modes of hydrodynamic instability. Characteristics of the most unstable self-oscillating mode in the initial region of the turbulent strongly swirling propane-air jet burning in the atmospheric air in the form of a lifted flame are determined. Analysis of data by principal component analysis and linear stability analysis revealed that evolution of the dominant self-oscillating mode corresponds to quasi-solid rotation with constant angular velocity of the spatial coherent structure consisting of a jet spiral vortex core and two spiral secondary vortices.     

A swirling jet with vortex breakdown: three-dimensional coherent structures

The paper reports on shape of a three-dimensional coherent structure in a velocity field of a high-swirl turbulent jet with the bubble-type vortex breakdown. A set of the 3D instantenous velocity fields was measured by using the tomographic particle image velocimetry (tomographic PIV) technique and processed by the proper orthogonal decomposition (POD) method. The detected intensive coherent velocity component corresponded to a helical vortex core of the swirling jet and two secondary spiral vortices. The entire coherent structure was rotating around the jet axis in compliance with the direction of the flow swirl. From the 3D data it is concluded that the dynamics of the strsucture can be described by a traveling wave equation: Re[A(y, r)·e ^{i(mθ + ky - ωt)}] with the number of the spiral mode m = +1 for positively defined k and ω.     

Scaling relations for the length of coaxial oxy-flames with and without swirl

An extensive experimental study is carried out to analyze scaling laws for the length of methane oxy-flames stabilized on a coaxial injector. The central methane fuel stream is diluted with N₂, CO₂ or He. The annular air stream is enriched with oxygen and can be impregnated with swirl. Former studies have shown that the stoichiometric mixing length of relatively short flames is controlled by the mixing process taking place in the vicinity of the injector outlet. This property has been used to derive scaling laws at large values of the stoichiometric mixture fraction. It is shown here that the same relation can be extended to methane oxy-flames characterized by small values of the stoichiometric mixture fraction. Flame lengths are determined with OH* chemiluminescence measurements over more than 1000 combinations of momentum ratio, annular swirl level and composition of the inner and outer streams of the coaxial injector. It is found that the lengths of all the flames investigated without swirl collapse on a single line, whose coefficients correspond to within 15% of flame lengths obtained for fuel and oxidizer streams at much larger stoichiometric mixture fractions. This relation is then extended to the case of swirling flames by including the contribution of the tangential velocity in the flow entrainment rate and is found to well reproduce the mixing degree of the two co-axial streams as long as the flow does not exhibit a vortex breakdown bubble. At higher swirl levels, when the flow features a central recirculation region, the flame length is found to also directly depend on the oxygen enrichment in the oxidizer stream.     

空心旋转液体射流初始阶段运动规律的研究

应用质量守恒定律和动量守恒定律,建立了描述空心旋转液体射流初始阶段运动的非线 性常微分方程组;该方程组可以用数值方法方便地求解。理论计算结果与实验拍摄到的射流照片吻 合很好。本结果表示射流受挑动失稳破碎成液滴前的最基本运动状态,是进一步从理论上研究空心 旋转射流破碎雾化机理的基础。     

旋转射流冲击换热液晶显示实验研究

采用热色液晶测温技术对以二氧化碳为工质的稳态射流冲击换热和管内插入扭转带方式的旋转射流冲击换热进行了实验研究。与普通射流相比,旋转射流导致驻点附近区域的换热趋于均匀化。其换热系数在大于某一半径之后高于普通射流,但在驻点附近相对较低。旋转射流对换热的此种影响随雷诺数的增大而减弱。     

旋流非预混燃烧火焰流动特性的实验研究

利用粒子图像速度场测量技术(PIV)对不同工况下的旋流非预混燃烧流场进行了测量,考察不同燃空速度比下旋流火焰的流动特性.结果表明,轴向截面上径向平均速度流场以燃烧器轴线呈中心对称,轴向平均速度、轴向脉动速度和径向脉动速度沿燃烧器轴线成对称分布,且轴向平均速度和轴向脉动速度的最大值出现在轴线处.随燃空速度比的增大,轴向平均速度和脉动速度增大,随着与燃烧器表面的距离增加,流场截面上轴向平均速度和脉动速度差异不断减小.     

影响空心旋转液体射流初始阶段运动的无量纲参数的分析

本文以描述空心旋转时流初始阶段运动规律的理论模型[1]为依据,重点讨论了影响射流运动的无量纲参数：韦伯数,液膜内外压力差,液膜初始内外环半径比,初始旋流度;初始径向速度等。从而进一步加深了对旋流喷嘴所产生的空心旋转射流初始阶段运动及发展规律的认识。     

Effects of fuel properties on self-ignition and flame-holding performances in SCRAM jet combustor for n-alkane fuels

Combustion characteristics of liquid hydrocarbon fuels are studied in a model combustor of SCRAM jet engines. The Mach number and total pressure of main flow in the combustor are 2.0 and 0.38 MPa, respectively, and the total temperature is varied from 1800 to 2400 K. Five kinds of n-alkane fuels such as n-heptane, n-octane, n-decane, n-tridecane and n-hexadecane are employed in experiments. Fuels are injected with a carrier nitrogen gas perpendicular to the mail flow in the combustor and the self-ignition behavior is investigated. The results show that the liquid fuels with lower carbon number have better self-ignition performance. This suggests that physical properties of liquid fuels such as volatility have a dominant effect on the self-ignition. The flame-holding behavior is investigated with the addition of pilot hydrogen to carrier nitrogen gas. The critical equivalence ratio at which the stable combustion keeps after cut-off of the pilot hydrogen is obtained. The relationship between the critical equivalence ratio and carbon number of fuel shows that fuels with the carbon numbers from 8 to 10 have the best flame-holding performance among the tested fuels. These experimental results can be expressed qualitatively by the simplified analysis with the concept of physical and chemical induction times.     

A study in the developing region of square jet

The developing region of a turbulent square jet is investigated using high-resolution particle image velocimetry (PIV). The mean velocity and turbulence stresses are presented in various horizontal planes, along the jet centerline covering the initial region of the jet as well as the transition to the self-similar region. To study the flow structure away from the central plane, velocity measurements in two additional horizontal planes, one located halfway from the jet central plane toward the edge and the other at the edge of the square jet, are also examined. Analysis of the instantaneous velocity fields reveal the presence of an arrow-like feature in the square jet due to the higher instability generated in the jet shear layer compared with a round jet. To elucidate the imprints of the vortex structures present in the jets, a swirling strength-based vortex identification methodology is applied on a large ensemble of instantaneous velocity fields. Statistical analysis of the number of vortex cores, and their size and rotational strength in the measurement plane is undertaken. Vortex population at the edge was found to be very different compared with that in the central plane.     

Effect of swirl on combustion dynamics in a lean-premixed swirl-stabilized combustor

The effect of inlet swirl on the flow development and combustion dynamics in a lean-premixed swirl-stabilized combustor has been numerically investigated using a large-eddy-simulation (LES) technique along with a level-set flamelet library approach. Results indicate that when the inlet swirl number exceeds a critical value, a vortex-breakdown-induced central toroidal recirculation zone is established in the downstream region. As the swirl number increases further, the recirculation zone moves upstream and merges with the wake recirculation zone behind the centerbody. Excessive swirl may cause the central recirculating flow to penetrate into the inlet annulus and lead to the occurrence of flame flashback. A higher swirl number tends to increase the turbulence intensity, and consequently the flame speed. As a result, the flame surface area is reduced. The net heat release, however, remains almost unchanged because of the enhanced flame speed. Transverse acoustic oscillations often prevail under the effects of strong swirling flows, whereas longitudinal modes dominate the wave motions in cases with weak swirl. The ensuing effect on the flow/flame interactions in the chamber is substantial.     

涡旋燃烧炉流动、传热和燃烧特性的研究

本文应用强旋湍流气一固两相流动和煤粉燃烧的数学模型，对新型涡旋燃烧炉内的流动、传热和燃烧过程进行了系统的模拟和分析，得到了与实验相符合的结果。结果表明，涡旋燃烧炉内的湍流空气动力场分布具有强旋、回流和正在发展流的特点。水冷壁总吸热量随燃烧热负荷的增大成比例地增加。煤粉颗粒在炉内的平均停留时间随初始粒径的增大而加长。炉内可实现煤粉的低温、强旋、高效率和高强度燃烧。     

旋流燃烧室内颗粒运动的数值模拟

本文应用流体相湍流脉动速度大小和方向均具有随机性的颗粒相随机轨道模型,对有直流一次风和旋流二次风的旋流燃烧室内的颗粒运动进行了数值模拟。得到的颗粒相轴向总质量流通量、轴向与切向速度分布与实验测量数据相符合,并比 Gosman 颗粒随机轨道模型的模拟结果有一定的改进。     

Counter rotating vortex pair structure in a reacting jet in crossflow

This paper analyzes the time averaged flow structure of a reacting jet in cross flow (RJICF), emphasizing the structure of the counter-rotating vortex pair (CVP) by using simultaneous tomographic particle image velocimetry (TPIV) and hydroxyl radical planar laser induced fluorescence (OH-PLIF). It was performed to determine the extent to which heat release, and the associated effects of gas expansion and baroclinic vorticity production, impact the structure of the CVP. These results show the clear presence of a CVP in the time averaged flow field, whose trajectory lies below the jet centerline on either side of the centerplane. Consistent with other measurements of high momentum flux ratio JICF in nonreacting flows, there is significant asymmetry in strength of the two vortex cores. The strength and structure of the CVP was quantified with vorticity and swirling strength (λ_ci), showing that some regions of the flow with high shear are not necessarily accompanied by large scale bulk flow rotation and vice-versa. The OH PLIF measurement allows for correlation of the flame position with the dominant vortical structures, showing that the leeward flame branch lies slightly above, as well as, in the region between the CVP cores.