非预混旋流火焰中部熄火及重新稳燃机理研究

本文利用5 kHz同步PIV/OH-PLIF实验装置,结合大涡模拟(LES)耦合PDF燃烧模型,系统研究燃料中心喷射下非预混旋流火焰中部熄火及再次稳燃机理,对比了两种旋流器出口结构的影响.结果 表明,LES-PDF模型可以准确地捕捉到旋流流场分布及火焰形态,包括中部熄火及火焰重新稳定.相比于直燃道的结构,扩张燃道的存在引导了更宽的回流区,改善了中部熄火现象,两种结构对其下游的主火焰再次稳燃高度影响不大.中央燃气射流带来的高标量耗散率使得热损失变大,从而引起中部熄火.同时进一步促进了CH4与空气的预先混合及部分反应前置物(CH2O),为非预混旋流火焰在下游重新稳燃提供了有利条件.通过对火焰再次稳燃处OH反应项及扩散项的分析,发现反应项占据主导地位,部分预混火焰传播为再次稳燃的主要机理.     

湍流强旋流预混火焰大涡模拟研究

本文采用合成涡方法生成湍流进口条件,对非受限强旋流甲烷/空气贫燃预混火焰进行了大涡模拟研究。模拟结果和实验以及DNS数据符合良好。基于此模拟方法,本文对甲烷火焰掺混氢气和水蒸气,进行了比较研究。研究发现,掺混氢气提高了火焰抗拉伸灭火能力,火焰在外剪切层强湍流作用下形成较厚火焰刷。掺混水蒸气,火焰降温明显,流场轴向速度加速较小,形成较大的当地旋流数。因而,相比于掺混氢气,其回流区更强、更宽。相比于掺混水蒸气,流场轴向速度加速较大的掺氢火焰,其轴向脉动在更强的Kelvin-Helmholtz不稳定性作用下,下游剪切层双峰结构得到了保持并增强。     

串列平板叶栅尾迹速度场的LDA实验研究

绕流叶栅的尾涡脱落是诱发水力机械振动噪声的重要因素。本文以串列布置平板叶栅为研究对象,进行雷诺数Re=5000与10000下的叶栅绕流尾迹速度场的LDA测量实验,分析不同雷诺数下绕叶栅流场速度分布,探究涡脱频率特性。实验结果表明:同一雷诺数下平板尾迹区中心线上速度分布可分为回流区、快速增长区、缓慢增长区三个区域;双平板绕流场,下游平板的存在明显抑制了上游平板尾迹的发展,与单平板模型比回流区流向长度减小;雷诺数从5000增大到10000时,平板尾迹回流区的流向长度变小,但最低流速分布升高;下游平板的存在抑制了上游平板的涡脱,使其频率降低,上、下游平板涡脱频率一致。     

旋风分离器内涡核摆动的特性研究

涡核摆动是旋风分离器流体运动的一种非稳定现象,会对流场造成扰动,导致颗粒返混进而影响分离效率.为了探究涡核摆动的规律和机理,本文对单入口蜗壳式旋风分离器内的流场进行研究,运用PDPA测量流场的基本数据,采用Reynolds应力输运模型对其流动进行数值模拟,并定义旋流不稳定性指数TⅡ来分析流动的不稳定性.结果表明,截面Z/D=0.71上瞬时切向速度波动经Fourier变换后,存在一个集中分布的频率,此频率由涡核摆动导致,而且数值模拟结果与实验测量数据吻合较好。通过对涡核中心的运动轨迹分析,发现在整个旋风分离器空间内部都存在涡核摆动现象;在同一时刻,各个截面的涡核中心偏离几何中心的情况并不相同,涡核区域也不一致;各个截面的涡核中心的连线并不在一个平面或者直线上,而是摆动的。此外,分析TⅡ曲线发现气流从不对称的入口结构进入旋风分离器内部,涡核中心开始偏离几何中心,不稳定性逐渐扩大;随着旋流运动向下,入口结构不对称对旋流的影响变小,不稳定性又逐渐减小;流动最终在Z/D=8之后趋于稳定.     

串列平板叶栅尾迹速度场的LDA实验研究

绕流叶栅的尾涡脱落是诱发水力机械振动噪声的重要因素。本文以串列布置平板叶栅为研究对象,进行雷诺数Re=5000与10000下的叶栅绕流尾迹速度场的LDA测量实验,分析不同雷诺数下绕叶栅流场速度分布,探究涡脱频率特性。实验结果表明:同一雷诺数下平板尾迹区中心线上速度分布可分为回流区、快速增长区、缓慢增长区三个区域;双平板绕流场中下游平板的存在明显抑制了上游平板尾迹的发展,与单平板模型相比回流区流向长度减小;雷诺数从5000增大到10000时,平板尾迹回流区的流向长度变小,但最低流速分布升高;下游平板的存在抑制了上游平板的涡脱,使其频率降低,上、下游平板涡脱频率一致。     

多点喷射燃烧室冷态流场特性研究

首先用可实现κ-ε湍流模型(RKE)和雷诺应力湍流模型(RSM)对双径向旋流杯下游流场进行数值模拟,并将计算结果与实验值对比.结果显示,在大部分区域计算值与实验值比较一致,RKE模型和RSM模型的最大误差分别为5%和3%,RKE模型可以用于旋流杯燃烧室计算.此后,用RKE模型对单头部燃烧室和多点喷射燃烧室的流场进行数值模拟.对于同样的燃烧室衬套结构,多点喷射燃烧室旋流器出口附近每个旋流器都有各自的回流区,并在下游逐渐融合,在主燃孔附近多点喷射燃烧室和单头部燃烧室的回流区结构类似.对于同样的多点喷射燃烧室头部,没有主燃孔时,回流区在原主燃孔位置前结束,有主燃孔时,单个的回流区在主燃孔附近又形成新的回流区.     

基于FGM的旋流预混燃烧室超大涡模拟

针对航空发动机燃烧室等复杂工程中的预混燃烧问题发展高精度、高效的数值预测方法,本研究发展了火焰面生成流型(FGM)详细化学反应建表方法结合超大涡模拟方法(VLES),对工程中的GE LM6000预混旋流燃烧室燃烧开展了高精度数值研究,并与实验结果进行了比较。计算结果表明,VLES-FGM方法可以较准确地预测出旋流预混燃烧室内的流场及温度场分布。为了进一步模拟航空发动机真实的燃烧工况,对原始单头部燃烧室使用周期性边界条件来类比全环燃烧室。计算结果表明,VLES-FGM方法计算得到的周期性燃烧室流场回流区相比较固壁边界燃烧室较小,并且固壁边界燃烧室温度场具有明显的颈部结构,燃烧室下游的高温区分布更为均匀。本文计算结果表明基于FGM燃烧模型的自适应湍流模拟方法VLES对于模拟复杂航空发动机相关的旋流预混燃烧具有很大的应用潜力。     

高温横流燃气喷雾掺混的数值研究

高温高压条件下受限空间内旋流喷雾与横向气流的掺混规律尚不清楚,本文对圆形通道内高温燃气与旋流喷雾掺混蒸发过程进行了三维数值模拟,探究流场特征与掺混规律,研究增强流场掺混效果的方法。研究结果表明,四喷嘴周向均匀布置可以使雾化液滴群较均匀的充满掺混空间。相邻喷嘴中间区域液滴群浓度较高,气相温降较大。靠近壁面区域出现多组小尺度的旋涡对结构,小尺度涡结构的发展是促进掺混蒸发过程的主要方式。喷嘴雾化锥角、轴向倾角、切向倾角影响雾化液滴滞留时间及喷雾轴向贯穿深度,进而对掺混效果有较大影响,选择适中的喷嘴雾化角及入射角有利于流场掺混均匀。     

燃气透平进口旋流对热斑迁移及动叶热负荷影响的研究

燃烧室出口旋流及热斑对下游透平叶栅气热特性产生显著影响,本文针对该现象,研究受不同方向和不同进口时序位置的旋流影响下的燃气透平内热斑迁移规律及动叶热负荷特性。结果表明:叶片压力面侧流速较低,流动方向较容易被改变,反向旋流与通道涡的联合作用使热斑核心区向上端壁迁移,而正向旋流作用相反,但两种旋流均使热斑核心区径向范围减小,周向范围增大;旋转的动叶通道使进入的残余旋流因掺混而迅速削弱,对热斑高温流体的径向驱动能力减弱;当热斑与旋流正对流道中心时,并未受到静叶前缘的切割,形态保持相对稳定,且核心区的径向范围较大;进口旋流的存在削弱了进口时序位置的影响,旋流方向对热斑迁移及动叶热负荷的作用更为显著。     

基于Omega涡识别方法的诱导轮内非定常空化流动分析

为精确捕捉诱导轮中空化的行为特征并揭示其非定常流动机理,采用大涡模拟(LES)对一凝水泵诱导轮在设计工况、临界空化数下单独进行数值分析。基于Omega涡识别方法,对诱导轮内空化流动的流动结构及空化泡的关系进行研究。结果表明,Omega涡识别技术可以很好反映流动结构(剪切流与旋转流)与空化泡的关系：设计流量、临界空化数下,诱导轮叶片入口附近空化泡主要集中在轮缘附近的回流区内,回流区内部流动主要以剪切流的形式存在;而在流道中的回流区内流动主要以旋转运动为主,且空化泡主要集中在回流区与主流的剪切层区。     

Flow and mixing fields of turbulent bluff-body jets and flames

The mean structure of turbulent bluff-body jets and flames is presented. Measurements of the flow and mixing fields are compared with predictions made using standard turbulence models. It is found that two vortices exist in the recirculation zone; an outer vortex close to the air coflow and an inner vortex between the outer vortex and the jet. The inner vortex is found to shift downstream with increasing jet momentum flux relative to the coflow momentum flux and gradually loses its circulation pattern. The momentum flux ratio of the jet to the coflow in isothermal flows is found to be the only scaling parameter for the flow field structure. Three mixing layers are identified in the recirculation zone. Numerical simulations using the standard k-? and Reynolds stress turbulence models underpredict the length of the recirculation zone. A simple modification to the C₁ constant in the dissipation transport equation fixes this deficiency and gives better predictions of the flow and mixing fields. The mixed-is-burnt combustion model is found to be adequate for simulating the temperature and mixing field in the recirculation zone of the bluff-body flames.     

Vortex shedding process investigation downstream a surface-mounted block

The laminar flow around a surface-mounted block is investigated by visualizations and PIV measurements. Flow topology and, especially, the vortex shedding dynamics are emphasized. The existence of two vortex shedding processes is highlighted by particles visualizations and instantaneous velocity fields analyses. On one hand, a dominant swirling mechanism with vortex matching process and a symmetrical topology sets up and, on the other hand, a non periodical evolution with a dissymmetric topology exist. In order to inquire about those processes, double velocity correlation functions are calculated and from those new data the space and time evolutions of the vortices are detailed.     

Study of vortex core precession in combustion chambers

The article presents the results of experimental investigation of swirling flow of lean propane/air flame in a model combustion chamber at atmospheric pressure. To study the unsteady turbulent flow, the particle image velocimetry technique was used. It was concluded that dynamics of high swirl flows with and without combustion was determined by a global helical mode, complying with a precessing double-spiral coherent vortex structure. The studied low swirl flame had similar size and stability characteristics, but amplitude of the coherent helical structure substantially oscillated in time. The oscillations were associated with intermittently appearing central recirculation zone that was absent in the nonreacting flow. It is expected that the low swirl flow without the permanent central recirculation zone should be more sensitive to an external active control. In particular, this result may be useful for suppression of thermoacoustic resonance in combustion chambers.     

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.     

旋转回流的一种多重交错网格计算方法

针对圆柱坐标系下原始变量Navier Stokes方程,在有限控制容积法和压力修正的基础上,引入多重交错网格算法及非线性方程的FAS全近似格式,并对封闭圆柱空腔内的旋转流动进行数值模拟.     

Helical modes in low- and high-swirl jets measured by tomographic PIV

We report on a parallel study on properties of large-scale vortical structures in low- and high-swirl turbulent jets by means of the time-resolved tomographic particle image velocimetry technique. The high-swirl jet flow is featured by a well-established bubble-type vortex breakdown with a central recirculation zone. In the low-swirl flow, the mean axial velocity, while intermittently acquiring negative values, remains positive in the mean but with a local velocity defect immediately downstream from the nozzle exit, followed by a spiralling vortex core system and its eventual breakdown. Measurements of the 3D velocity fields allowed direct analysis of the azimuthal/helical modes via Fourier transform over the azimuthal angle and proper orthogonal decomposition (POD) analysis in the Fourier space. A precessing vortex core is detected for both swirl cases, whereas the POD analysis showed that the one originating in the bubble-type vortex breakdown is much more energetic and easier to detect.     

Large eddy simulation of bluff-body stabilized swirling non-premixed flames

Large eddy simulations (LES) using a subgrid mixing and combustion model are carried out to study two bluff-body stabilized swirling non-premixed flames (SM1 and SMA2). The similarities and differences between the two flames are highlighted and discussed. Flow features, such as, the recirculation zone (RZ) size and the flame structure are captured accurately in both cases. The SM1 flame shows a toroidal RZ just behind the bluff body and a vortex breakdown bubble (VBB) downstream. In addition, a highly rotational non-recirculating region in-between the RZ and VBB is observed as well. On the other hand, the SMA2 shows a single elongated recirculation zone downstream the bluff body. Flame necking is observed downstream the bluff body for the SM1 flame but not for the SMA2 flame. The time-averaged velocity and temperature comparison also shows reasonable agreement. The study shows that the sensitivity of the flame structure to inflow conditions can be captured in the present LES without requiring any model changes.     

Vortex identification in the wake of a model wind turbine array

Vortex identification techniques are used to analyse the wake of a 4 × 3 array of model wind turbines. The Q-criterion, Δ-criterion, and λ₂-criterion are applied to particle image velocimetry data gathered fore and aft of the last row centerline turbine. The Q-criterion and λ₂-criterion provide a clear indication of regions where vortical activity exists, while the Δ-criterion is not successful. Galilean decomposition, Reynolds decomposition, vorticity, and swirling strength are used to further understand the location and behaviour of the vortices. The techniques identify and display the high-magnitude vortices in high-shear zones resulting from the blade tips. Using Galilean and Reynolds decomposition, swirling motions are shown encapsuling vortex regions in agreement with the identification criteria. The Galilean decompositions selected are 20% and 50% of a convective velocity of 7 m/s. As the vortices convect downstream, the strength of the vortices decreases in magnitude, particularly in the far wake of the array, to approximately 25% of those present in the near wake. A high level of vortex activity is visualised as a result of the top tip of the wind turbine blade -- the location where the highest vertical entrainment is present. Analysing the full frame set, the Q-criterion, λ₂-criterion, and swirling strength prove comparable, while the Δ-criterion under-performs in regions of high turbulence activity, namely in the back of the turbine. Entraining flow into the turbine canopy interacting with the turbine generates high-magnitude vortices concentrated at the blade tips. The count of vortices decreases when moving from the top tip down to the wall, as well as their strength for each Galilean technique when a non-zero threshold is applied. Vortex sizes in the near wake are found comparable to turbine blade, hub, and mast dimensions. In the far wake, the resulting size of the vortices is approximately 30% of those in the near wake. These vortices increase in velocity as they convect downstream, following the mean velocity behaviour. The lowest magnitude vortices reside at the hub height in the near-wake region, where they convect at nearly half the speed of those at the blade tips.     

不连续双斜向内肋管纵向涡流动显示

利用染色线流动显示方法,显示并分析了不同雷诺数下不连续双斜向内肋管管内流场。流动显示实验结果表明不连续双斜向内肋管管内存在较强的贴近壁面的纵向涡,纵向涡沿流向可持续约一倍管径距离,并且纵向涡的强度随雷诺数的增加而增强。     

PSP measurement of flow fields in a supersonic combustor with a backward-facing step

The mixing fields within a SCRAM-jet combustion chamber are visualized using pressuresensitive paint (PSP) as an oxygen sensor. The experiments are performed in a small supersonic wind tunnel at the National Aerospace Laboratory — Kakuda Research Center (NAL-KRC). The main stream Mach number is 2.4, and the dynamic pressure ratios between the injected gas and the main flow are 0.3, 0.7, 1.1 and 1.5. Three fuel injection nozzles are used; oxygen is injected from the central nozzle and air from the two nozzles at either side. The spread of the injected gas is measured to observe the effects of placing the nozzles in different positions. The results show that the jet has its own independent flow structure, and that little mixing of gases occurs between the flow structures created by each nozzle. When the injection dynamic pressure ratio is increased, the oxygen fraction rises in the recirculation zone and falls in the separation zone downstream of the injection.